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ABSTRACT 


The  behavior  of  an  electron  beam  in  a  buncher  gap  is  investigated 
using  a  graphic  analysis  based  on  the  results  given  by  a  Burroughs  Datatron 
220  digital  computer.  The  velocity  and  current  distribution  emerging  from 
this  gap  are  found  The  graphic  results  are  approached  by  successively 
approximating  the  transit  time  correction  factor. 

The  results  obtained  by  successive  approximation  are  then  used 
as  initial  conditions  in  analyzing  the  behavior  of  the  electron  beam  in  the 
drift  space. 
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I.  INTRODUCTION 


The  ballistic  approach  of  analyzing  an  electron  beam  in  a  velocity- 
modulated  tube  was  first  introduced  by  Webster.  ^  Although  this  theory 
lacks  validity  at  or  after  crossover  because  of  space-charge  forces  and 
does  not  apply  to  electron  beams  of  finite  radius  because  of  fringing  of  the 
space-charge  fields,  it  has  two  advantages:  (1)  It  is  simple  and  provides 
an  insight  into  the  physical  phenomena,  and  thus  serves  as  a  guide  for 
more  complicated  theories.  (2)  It  is  fairly  accurate  at  large-signal 
levels  for  low-perveance  beams. 

The  current  tendency  to  demand  higher  power  levels  from  klystrons 

2 

increases  the  importance  of  ballistic  theory.  In  an  experiment,  Mihran 
showed  that  for  large  signals,  the  electron  beam  shows  ballistic  behavior; 
i.e.  ,  at  large  signal  levels,  the  debunching  effect  of  the  space-charge 
forces  becomes  less  important.  At  the  same  time,  in  finite  beams,  the 
space-charge  forces  act  to  enhance  bunching  at  large  signals  by  debunching 
inner  and  outer  electrons  differently. 

The  purpose  of  this  study  is  to  analyze  the  electron  beam,  both  in 
the  first  gap  and  in  the  drift  space,  using  a  ballistic  approach,  and  to 
express  the  functional  relationships  among  the  various  parameters  in  ana¬ 
lytic  form.  To  simplify  this  analysis,  the  following  assumptions  are  made: 
(1)  beam  of  infinite  cross  section,  (2)  negligible  space  charge,  (3)  gridded 
finite  gaps,  (4)  nonrelativistic  velocities. 
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II.  ANALYSIS  OF  THE  FIRST  GAP 


The  first  or  buncher  gap  modulates  the  velocity  of  an  electron  beam. 
This,  in  turn,  causes  density  modulation  in  the  drift  space.  In  a  finite  gap, 
however,  density  modulation  also  takes  place.  This  factor  has  been 
neglected  in  previous  analyses,  and  including  it  extends  Webster's  analysis. 

Analytical  expressions  of  functional  relationships  will  be  based  on 
an  exact  graphic  method  with  the  assumptions  mentioned.  Thus,  in  this 
chapter,  a  graphic  analysis  is  first  made  and  then  analytical  formulas  are 
obtained  by  successive  approximations. 

A.  GRAPHIC  ANALYSIS  (Computer  Problem)^ 

The  equation  of  motion  of  electrons  in  the  buncher  gap,  with  a 
sinusoidal  input  voltage  and  no  space  charge,  is  a  simple  second-order 
differential  equation, 


eV 

1  •  * 

z  =  .  »'  sin  wt 

ma 


(2.1) 


Using  the  notation  of  Figure  1,  integrating  Equation  (2.  1)  twice, 
and  substituting  boundary  conditions  at  t  =  ^>^d  at  t  =  t^  ,  one  obtains 


.  ^^1 
z  =  V  +  ■  ■  j 

o  mujd 


(cos  wt  -  cos  wt)  , 
3i 


(2.2) 


eVi  \  eVj 

d  =  Iv^  +  — cosut^l  (t^  -  t^)  +  - (sinwt^  -  sinwt^) 

'  /  mw  d 


(2.3) 


Defining  the  d-c  transit  angle  and  depth  of  modulation  gives 


% 
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(2.4) 
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and  normalizing  with  respect  to  d-c  velocity  gives 
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7^  (coswt^  -  coswtj^) 

g 


(2.6) 
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g 


cos  ult 
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(wtj^ 


ojt  )  +  ^  ^  (sin  u)t 
a  20  a 
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sin  wt,  ) 
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(2.7) 


Equations  (2.6)  and  (2.7)  completely  define  the  motion  of  electrons  during 
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Figure  1.  Schematic  Diagram  of  Two-Cavity  Klystron  with  Parallel 
Gridded  Gaps. 
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the  passage  through  the  gap  and  give  the  implicit  relation  between  entrance 
time  and  exit  time  with  gap  transit  angle  and  depth  of  modulation  as  para¬ 
meters  . 

The  numerical  computations  of  Equations  (2.  6)  and  (2.7)  have  been 
carried  out  by  a  Burroughs  Datatron  220  digital  computer,  and  are  pre¬ 
sented  in  a  series  of  graphs  of  normalized  exit  velocity,  •  versus 

exit  time,  ut,  ;  normalized  exit  current,  i,  /I  ,  versus  exit  time,  wt,  ;  and 
b  bo  b 

the  transit  time  correction  factor,  wt.  -  wt  -  0  ,  versus  entrance  time, 

bag 

wt 

a 

Exit  current  is  obtained  by  applying  the  principle  of  conservation 
of  charge  through  the  gap;  it  can  be  determined  from 


o 


dwt 
_ a 

dwt, 

b 


(2.8) 


Numerical  values  of  the  normalized  current  are  found  by  measuring  the 

slope  of  the  curves  of  the  transit  time  correction  factor  versus  entrance 

time  at  each  point.  Figures  2a,  b,  c  show  the  results  of  this  method  for 

three  different  gap  transit  angles  with  depth  of  modulation  as  a  parameter. 

Additional  results  will  be  presented  in  later  sections. 

It  can  be  seen  from  Figures  2a,  b,  and  c  that  the  normalized  exit 

current  gradually  becomes  peaked  as  the  d-c  gap  transit  angle  (6  )  is 

8 

increased,  and  that  there  is  a  phase  difference  between  the  normalized 
exit  current  and  the  normalized  exit  velocity.  Since  deceleration  and 
acceleration  of  electrons  during  the  corresponding  half  cycles  of  the 
input  voltage  counteract  each  other,  there  is  a  limit  to  this  peaking  of 
the  exit  current  as  the  d-c  gap  transit  angle  is  further  increased;  0  =  it 
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is  a  practical  gap  transit  angle  and  it  is  used  in  the  computations  through¬ 
out  the  rest  of  this  study. 

The  wave  forms  of  the  exit  velocity  and  the  exit  current  show  that 
harmonics  are  present  in  both  and  that  transit-time  effects  are  nonlinear 
in  the  buncher  gap. 

B.  FIRST-ORDER  ANALYSIS 

When  the  results  of  graphic  analysis  are  known,  one  can  proceed 
with  the  analytical  study  by  making  approximations  in  the  derivations,  and 
comparing  the  results  thus  obtained  writh  those  of  the  exact  graphic  method. 
The  discrepancies  will  show  the  validity  and  the  range  of  the  analytical 
forms . 

1.  Velocity  Modulation 

One  must  approximate  the  transit  time,  since  it  plays  an  important 
role  in  the  physical  phenomena  within  the  buncher  gap.  Assuming  a  cor¬ 
rection  term. 


cat.  =  cat,  +  0  , 

b  a 


(2.9) 


where 


0  »  0  (1  +  6)  0  6  =  cat.  -  cat  -  0  .  (2.  10) 

g  g  bag 


and  substituting  Equation  (2.  9)  into  Equation  (2.  2),  one  obtains  for  f  > 


,e  1~.0.  /.  0\ 

^  pb  ■  i) 


o  m 


.  0 

^  1  . 

=  V  +  »  V  o  -K - sin 

o  2  o 


! 


(2.11) 
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Expanding  sin  (ut^  -  6/2)  by  trigonometrical  identities  and  neglecting  6 
terms  of  higher  order  than  the  linear  term,  one  can  write  Equation  (2.  11) 
as 


—  =  1  +  ^  pa  sin  ((«t^  -  ^  a6  cos-^  sin  (a>t^^  - 


^pa0g  cos 


(2.12) 


where 


■i 


(2.13) 


Integrating  Equation  (2.  2)  and  substituting  in  it  the  boundary  con¬ 
ditions  at  t  =  t  and  at  t  =  t,  ,  and  Equation  (2.  9),  one  obtains 

2l  D 

V  0  a0v  Qv 

<"‘b  ‘  IS?”  <"‘b  •  • 

g  g 

(2.14) 

Assuming 

cos  Cwtb  -  (1  +  5)11^  cos(u>t,  -  0  )  > 


sin  Qwtj^  -  0g  (^  ■•■  ^)II  ^  sin(«t^  -  0^)  , 


sin  Cwtj^  -  0  6]]  sinut^  , 


and  neglecting  the  a6  term  gives 


6  = 


20 


-j  QA  sin(*>t^  +  B  cos(jjt^J 


(2. 15) 


where 
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(2.  16a) 


A  =  1  -  cos  0  -  0  sin  0  , 

g  g  g 


B  =  sin  0  -  0  cos  0 

g  g  g 


(2. 16b) 


Substitution  of  Equation  (2.  15)  into  Equation  (2.  12)  leads  to 


1  /  ^  \ 

=  1  +  A  +  -j-  sin  jut^  -  +  M  sin  2  (ut^  +qp)  ,  (2.  17) 


where 


^  ®g  •  * 


(2.  18a) 


M  =~\/P^  + 


(2.  18b) 


P  = 


=  <*^'^2®g-«i''®g-®g-^20gSin^6g)  ,  (2.18c) 


-  2, 


^  "  7  (■^)  ®g  ■  ^  co8^0g  -  0^  sin  20  J  , 


(2.  18d) 


2  (f  =  tan’^  ^ 


(2. 18e) 


Equation  (2.  15)  is  plotted  in  Figure  3  for  0^  =  it  .  It  can  be  seen 
that  while  Equation  (2.  15)  approximates  well  the  graphic  method  for  small 
signals,  it  is  a  less  satisfactory  approximation  for  large  signals. 


2.  Density  Modulation 

Substituting  Equation  (2.  15)  into  Equation  (2.  10)  and  rewriting  it, 
one  obtains 

wt^  =  wt^  "  ®g  ~  7^  sinwt^  +  B  cosut^3  •  (2-  19) 
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As  in  the  graphic  method,  by  using  Equation  (2.  8),  one  can  express  the 
normalized  exit  current  as 

i, 

r  =  ^  sinwt^  -  A  coswt^^)  (2.  20) 

o  g 

Equation  (2.  20)  can  also  be  written  as 


y—  =  1  +  N  sin  (wt^  +  i)')  , 

o 


where 


KT 

^  -  IT 

g 


=  Ja^  +  b' 


»  I  -  A 

tan  4(  = 


(2.21) 

(2.  22a) 

(2.  22b) 

(2. 22c) 


From  Figure  3  one  can  anticipate  that  Equation  (2.  21)  will  approach 
the  results  of  the  graphic  method  for  small  signals,  but  the  discrepancy 
will  be  considerable  for  large  signals. 

It  should  be  observed  that  the  first-order  analysis  introduces  a 
second  harmonic  term  into  the  velocity  expression,  while  only  the  funda¬ 
mental  is  involved  in  the  current  expression.  One  can  conclude  that  the 
first-order  analysis  results  in  analytic  forms  which  approximate  those 
of  the  graphic  method  for  the  small-signal  case,  and  that  the  method  of 
approximation  for  the  transit-time  correction  factor  appears  promising 
for  the  large-signal  case.  In  the  following  analysis,  therefore,  this 
approach  will  be  followed. 
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C.  SECOND-ORDER  ANALYSIS 


Extending  the  transit-tinne  correction-factor  approximation  one 
step  further,  assume  that 

e  =  (1  +  6j  +  62)  ,  (2.  23) 

where  is  given  by  Equation  (2.  15).  Substituting  Equation  (2,  23)  into 
Equation  (2.  14),  assuming 


cos  (wt,  -  0)  %  cos(wt,_-  0  -0  6,)+0  6,  sin  (wt,  -  0 
b  bggl'g2  b  g 

sin  (wt,  -  0)  ^  sin(wt,  -  0  -  0  6,)-0  6,  cos  (wt, 

b  'b  g  gl'  g2  'b 

and  neglecting  a62  and  &2  terms,  one  obtains 


-  Vi> 


- 


6-  =  -  6  + Csinwt.  -  sin(wt.  -  6  -  6j)co8  (wt,  -  0  -0  6  )H, 

D&O  60 


20 


(2.  24a) 


and 

wt,  -  wt  -  0  =  Fsin  wt,  -  sin  (wt,  -  0  -0  6,)-0(l  +  6,  )cos(wt,  -  0  -0  6, )”]. 

bag  2^^  b  '  b  g  g  1'  g'  1'  '  b  g  g  1'-* 

(2.  24b) 


Equation  (2.  24b)  is  plotted  in  Figure  3  for  0  =  u  .  It  can  be 

S 

seen  from  this  figure  that  Equation  (2.  24b)  gives  a  fair  approximation  to 
the  gr  aphic  results  for  larger  signal  values. 

1.  Velocity  Modulation 

Substitution  of  Equation  (2.  24a)  into  Equation  (2.  13)  leads,  after 
some  rearrangement,  to 
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Ik  - 1  + W  f 

o  '8/ 


(cosGg-  1) 


-  2  2 
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JTj 


61  L 


0  /  \2 

1  -  cos  9g  +  (2B  -  A  sin  20^  -  B  cos  20^) 

6' 


6' 


cos  cat. 


4 


pos  0  -  A  sin  0  +T 

6  g  2 


i  (irf  C'A  B^)sin  20g-  2AB3]|  s  in  2cat 


^  +  B^)cos  20^+ -B^D^cos  2cat^ 


0  /  v3 

(A  cos  20^  +  B  sin  20^)  sin  3catj. 


0  /  a  /  \4 

■  ■'■^  [^j  sin  20^ 


+  2AB  cos  20  "]sin4wtj^  4(z?-)  Qb^-A^)cos20  -  2AB  sin  20  ]]cos  4cat, 


(2.25) 

Equation  (2.  25)  is  plotted  in  Figure  4  for  0^  =  it  for  comparison 
with  the  graphic  method  and  in  Figures  5  and  6  with  the  d-c  gap  transit 
angle  as  a  parameter.  It  can  be  seen  from  Figure  4  that  Equation  (2.  25) 
approximates  the  graphic  method  fairly  well  up  to  large  signals.  Figures 
5  and  6  show  the  effect  of  the  d-c  gap  transit  angle  on  the  velocity  modu¬ 
lation.  It  can  be  seen  that,  as  expected,  the  modulation  effect  decreases. 
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Figure  4.  Normalized  Exit  Velocity  versus  Exit  Time  ut^^, 

Comparing  Graphic  and  Second-Order  Methods. 


ao 
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Lti. 
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Figure  5.  Normalized  Exit  Velocity  versus  Exit  Time  wt^  for  Different  D-C  Gap 

Transit  Angles  and  Depths  of  Modulation. 


because  of  partial  cancellation,  as  the  d-c  gap  transit  angle  is  increased. 
2.  Density  Modulation 

Differentiating  Equation  (2.  24b)  and  neglecting  third-order  terms, 
one  obtains 


^  ‘  ^ 


=  1 


a  r  / 

TT  "‘“•“'b  -  K  ■  V  ^  ®g  ^  •  ®g  TO 

B  u  •  C 


^^1  \  Z  (  2  \ 

■  ®g®i  ■  V  ■  ®g  n  "  ®g^iTO)^°®  <“‘b  ■  V 


b  g' 

(2.  26) 


Substituting  Equation  (2.  15)  and  its  derivative  with  respect  to 
ut^  and  rearranging  terms  gives 

.2 
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g 


0_sin0  >coswt 


_ -  — ■  ;  '  o_  sin  o 
8  Z  g  g 


1! 


+  0g  (A  cos  0g  +  B  sin  0^)  sin  2«t^  -  1  ®g(A  sin  0^-  B  cos  0^)cos  2ut^ 

'\  8  \  ff/ 


a  \^r  (A^  (B^-A^.0  AB) 

AB  cos  0g  -  ' - ^ - !  0^  cos  0^  +  - -  g  g -  sin  0 


sin  3wt. 


[  ® 

\ir 

g 


(B^  -A^  -  0  AB) 

AB  sin  0^  -  ^  ^  ^  ®g  sin  0^  -  ^  ^  cos  0 


(A^-B^) 


cos  3wtj^ 
(2.27) 
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Equation  (2.  27)  is  plotted  in  Figure  7  for  0^  =  ir  for  comparison 
with  the  graphic  method  and  in  Figures  8  and  9  with  the  d-c  gap  transit 
angle  as  a  parameter.  Again  it  can  be  seen  from  Figure  7  that  Equation 
(2.  27)  approximates  the  graphic  method  fairly  well  up  to  large  signals. 
Figures  8  and  9  show  the  effect  of  the  d-c  gap  transit  angle  on  density 
modulation.  It  can  be  seen  that  as  in  the  graphic  method,  density  modu¬ 
lation  first  increases  with  the  d-c  gap  transit  angle  and  then  decreases. 

It  should  be  observed  that  the  second-order  analysis  introduces 
harmonics  up  to  the  third  harmonic  in  the  current  expression,  and  up  to 
the  fourth  harmonic  in  the  velocity  expression,  and  that  it  extends  the 
approximation  to  larger  signals. 
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GRAPHIC  - GRAPHIC 

SECONO-OROER  I  *  SECOND-ORDER 


Figure  7.  Normalized  Exit  Current  versus  Exit  Time  u)t^  Comparing  Graphic  and 

Second-Order  Methods. 


Normalized  Exit  Cur¬ 
rent  i,  /I  versus  Exit 
D  o 

Time  wt^  for  Different 
D-C  Gap  Transit  Angles 
and  Depths  of  Modulation 


III.  ANALYSIS  OF  THE  DRIFT  SPACE 


The  electron  beam  enters  the  drift  space  with  both  velocity  modu¬ 
lation  and  density  modulation,  and  it  drifts  in  a  field-free  space.  This 
drift  action  produces  further  bunching  of  electrons,  thereby  increasing  the 
harmonic  content  of  the  beam  current.  In  his  analysis,  Webster  neglects 
the  density  modulation  produced  in  the  gap,  since  his  analysis  deals  with 
small  signals.  It  was  shown  in  the  previous  chapter  that  the  density  modu¬ 
lation  in  the  first  gap  was  considerable  for  large  signals.  The  present 
analysis  will  consider  the  density  modulation  in  the  gap.  It  will  include 
only  first-order  terms  and  will  not  consider  space -charge  debunching. 
Although  the  present  analysis  is  not  valid  after  crossover,  it  will  be 
extended  beyond  crossover;  for  qualitatively,  it  anticipates  important 
trends. 

The  time  at  plane  c  (see  Figure  1)  can  be  expressed  as 


Shifting  the  reference  of  time, 

.  e 

“'‘b  ^  “S  ■  z  • 

Dropping  the  prime,  and  neglecting  the  second-harmonic  term  in  Equation 
(2.  17),  one  obtains  from  Equation  (3.  1), 

wt^  -  <j)  =  -  X  sinwt^  ,  (3.3) 


-21- 


where 


<j>  =  0^(1  -  A)  +-^  , 


(3.  4a) 


0  -  ~ 

®o  =  ~  ' 

o 


(3.4b) 


X  s  ^  |1CL  ' 


(3.4c) 


Equation  (3.  3)  shows  the  functional  relationship  between  the  entrance  and 
exit  tinnes  of  the  drift  space. 


A.  DENSITY 

The  principle  of  conservation  of  charge  for  the  drift  space  states 


i.  db>t,  =  i  dwt 
b  b  c  c 


(3.5) 


One  can  express  the  current  at  plane  c  as  a  Fourier  series  given  by 


i  =  A  + 
c  o 


y  A^  cosn  (wt^  -  4»)  +  B^8inn{«t^  -  ^)  , 


(3.6) 


where 


A  =  I 
o  o 


(3..  7a) 


i  =  —  I  i  cos  n  (wt  -  (b)  dut 
n  ir  c  c  ^  c 


(3.  7b) 


B  =  —  I  i  sinn  (wt  •  d>)  dwt 
n  IT  I  c  '  c  c 


(3.7c) 
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If  one  uses  Equations  (2.  20),  (3.  2),  (3.  3)  and  (3.  5),  then  Equations  (3.  7b 
and  c)  become 


A  =  — 
n  IT 


1  +  N  sin  (ut^  +  4»  + 


-^jljcos  n  Qwt^  -  X  sinwt^^dwt^  ,  (3.8a) 


®n  =ir 


1  +  N  sin  <i)t^  +  4^  +  -f  sinn  -  X  sinut^^J  (3* 


Using  Bessel  function  expressions^  for  trigonometric  functions,  integrating, 
and  using  recursion  formulas  for  Bessel  functions,  one  obtains 


oo 


1  +  2 


1  + 


^  sin  J^(nX)  cosn  (wt^  -  ^)  , 

n  =  1 


/  0  \  ^  j’  (nX) 

+  2N  cosi4<  .  sinn  (ojt^  -  ^)  .  (3.9) 


It  can  be  seen  from  Equation  (3.  9)  that  the  beam  current  at  plane  c 
is  rich  in  harmonics  and  that  each  harmonic  can  be  calculated  for  any 
specified  condition  from  the  equation.  Equation  (3-9)  reduces  to  that  of 
Webster  for  small  signals.  It  should  also  be  noted  that  an  additional 
phase  angle,  ,  is  introduced  by  the  two  independent  components  of 
each  harmonic,  given  by 


N  coslv|<  + 


tan 


n 


N  a 

1  +  ^  .inU  +  -f 


(3.10) 


J  (nX) 
n  '  ' 
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Figures  10a,  11a,  and  12a  show  the  first  three  harmonics  and 
Figures  10b,  11b,  and  12b  show  the  phase  angles  of  these  harmonics  as 
a  function  of  the  bunching  parameter  (X),  with  the  depth  of  modulation  as 
a  parameter,  for  0  =  ir  .  The  maximum  amplitudes  of  the  fundamental 

current  for  different  depths  of  modulation  (see  Figure  10a)  are  seen  to  be 
greater  than  the  usual  value  of  1.  16,  which  was  also  predicted  in  a  bal¬ 
listic  analysis  by  Webber^  that  included  space -charge  effects. 

The  first  three  harmonics  are  also  plotted  in  Figure  13  as  functions 
of  the  drift  angle  9^  for  6^  =  ir  .  Qualitatively  Figure  13a  resembles 
the  experimental  curves  of  Mihran;  he  observes  a  saturation  of  the  maxima 
of  the  fundamental  current,  although  he  did  not  explain  the  reason  for  this. 

Although  space-charge  effects  were  neglected  in  this  study,  the 
results  obtained  parallel  those  of  others  within  the  range  of  validity.  A 
comparison  with  Solymar's^  results  is  shown  in  Figure  14,  where  it  was 
assumed  that 


B.  VELOCITY 

The  transit  time  through  the  drift  space  can  be  assumed  to  be 


where 


wtb  -  (wt^  -  4>)  = 


sin  1  (ut 

c 


-  4*) 


(«t^  -  40  1  d  («t^  - 


(3.12) 


(3.13) 
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t*  •  <!> 


re  10.  (a)  Absolute  Value  of  Normalized  Fundamental  Current 

|ij|/I^  versus  Bunching  Parameter  X  ,  (b)  Phase  Angl 
P  ,  versus  Bunching  Parameter  X  . 


MMM  AN*C£ 


Figure  11.  (a)  Absolute  Value  of  Normalized  Second -Harmonic  Current 

|i  I /I  versus  Bunching  Parameter  X  ,  (b)  Phase  Angle 
p  versus  Bunching  Parameter  X  ,  in  Drift  Space. 


(a)  ■UNCHWO  MNAMCTEH  X 


Figure  12.  (a)  Absolute  Value  of  Normalized  Third-Harmonic  Current 

UjI/Iq  versus  Bunching  Parameter  X  ,  (b)  Phase  Angle 
versus  Bunching  Parameter  X  in  the  Drift  Space. 
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DRIFT  ANGLE  •«  (RADIANS) 


Figure  14.  Normalized  Fundamental  Current  ij/I^  versus  Normalized  Drift  Distance  p^z 
for  Different  D-C  Gap  Transit  Angles.  (The  curves  shown  are  those  calculated 
by  Solymar.  The  solid  dots  are  the  points  calculated  from  Equation  (3.9)  of  this 
study. ) 


Integrating  Equation  (3.  13)  by  parts  and  substituting  Ecuation  (3.  3),  one 
obtains 


wt,  =  wt  -  ()>  + 


y- Jj(lX)  Sin  1  (wt^  -  <j>) 


(3.14) 


Equation  (3.  14)  shows  the  functional  relationship  between  the 
entrance  time  and  the  exit  time  of  the  drift  space.  Figure  15  shows  a 
comparison  of  Equation  (3.  3)  and  (3.  14),  where  only  the  first  three 
harmonics  are  included  in  Equation  (3.  14). 

Since  the  drift  space  is  assumed  to  be  field -free,  the  equation  of 
motion  becomes 

z  =  0  .  (3.15) 

Integrating  Equation  (3.  15)  and  applying  the  boundary  conditions  at 
t  =  tjj  ,  substituting  Equation  (2.  17)  with  second  harmonic  terms  neglected 
and  Equations  (3.  2)  and  (3.  14),  one  obtains 


V 

c 

V 

o 


1  +  A  + 


1 

.y  fia  sin 


-  + 


2Jj(lX) 
- 1 - 


sin  1  (ut 
'  c 


(3.16) 


Equation  (3.  16)  is  plotted  in  Figure  16  with  depth  of  modulation  as  a 

parameter  for  9  =  tt  and  0  =  ir/Z  . 

go 
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Figure  15.  Arrival  Phase,  -  L(1  -  A)  ,  versus  Departure  Phase 

for  the  Drift  Space  for  Different  Bunching  Parameters. 


Figure  16.  Normalized  £xit  Velocity  v  /v  versus  Exit  Time  wt 
®  '  c  o  c 

at  End  of  Drift  Space  for  Different  Depths  of  Modulation. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


Although  this  study  was  based  on  a  very  simple  model,  many 
interesting  and  important  results  were  obtained  that  may  serve  as  guides 
for  future  investigations. 

Exact  graphic  analysis  reveals  the  physical  relationship  in  the  first 
gap  and  in  the  drift  space,  but  it  is  tedious  and  time  consuming  and  does 
not  give  the  functional  relationships  among  various  parameters.  Approxi¬ 
mating  the  graphic  method  analytically,  on  the  other  hand,  is  demonstrated 
to  be  valid:  the  second-order  method  both  approximates  the  large-signal 
cases  well  and  also  gives  analytical  formulas  which  reveal  information 
about  the  harmonics. 

The  conclusions  from  this  ballistic  analysis  can  be  grouped  under 
the  first  gap  and  the  drift  space. 

A.  FIRST  GAP 

1.  Both  velocity  and  curre  .t  modulations  are  produced  during  the 
passage  of  the  beam  through  the  gap,  and  the  transit  time  plays  an 
important  role  in  these  phenomena. 

2.  The  velocity  modulation  is  proportional  to  the  depth  of  modu¬ 
lation,  but  is  inversely  proportional  to  the  d-c  gap  transit  angle. 

3.  The  difference  between  the  minimum  velocity  and  the  d-c 
velocity  is  greater  than  the  corresponding  difference  between  the  maximum 
velocity  and  the  d-c  velocity. 

4.  The  current  modulation,  on  the  other  hand,  is  also  proportional 
to  the  depth  of  modulation;  but,  first  it  increases  and  then  it  decreases  as 
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the  d-c  gap  transit  angle  is  increased  continously. 

5.  Both  velocity  and  current  expressions  contain  harmonics  at 
the  exit  of  the  first  gap,  but  this  effect  is  more  pronounced  in  the  current 
expression. 

B.  DRIFT  SPACE 

1.  In  a  field-free  space,  density  modulation  is  produced  in  a 
velocity-modulated  beam  through  drift  action.  The  beam  current  develops 
harmonics,  and  the  amplitudes  of  harmonics  are  proportional  to  the  depth 
of  modulation. 

2.  The  maximum  value  of  the  fundamental  current  is  seen  to  be 
greater  than  the  value  of  1.16  predicted  by  Webster.^  The  point  in  the 
drift  space  at  which  the  maximum  value  of  a  harmonic  occurs  is  inversely 
proportional  to  the  depth  of  modulation.  These  results  have  been  shown, 
at  least  qualitatively,  in  the  work  of  others. 

3.  An  additional  phase  angle  is  introduced  into  the  current  by 
the  two  independent  components  of  each  harmonic. 

4.  The  velocity  becomes  more  nonlinear  at  the  end  of  the  drift 
space;  the  extreme  values  do  not  change  considerably. 

C.  RECOMMENDATION 

This  report  has  dealt  with  the  simplest  case.  Investigation  of  the 
second  gap  and  inclusion  of  the  effects  of  space  charge  and  of  gridless  gaps 
for  finite  beams  would  be  a  logical  continuation  of  this  study. 
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ABSTRACT 


A  general  theoretical  expression  is  derived  for  an  electron  trajectory 
in  an  oscillating  electric  field.  The  limits  upon  the  validity  of  the  analysis 
are  found,  and  application  is  made  to  the  ballistic  theory  of  gridded  klystron 
gaps.  In  the  case  of  the  input  gap,  exact  results  for  exit  current  and  exit 
velocity  are  obtained  for  large  signals,  the  only  restriction  being  that  the 
normalized  modulation  voltage  magnitude  be  less  than  the  normalized  gap 
length.  Large-signal  analytical  expressions  for  exit  current  and  exit 
velocity  are  obtained  for  an  output  gap,  and  results  of  numerical  calculations 
are  presented  graphically  for  several  arbitrary  cases.  It  appears  from 
these  results  that  an  r-f  structure  placed  beyond  the  output  gap  could  be 
used  to  reduce  the  velocities  of  high-velocity  electrons  in  the  beam,  and 
help  reduce  production  of  X-rays  at  the  collector. 


I.  INTRODUCTION 


With  the  invention  of  the  klystron  by  the  Varian  brothers  in  1939 
came  the  first  theoretical  analysis  of  the  device  by  Webster.  ^  His  analy¬ 
sis  used  ballistic  theory,  but  he  limited  his  treatment  with  a  restriction 
to  small  input  signals.  Many  papers  have  recently  been  written  about 
klystron  operation,  but  theoretical  analyses  have  been  almost  invariably 

limited  by  this  same  assumption.  The  space -charge  wave  theory,  proposed 
2  ? 

by  Hahn  and  also  Ramo,  was  also  limited  by  the  small-signal  assump¬ 
tions,  but  since  it  included  the  effects  of  space  charge  in  the  drift  region, 
it  was  presumably  more  accurate  and  has  recently  enjoyed  much  popularity. 

Both  of  these  approaches  dealt  primarily  with  the  behavior  of  the 
beam  in  the  drift  region,  and  relatively  little  was  said  about  the  behavior 
of  the  beam  in  the  two  gap  regions.  Usually,  the  gap  was  assumed  to  be 

gridded,  with  a  very  short  length,  and  it  was  assumed  that  no  power  was 

« 

added  to  the  beam  by  the  input  gap.  The  input  gap  was  supposed  to  pro¬ 
duce  pure  velocity  modulation  with  no  current  modulation,  and  these  con¬ 
ditions  were  used  as  boundary  conditions  on  the  treatments  of  the  drift 
space.  Early  treatments  of  the  input  gap  region  were  similar  to  the 
method  given  here,  in  that  they  considered  the  electron  passing  through 
a  time -varying  electric  field  that  was  constant  in  the  Z  direction,  the 
direction  of  beam  flow.  The  velocity  was  calculated  at  the  end  of  this 
gridded  gap.  The  more  recent  approach  has  been  to  use  kinetic  energy, 
and  assume  that  the  gap  is  so  short  that  the  difference  between  entrance 
and  exit  time  is  negligible.  Thus  the  electric  field  would  be  essentially 
a  constant  in  time  for  each  electron.  If  the  gap  angle  =  wd/2^  is  not 
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small,  this  procedure  is  incorrect,  because  the  kinekic  energy  is  based 
upon  a  time  integral  of  Newton's  second  law  where  the  force  is  not  a 
function  of  time.  This  will  be  illustrated  more  clearly  in  the  exact  bal¬ 
listic  treatment  in  Section  II. 

Modern  klystron  designers  are  rersponding  to  the  demand  for  higher 

and  higher  power  outputs,  and  several  theoretical  analyses  have  treated 

this  problem,, for  in  high  power  operation,  the  small-signal  assumptions  of 

previous  theories  are  invalid.  Perhaps  the  most  fruitful  approach,  al- 

4  5  6 

though  a  tedious  one,  is  that  of  Paschke.  ’  '  Paschke  writes  the  general 
nonlinear  space -charge  wave  equation  for  the  drift  space,  and  solves  it 
by  successive  approximation:  that  is,  he  solves  for  a  linear  solution 
similar  to  that  given  by  Hahn  and  Ramo  and  substitutes  this  into  the 
nonlinear  equation  to  get  a  second -order  solution.  This  process  is  re¬ 
peated  to  get  a  third-order  solution.  His  results  have  been  verified  by 

7  8 

experimental  work  of  Mihran.  ’  Paschke,  however,  uses  the  traditional 
small-signal  assumption  at  the  input  gap,  where  he  applies  the  boundary 
condition  of  pure  velocity  modulation  to  his  drift-space  analysis. 

A  clear  exposition  of  the  problem  faced  in  ballistic  gap  theory  is 
given  by  Beck.  ^  Beck  extends  the  analysis  of  the  gridded  gap  to  the  case 
of  a  gap  length  which  is  nonzero,  but  small.  His  first  approximation 
involves  the  assumption  that  the  transit  angle  of  the  electron  going  across 
the  gap  is  the  d-c  transit  angle  ,  and  subsequently  he  allows  the  actual 
transit  angle  to  deviate  from  this  by  a  small  amount  5  .  By  several 
approximations,  a  fairly  simple  expression  for  6  is  derived.  This  ex¬ 
pression  is  useful,  but  no  limits  upon  its  validity  have  been  established, 
and  whether  it  is  valid  for >extremely  large  signals  is  not  known.  A 
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similar  approach  to  this  problem  has  been  given  by  Blair,  except 
that  Blair  extends  these  approximation  methods  to  higher  orders.  The 
simplest  large -signal  gap  theory  has  undoubtedly  been  that  presented  by 
Warnecke  and  Gurnard.  Their  graphic  methods  are  applicable  to  very 
large  signals  and  give  rather  rapid  results.  There  seems  to  be  no  limit 
to  the  applicability  of  graphic  techniques,  but  they  have  the  disadvantage 
of  only  moderate  accuracy  and  must  be  repeated  for  every  choice  of 
parameters.  The  problem  has  been  very  well  formulated  by  Warnecke 
and  Gurnard,  however,  and  has  been  clearly  restated  by  Kleen.^^  The 
methods  they  give  are  definitely  the  best  available  to  date  for  the  extremely 
large  signal  treatment  of  the  gap  region. 

All  of  the  methods  of  analysis  mentioned  above  with  the  exception 
of  the  graphic  method  share  the  same  problem  of  small-signal  limitation. 

It  would  be  useful  to  obtain  an  analytic  expression  that  is  valid  for  large 
signal  levels.  Almpst  as  important  as  obtaining  the  expression  is  finding 
the  limit  upon  its  validity;  this  requirement  has  been  almost  universally 
neglected. 

Such  an  analytic  expression,  when  found,  could  be  used  to  establish 
rigorous  boundary  conditions  on  the  large -signal  treatments  of  the  drift 
region.  Furthermore,  the  gap  theory  coupled  with  the  drift-space  theory 
can  give  a  fairly  complete  theoretical  description  of  the  operation  of 
modern,  high-power  klystrons.  This  would  include  answers  to  such 
problems  as  (1)  analyzing  the  beam  current  for  harmonic  content  that 
may  produce  spurious  outputs  and  (2)  finding  the  velocity  and  current 
distribution  of  the  electrons  in  the^beam  after  passing  through  the  output 
gap  to  give  information  about  the  production  of  dangerous  high-intensity 
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X-rays  at  the  collector.  In  this  treatment,  application  of  the  theory  is 
made  to  the  second  problem.  Finally,  it  would  be  possible  ultimately 
to  apply  the  theory  for  the  2 -cavity  klystron,  to  the  multicavity  klystron, 
which  is  becoming  so  useful  in  high-power  applications. 

II.  LARGE  SIGNAL  BALLISTIC  THEORY  OF  GRIDDED  GAPS  WITH 
ARBITRARY  INPUT  VELOCITY 

A.  SYMBOLS  USED  IN  THIS  SECTION 

magnitude  of  applied  voltage  _  V 
d-c  beam  voltage 

A  =  small  increment  of  any  quantity  following, 
c  s  eccentricity,  defined  in  Equations  (13),  (32). 

s  mean  anomaly,  defined  in  Equations  (14),  (33). 

6  s  phase  of  applied  voltage  relative  to  time  t  =  0  . 

V  s  summation  index,  equal  to  0,  1, 2,  3,  4  ...  . 
p  =  total  charge  density  at  entrance  plane  a  . 

=  total  charge  density  at  exit  plane  b . 

=  normalized  length  of  the  gap,  equal  to  -j- 
<i)  s  angular  frequency  of  applied  voltage. 

A  s  unit  cross-sectional  area  of  a  beam, 
b^  s  coefficient  in  Fourier  series  expansion, 
d  =  length  of  the  gridded  gap. 
e  =  electronic  charge. 

E  =  electric  field  in  the  gridded  gap. 

g  =  normalized  time -varying  part  of  the  entrance  velocity, 
defined  in  Equation  (30). 

I^  =  total  conduction  current  at  entrance  plane  a. 
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=  total  conduction  current  at  exit  plane  b. 

=  d-c  beam  current  before  any  gap  modulation. 

=  total  conduction  current  at  any  arbitrary  plane  Z  within  the 

gap. 

=  total  conduction  current  density  at  entrance  plane  a. 

=  total  conduction  current  density  at  any  arbitrary  plane  Z 
within  the  gap. 

=  Bessel  function  of  the  first  kind  of  the  order. 

=  Bessel  function  of  the  first  kind  of  the  order. 

=  summation  index,  equal  to  1 , 2,  3,  4,  ...  . 

=  electronic  mass. 

=  arbitrary  ternri  index  in  the  series,  Equation  (41)  only. 

=  number  of  particles  per  unit  volume  at  an  arbitrary  plane  Z 
in  the  gap. 

=  number  of  particles  per  unit  volume  at  the  entrance  plane  a 
of  the  gap. 

=  total  charge  at  plane  a. 
s  total  charge  at  plane  b. 

=  time  that  an  electron  arrives  at  an  arbitrary  plane  Z  within 
the  gap. 

=  time  that  an  electron  passes  entrance  plane  a. 

=  time  that  an  electron  passes  exit  plane  b. 

=  eccentric  anomaly,  defined  in  Equation  (12). 

=  magnitude  of  applied  voltage. 

=  magnitude  of  d-c  beam  voltage. 

=  total  velocity  of  an  electron  at  arbitrary  time  t  within  the 
gridded  gap. 

=  total  velocity  of  an  electron  at  entrance  plane  a,  at  time  t^ . 

=  total  velocity  of  an  electron  at  exit  plane  b. 

=  arbitrary  position  co-ordinate  of  an  electron  within  the  gap. 
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Z  =  total  velocity  of  an  electron  at  an  arbitrary  plane  within  the 
gap.  Dot  indicates  time  derivative. 

=  d>c  beam  velocity,  defined  in  Equation  (29). 

Z  =  acceleration  of  an  electron  at  an  arbitrary  plane  in  the  gap. 

B.  THE  ANALYTICAL  PROBLEM 

The  ballistic  theory  developed  for  electron  beams  neglects  any 
forces  on  an  electron  resulting  from  the  presence  of  other  electrons. 
Furthermore  it  usually  assumes  that  all  quantities  vary  in  the  Z  direction 
only;  i.  e. ,  in  the  direction  of  beam  flow.  We  are  adopting  these  assiimp- 
tions.  This  is  tantamount  to  assuming  that  the  time -varying  part  of  the 
space-charge  density  function  is  negligibly  small  and  that  we  are  con¬ 
sidering  only  confined  flow.  That  is,  the  electrons  are  constrained  to 
move  in  the  Z  direction  only,  as  if  the  beam  were  of  infinite  cross  section 
and  immersed  in  an  infinitely  strong  magnetic  focusing  field.  These 
assumptions  will  be  discussed  in  more  detail  in  Section  IV,  when  we 
apply  the  theory  to  a  specific  problem. 

Furthermore,  we  assume  the  region  of  electron  flow  to  be  gridded, 
and  consider  the  electric  field  between  these  grids  to  be  constant  in  Z 
and  to  be  a  sinusoidal  function  of  time.  The  electric  field  is  assumed  to 
be  due  only  to  the  voltage  across  the  grids.  The  amplitude  and  phase  of 
the  gap  voltage  is  always  assumed  to  be  known,  whether  supplied  externally 
or  induced  by  a  bunched  beam.  This  will  be  discussed  further  in  Section  IV. 

As  has  been  pointed  out  by  Kleen^^  and  others, the  successive 
integration  of  Newton's  law  to  get  velocity  and  displacement  involves  no 
special  difficulty.  The  difficulty  arises  in  carrying  through  the  method 
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of  quadrature*,  which  requires  the  inversion  of  the  distance  as  a  function 
of  time  to  get  the  time  as  a  function  of  distance.  This  is  then  substituted 
into  the  velocity  function  of  time  to  yield  finally  the  velocity  as  a  function 
of  distance,  i.  e, ,  the  path  of  the  particle.  This  reduction  of  the  problem 
to  quadratures  is  of  frequent  use  in  classical  mechanics,  but  can  be 
applied  successfully  elsewhere  whenever  the  inversion  can  be  carried 
out. 

In  this  case,  the  inversion  process  requires  the  solution  of  a 
transcendental  equation,  and  all  previous  theories  using  this  approach 
approximated  the  inversion  of  the  transcendental  equation  in  various 
ways.  Because  of  its  very  special  form,  however,  this  equation  can 
be  readily  inverted  so  that  a  complete,  exact  solution  can  be  obtained. 

C.  MATHEMATICAL  SOLUTION  OF  THE  TRAJECTORY 

Following  our  previous  treatment,  we  consider  a  single  electron 
incident  from  the  left  into  the  gridded  gap  as  shown  in  Figure  1 . 

V - - - ^ 


Figure  1.  Diagram  of  Ideal  Gridded  Gap. 
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For  a  field  distribution  in  the  Z  direction  of 


E  s  ^  sin (wt  -f  6)  .  (1) 

d 

we  have,  for  Newton's  second  law, 

mZ  -  ^  8in(cdt  +  8)  ,  (2) 

or 

'i  =  ~  8in(s»t  +  0)  (3) 

This  second -order  differential  equation,  together  with  its  boundary 
conditions,  completely  describes  the  path  of  the  particle  in  the  gap.  The 
voltage  is  allowed  to  have  any  magnitude  V  and  any  arbitrary  phase  6 

with  respect  to  the  time  t  =  0  . 

\ 

Integrating  odk|ce  gives  us  the  velocity  equation: 
v(t) 

dZ  s 

v(t^) 


8in(wt  +  0)  dt 


('4) 


V(t)  *  v(t^)  -  [cos  (wt  +  0)  -  cos  (wt^  +  0)j  .  (5) 

In  carrying  out  this  integration,  we  have  chosen  as  lower  limits  the 
initial  conditions  on  the  electron  as  it  entered  the  gap,  and  as  upper 
limits  any  arbitrary  time,  with  its  corresponding  velocity,  that  the 
electron  is  in  the  gap.  For  the  input  gap  of  a  klystron,  the  initial 
velocity  v(tj^)  would  be  merely  the  d-c  velocity  ,  For  the  moment, 
however,  we  wish  to  keep  the  theory  as  general  as  possible,  and  later 
specialize  it  in  its  applications. 
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Integrating  Equation  (5)  again  gives  us  the  distance  equation: 


[ 


cos  (fait  +  8) 


(6) 


Z(t)  =  [v(ta)  +  cos  (wta+  8)]  (t  -  t^) 

eV  r  1 

- - r“  sin  (e»t  +  8)  -  sin  (wt  +8)  (7) 

md«^  '■  »  J 

In  this  integration  we  have  chosen  the  entrance  time  t^  and  the  entrance 
position  Z(t^)  =  0  as  lower  limits  with  arbitrary  upper  limits.  In  this 
notation,  Z(tj^)  means  the  position  at  a  particular  entrance  time  t^ 
for  tl^e  particular  electron. 

Knowing  the  position  and  velocity  as  a  function  of  time  for  this 
electron,  we  can  now  invert  the  distance  relation,  substitute  into  the 
velocity  relation,  and  find  the  velocity  as  a  function  of  distance.  Re¬ 
arranging  Equation  (7),  we  have 

Z(t)  +  [v(t^)  +  cos  (wt^  +  8)]  sinMj^  +  8) 

mdu 


eV 

mdf;^ 


cos  («t^  +  8)j 


t  - 


sin  (wt  ■¥  8) 


(8) 


Ass  tuning  that 

^  <"*a  +  ®)  >  0  •  (9) 

we  divide  both  sides  of  Equation  (8)  by  this  quantity.  We  investigate  this 
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assumption  later.  We  then  have 


eVsin(«tj^+  0) 

md«2[v(g  +  ^  cos  («t^  +  0)] 


(10) 


Before  this  equation  is  amenable  to  inversion,  we  must  perform  some 
algebraic  manipulation.  Multiply  both  sides  of  the  equation  by  u  ,  and 

i 

add  0  to  both  sides.  Since  we  now  are  used  to  the  convention  of  writing 
the  distance  Z  at  some  arbitrary  time  t  as  Z(t)  ,  we  drop  the  t  and 
merely  call  it  Z  ,  although  retaining  exactly  the  same  meaning.  Equation 
(10)  then  becomes 

eVsin(s#t-  +  0) 

(•t  +  0)  +  r - ”  - r - eV - T 

[''<*»>  *  SS  “*  *  *"1  L''"a>  +  'o*  '"•»  *  »>] 


=  («t  +  0)  - 


eV 

■eV 


md«rv(t^)  +  —2 —  cos  (sst  +  0)1 
md(t>  ■* 


sin  (u>t  +  0) 


(11) 


Define  the  following  symbols: 


u  =  (wt  +  0) 


<  = 


mdw 


eV 

r7 


'  eV  T 

v(t^)  +  cos  (srtj^  +  0)J 


(12) 

(13) 
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(14) 


4  =  (mt^  +  6)  + 


wZ 


H'a) 


eV 

mdfa) 


CO 8  (88t^  +  0 


)] 


eV  sin  (ut^  +  0) 


md»rv(t  )  +  cos 

I-  a  md(i) 


(wt^  +  0)j 


Then,  Equation  (11)  may  be  written  in  the  compact  form 


^  =  u  -  c  sinu 


(15) 


In  this  equation,  the  only  term  containing  t  is  u  ,  and  the  only 
term  containing  Z  is  ^  .  Thus,  if  we  can  solve  this  transcendental 
equation  for  u  in  terms  of  ^  and  c  ,  we  will  have,  by  a  short  algebraic 
step,  the  quantity  t  as  a  function  of  Z  and  the  various  Initial  conditions. 

Several  methods  of  approximation^ have  been  used  to  invert 
Equation  (15).  Here  we  shall  develop  an  exact  inversion.  Equation  (15) 
is  Kepler's  equation  from  celestial  mechanics.  It  occurs  often  in  the 
classical  astronomical  literature.  Its  inversion  is  usually  attributed  to 
Bessel,  but  was  originally  developed  by  Lagrange. The  classical 
treatise  on  the  subject  of  celestial  mechanics  is  tl^at  of  Tisserand, 
who  treats  this  equation  in  considerable  detail.  A  more  popular  alter¬ 
nate  derivation  is  given  by  Whittaker, but  we  derive  the  identical  result 
in  a  manner  slightly  different  from  the  classical  authors. 

From  inspection,  one  can  see  that  the  quantity  u  -  4  S'  periodic 

function  of  4  with  period  2ir  .  Expanding  in  a  Fourier  series  give* 

oo 

u  -  t  ^  b^  sin  /  ^  ,  (16) 

/  =  1 


where 
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Integrating  by  parts  gives 


•If 


tr 

-ij  J  cosi4d(u-C) 


-•rr 

but  the  first  term  drops  out  when  evaluated,  so  we  have 

ir  IT 


'1 


■h  J 

-IT  -IT 


COS  /  ;d4 


(18) 


(19) 


The  lart  term  in  Equation  (19)  drops  out  also,  so  that  in  the  final  form 

TT 


1 


hi  -  I  cos  /  4  du  , 


(20) 


or 


IT 

-IT 


cosQ  (u  -  c  sin  u)l  du 


(21) 


Now,  recall  that  the  Bessel  function  of  the  first  kind  is  defined  by 

IT 

J^(X)  =  2^  {  cos  (n  0  -  X  8  in  0)  d  0  , 


(22) 


so  that  in  this  case 


b,  =f  j,(/.) 


(23) 


and  the  inversion  series  is 


00 


u  =  4  + 


^  {U)8inli 
I  =  1 


(24) 


From  the  definition  of  u  in  Equation  (12),  we  have  the  final  result  for 
the  inversion: 


00 


t  =  i  |-0  +  4+  ^  f- (/€)sin/4j  . 

I  9  1 

Note  that  Equation  (24)  can  be  substituted  directly  into  the  velocity 
Equation  (5)  to  give  the  velocity  as  a  function  of  distance: 

00 


(25) 


'^<Z)  =  v(t^).j^ 


cos 


^  J-  (1 « )  »in  1  tj  -  cos  (wtj^  +  e)j 


I  =  1 


(26) 


In  order  to  investigate  assumption  (9)  ,  we  shall  define  4  ^nd 
<  in  a  more  conventional  way.  Let 


signal  voltage  magnitude  _  V 
beam  voltage  ma^itude  “ 


6  = 
^o 


ud 

To 


=  d-c  transit  angle 


Z 

o 


(27) 

(28) 

(29) 
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(30) 


1  +  *(t,) 


Thus  we  have  the  relation  that 


m^u 


(31) 


This  relation  enables  us  to  introduce  convenient,  normalized 
quantities  into  our  relationships,  in  terms  of  easily  measurable  para> 
meters.  Substituting  Equations  (30)  and  (31)  into  definitions  (13)  and 
(14),  we  have 


«  s 


t  s  +  0)  +  f ^ - - - 1 

^o  Za”  <"‘a 


(32) 


a  sin  (wt^  <(■  9) 

*♦0  [*  *  ^  '=*• *  *’] 


(33) 


In  terms  of  these  normalized  parameters,  the  assumption  (9)  is 

K  *  «^‘a^  *  z|“  <"*a  +  ®)]  >  0  (34) 

Since  >  0  ,  we  can  divide  by  without  changing  the  inequality, 
and 


1  * ,(..) .  ^ 


cos  (wt_  +  0)  >  0 

ck 


(35) 
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This  is  the  first  restriction  on  the  validity  of  the  analysis.  Actually, 
we  will  see  shortly  that  a  more  stringent  requirement  than  this  exists 
for  convergence  of  the  series  in  Equations  (24) -(26),  but  at  least  this 
is  a  definite  upper  bound  on  the  analysis.  If  inequality  (35)  is  violated, 
then  it  would  be  possible  to  have  the  time  t  at  the  point  Z  be  infinite, 
meaning  that  the  electron  velocity  would  become  negative  and  the  electron 
would  never  reach  that  point.  When  such  a  condition  is  reached  in  a 
practical  device,  inefficiency  of  operation  would  radically  increase 
because  of  electron  interaction  within  the  gap  itself.  Signal  levels  large 
enough  to  violate  Equation  (35),  then,  are  of  no  interest  at  this  point. 

The  special  case  of  the  input  gap  simplifies  Equation  (35),  for 

then 

''(‘a)  =  K  *  «<‘a)  =  0  •  (36) 

and  Equation  (35)  reduces  to 

•  <  1  (37) 

¥fhen  actually  evaluating  the  series  expressions  in  Equations 
(24) -(26),  it  is  found  that  the  terms  have  different  signs  -  some  positive, 
others  negative.  According  to  Tisserand^^  (p.  262),  the  series  remains 
convergent  for  €  between  zero  and  one,  whatever  the  values  of  may  be. 
That  is,  the  criterion  for  convergence  is 

0  <  €  <  1  .  (38) 

In  the  special  case  of  the  input  gap,  where  Equation  (36)  holds, 
inequality  (38)  reduces  to  the  condition, 

«<6o  (39) 
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In  general,  however,  if  relation  (38)  ie  eatiefied,  the  convergence  holda,  and 


0  <  - = - ^ - -  <  1  ,  (40) 

In  the  aatronomical  literature,  the  quantity  c  is  called  the  ec¬ 
centricity,  and  is  called  the  mean  anomaly.  Much  work  has  been 
done  in  this  field,  some  of  which  is  directly  applicable  to  our  problem. 
The  inversion  problem,  as  done  here,  gives  the  series  with  the  greatest 
applicability.  For  small  values  of  c  ,  however,  the  development  of  u 

as  a  power  series  in  t  is  possible.  This  series  has  been  given  first 
14  17 

by  Lagrange,  also  by  Laplace,  and  follows: 

2 

u  s  1^  4-  c  sin(  -I-  ~  ain2  ( 

3  , 

4  -  1.  ■  -  (3^sin3(  -  3sin() 

3‘2*2^ 

4  —  - - -  (4  sin4|^  •  4*2'^sin2() 

4«3*2*1‘2^ 


4 - * - (5^  sin  -  5  •  3^  sin  3^  4  ain^) 

5*4'3*2‘1«2^  1-2 


4 - ^ - r  (6^  sinbt  -  6’ 4^  a  in  4^  4  2^  sin  2 

6*5»4*3*2»1»2’  1-2 


4  .  .  . 


f  m-1  ,  m  ,  ,vm-l  .  , 

4 -  Im  sinm(^  - -p  (m  -  2)  sin  (m- 2)  4 
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+  (m  -  4)“  ^  sin  (m  -  4)  ;  -  ...  j 


+...  .  (4i) 

The  derivation  of  this  series,  as  given  by  Tisserand^^  (pp.  Z6Z-Z69),  is 
found  in  the  Append  x. 

As  wo(;ld  be  expected,  this  series  is  useful  only  for  relatively 

1  d 

small  values  of  <  .  Indeed,  Laplace  was  the  first  to  shuw  that  this 
aeries  converges  only  if 

I  <  0.6627 _  (42) 

A  simplified  proof  of  this  limit  of  convergence  is  also  given  in  the 
Appendix.  Si.xce  the  application  of  this  theory  will  undoubtedly  be  to 
large-signal  cases,  the  expression  (41)  will  be  of  limited  value  in  com¬ 
putations.  However,  it  may  provide  an  easier  adaptation  to  the  boundary 
conditions  on  a  drift  region  analysis,  and  for  that  reason  it  is  included. 

D.  REDUCTION  TO  SMALL  SIGNALS 

The  expression  (25)  can  be  easily  reduced,  m  the  limit  of  small 
a  ,  to  the  small -signal  expression.  Writing  out  Equation  (25)  in  full, 
gives 

^o  +  «<*a>  +  if"  <^a  + 

a  s'n  («t  +  0) 

_ Ca 
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I  =  1 


2*0  +  «(*a>  ♦  2^ 


a  a  in  (lat^  -f  8) 

^♦o[^  +*<V  +  l5-  co.(#t^ 


(43) 


To  simplify  this  expression,  we  choose  the  following  assumptions: 


e  «  0 


(44) 


g(t  J  *  0 


(45) 


(46) 


(47) 


Then,  setting  Z  «  d  and  t  s  t^  ,  for  the  arrival  time  at  the  exit  plane 
of  the  input  gap,  we  have,  neglecting  order  , 


«tb  -  •«»  =  *0  -  4”  ■*■4" 


-|cos«t^+ O(a^) 

Therefore,  the  exit  velocity'  given  from  Equation  (5)  is 
'^out<'al  "  ^o[*  -  4"  •"*>  *  *o’  ■  '"'a>] 


(48) 


(49) 
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Q 

This  is  essentially  the  small-signal  result  obtained  by  Beck' (p.  353). 

The  assumptions  made  in  Equations  (44) >(47)  are  very  radical;  some 
of  them  can  be  justified  under  all  conditions,  and  some  can  be  justified 
only  under  small-signal  conditions.  A  complete  discussion  of  all  these 
points  is  given  in  Section  III. 

E.  EXIT  CURRENT  AND  KINETIC  POWER 

When  using  ballistic  theory  and  considering  only  current  carried 
by  the  motion  of  electrons,  it  is  common  to  use  the  law  of  conservation 
of  charge  to  compute  current.  In  order  to  justify  this,  we  recall  that 
the  current  crossing  a  surface  S  in  space  is  defined  as  the  rate  that 
charge  flows  across  S  .  We  will  use  this  definition  in  the  following 
derivation. 

Consider  two  planes,  1  and  2,  as  shown  in  Figure  2.  Now  con¬ 
sider  some  small  finite  area  on  planes  a  and  Z  ,  called  AA.  For 
electrons  incident  from  the  left,  the  charge  that  passes  plane  a  in  some 


I 

I 

1 


I 

0 


Figure  2.  Illustration  of  Charge  Conservation. 
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small  time  interval  dt  is  merely  the  number  of  particles  which  pass 

A 

plane  a,  multiplied  by  the  charge  of  each  particle.  But  the  number  which 
pass  plane  a  in  the  infinitesimal  time  dt  is  equal  to  the  number  contained 
in  the  small  volume 

dV  =  v.dt  AA  .  '50) 

a  a  a 

This  is  given  by  n^dV^  ,  and  the  charge  is 

Aq^  =  n^ev^AA  ,  (51) 

where  n  is  the  number  of  particles  per  unit  volume  at  time  t  ,  and 
V  is  the  average  velocity  of  the  particles  at  time  t  .  The  interval 

21  cl 

dt  is  assumed  to  be  so  short  that  n  and  v  are  constant  during  the 

A  A  A 

interval. 

Following  the  same  particles  to  the  plane  Z  ,  we  can  say  in  a 
similar  way  that  the  charge  crossing  plane  Z  in  a  short  time  interval 
dt  is  the  electronic  charge  multiplied  by  the  number  of  electrons  in  the 
small  volume  dV  ,  that  is 

Aq=nevdtAA  .  (52) 

Here,  n  is  the  number  of  particles  per  unit  volume  at  plane  Z  ,  and  v 
is  their  velocity.  These  are  assumed  constant  during  the  interval  dt  , 
and  are  in  general  different  from  n  and  v  . 

From  conservation  of  charge,  these  two  charge  increments  are 
the  same;  therefore 

n  ev  dt  AA  =  nevdtAA  .  (53) 

AAA 

But  current  density  J  is  defined  as  nev,  so  Equation  (53)  may  be 

z 
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written  as 


(54) 


J-  AA  dt  =  J,  AAdt 
a  a  z 

Integrating  over  a  beam  cross  section,  we  have  the  total  current.  If 
the  current  density  is  assumed  not  to  vary  with  transverse  dimensions, 
then 

la^^^a  =  «  (55) 


that  is. 


wherever  the  derivatives  in  Equation  (56)  are  defined.  A  full  discussion 

of  expression  (56)  is  given  in  Section  III.  These  expressions  for  current 

implicitly  assume  an  electron  stream  composed  of  many  electrons  with 

different  entrance  times  t^  ,  so  that  actually  t^  becomes  a  time  vari> 

able  instead  of  a  constant  for  a  single  electron.  The  Equation  (56)  can 

be  written  out  in  full  for  the  current  at  any  plane  Z  as  a  function  of 

entrance  time  t_  : 

a 


I*  la 


1  +  g(ta)  +  2^  cos  (wt  +  0) 


1  +  2g(t^)  +  g^(t J  +  cos  (erta  +  0)  +  2^  g(tj^)  cos  ( ut^  +  0) 

™o  ™o 


sin(ut  ■¥ 
'  a 


e) 


sin^  («t^  + 


e) 


sin  (wtg^  +  0) 
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a>|3 


(v  +  1 )  a 


+  g<t^)  +  C08  («t^  +  0)J 


'v  +  Zi 


(y  ■<•  l  )a  _ 

2<^o  *  4”  *  ®^] 


sin(v  +  l)ci<t^  +  -j — 


1 


[l  +  *(‘a» 
a  sin  (ut3^  ■¥  6) 


y  I  "  COS 


(wt^  +  e)] 


[l  +  g(t^)  +  2^  COS  (<i»t^  +  0)j 


{■ 


♦  2^1+  2g(t^)  +  g^(t^)  +  CO*  +  B) 


+  TT-  g(t  )  cos  (wt  +  0)  +  (55.)  sin(irt-  +  0) 

^  ^  \ZoJ  ^  ‘ 


/«Z  \r£|<!a) 

KA  "•.j 


"  'W 


sin^  (srt-  +  0) 


a  rd*(ta) 


i  sin(wt^  +  0)  >  . 


Jy  +  !< 


(V  +  l)a 


24»o  1  +  g(t^)  +  cos  (wtj^  +  0)j 


cos  {v 


i  1){  wt  +(^1  - - 1 - 

COB  (wt^  + 


»] 


+  0  - 


asin((i)t-  +  0) 


2^^[l  +  g(V  +  2^  COS  +  0)] 


-1 


(57) 
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This  Equation  (57)  is  the  most  general  ballistic  expression  for  the  current 
at  any  point  Z  in  a  gridded  gap.  The  derivation  arises  merely  from 
taking  the  derivative  of  the  expression  (25)  with  respect  to  tj^  . 


At  plane  a  ,  we  have  the  number  of  particles  n^^v^^dt^AA  ,  each 
mv^ 

with  kinetic  energy  — y —  .  We  have  the  same  number  of  particles 


at  plane  Z  ,  but  here  each  has  a  kinetic  energy  of 
gain  in  kinetic  energy  for  these  particles  is 


Thus  the 


AE  =  (nj^Vj^dtj^AA)  ^  (v^  -  vj) 


(58) 


Inserting  the  definition  of  current  density  and  integrating  over  a  beam 
cross  section,  we  have  the  total  kinetic  energy  change  of  the  beam  in 
a  small  interval  of  time,  dt  :  - 

A 


beam 


(v^ 


(59) 


Since  power  is  rate  of  change  of  energy,  we  have  for  the  kinetic  power 


Pa  = 


dE  ,  « 

dT  '  ‘a  n-  I*  - 


m 


(60) 


This  expression  will  be  useful  in  Section  IV,  when  we  wish  compute  the 
average  kinetic  power  gained  or  lost  by  the  beam  over  a  period.  This 
is  found  by  integrating  this  instantaneous  power  function  over  a  whole 
period  and  dividing  by  2ir  . 

In  the  two  sections  following, we  show  the  meaning  of  this  func¬ 
tional  form,  discuss  the  effect  of  varying  some  of  its  parameters,  and 
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specialize  Equations  (26)  and  (57)  to  the  input  and  output  gaps  of  a  two> 
cavity  klystron. 


III.  DISCUSSION 

A.  VARIATION  OF  ENTRANCE  TIME 

In  the  expression  that  has  been  derived  for  the  current,  we  have 
assumed  implicitly  that  more  than  one  electron  is  crossing  the  gap.  It 
is  entirely  reasonable  to  assume  this,  since,  in  the  klystron,  electrons 
are  continually  being  emitted  from  the  cathode  and  are  focused  into  a 
uniform  beam. 

Let  us  first  assume  that  the  electrons  enter  the  gap  with  a 
constant  velocity,  the  d-c  beam  velocity  ,  thus  assuming  that 
g(t^)  -  0  •  This  assumption  is  valid  for  the  input  gap.  The  velocity 
and  current  expressions  are  simplified  by  this  assumption  and  become 


(••t,  +  0  + _ 

b>  Z 

asin(ut  +0) 

1 

Z+zi^cos  («t^  +  0)j 

^♦o[' *  4;“” 1 1 

-  cos  (OJt^  +  0) 


(61) 
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and 


1  +2f-  ':os(wt^+ 


.f 


a  / wZ \  a 

1  +  2^  cos  (wt^  +  0)  +  p in  (<jt^  +  6) 


00 


E[(^j  sin(«tg^+  0)  j 


■9  =  0 


/ 

\ 


(v  +J )  g _ 


T  )  (v  •*•  l)a 


•  sin(i»+  1)  <  wt^+  — 


u)Z 


a  sin  ((i)t  -f  0) 
+  0 - ^ - 


+  2 


K  +  2?;  ("^a  +  ®)]  ^  <«‘a  +  ®)] 

1+j^  «o*  (••tj^  +  ®)+^!|Zj  ^  sin(«t^+0) -|^j|^^sin^(«t^+0)j 


(v-t-  l)a 


*  ^  if"  ^*'*a  *  ®)] 


cos  (»+  l)<wt_  +  0  + 

d 


b>Z 


a  s  in  (wt 


-1 


(62) 


It  may  be  assumed  without  loss  of  generality,  that  for  the  input 
gap,  the  phase  0  is  zero.  This  assumption  is  based  upon  the  fact  that 
there  are  a  uniform  number  of  electrons  passing  plane  a  per  unit  time 
t  That  is,  it  makes  no  difference  whether  the  input  gap  voltage  is  zero 
or  any  other  value  at  t  =  0  ,  insofar  as  the  periodic  behavior  of  the  output 
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velocity  and  current  is  concerned.  This  further  simplifies  expressions 
(61).  (62). 

With  more  than  one  electron  to  consider,  we  naturally  have  more 
than  one  entrance  time  t^^  .  First  of  all,  for  a  particular  choice  of  t^  , 
we  may  see  a  large  number  N  of  electrons  at  the  plane  a  .  Now,  at 
some  time  an  instant  later,  we  shall  have  another  group  of  N  electrons 
at  plane  a  ,  with  a  different  number  for  their  initial  time  t^  .  As  we 
continue  to  increase  real  time,  the  entrance  time  t^^  likewise  increases, 
as  long  as  we  consider  the  different  electrons,  each  group  of  which  has 
its  own  entrance  time  t^  .  In  this  way,  we  can  use  the  time  of  entrance 
of  the  electron  t^^  as  our  real  time  variable.  There  are  an  equal  number 
of  electrons  passing  plane  a  per  unit  time,  for  the  constant  current  and 
constant  velocity  input.  Thus,  we  have  established  the  meaning  of  using 
the  entrance  time  t^  as  a  variable,  and  expressing  functions  and  deri¬ 
vatives  of  functions  In  terms  of  this  variable. 

B.  EXPRESSING  FUNCTIONS  IN  TERMS  OF  EXIT  TIME 

In  a  manner  similar  to  that  indicated  in  Section  A,  the  velocity 
and  current  functions  may  also  be  expressed  in  terms  of  exit  time  tj^  . 

We  have  derived  a  functional  relationship  between  the  exit  time  and  the 
entrance  time  of  an  electron.  If  one  wishes  to  express  the  velocity  and 
current  as  functions  of  exit  time,  it  is  only  necessary  to  solve  Equation 
(11)  for  t^  as  a  function  of  t  .  This  is  extremely  difficult,  and  perhaps 
impossible  for  a  transcendental  equation  as  complex  as  this.  Fortunately, 
it  is  not  iiecessary  to  express  the  current  and  velocity  in  terms  of  the 
exit  time  tj^  ,  as  the  entrance  time  t^^  is  entirely  adequate  for  our 
purposes.  It  should  be  noted,  however,  that  if  the  exit  time  is  desired. 
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it  can  be  easily  calculated  numerically  for  several  values  of  the  entrance 
time  from  expression  (25)  with  Z  =  d  .  The  expression  of  current  and 
velocity  as  functions  of  exit  time  is  thus  quite  possible  in  graphic  or 
numerical  form.  For  a  choice  of  t^^  ,  the  corresponding  tj^  can  be 
found,  and  then  the  current  and  velocity  can  also  be  found  for  that 
particular  t  .  Then  a  plot  can  be  readily  made  of  current  and  velocity 
versus  t^,  . 

There  may  be  an  apparent  anomalous  relationship  between  t 

Ck 

and  t|^  that  makes  it  difficult  to  express  current  and  velocity  in  terms 
of  t|j  .  This  difficulty  is  that  for  very  large  signals  it  would  seem  that 
crossover  could  occur  within  the  gap  to  some  extent,  so  that  for  a  given 
exit  time  there  could  be  several  corresponding  entrance  times.  That 
is,  several  groups  of  electrons  that  entered  the  gap  at  different  times 
could  conceivably  leave  the  gap  at  the  same  time.  This  would  involve, 
then,  the  entrance  time  being  a  multivalued  function  of  the  exit  time,  so 
that  for  a  given  exit  time,  there  would  be  several  values  for  the  entrance 
time  and  no  unique  entrance  time. 

When  this  difficulty  occurs,  as  it  undoubtedly  does  at  some  ntodu- 
lation  level,  the  rigorous  current  relationships  derived  in  Equations  (50)  • 
(57)  are  no  longer  true.  One  cannot  simply  take  a  derivative  of  a  multi¬ 
valued  function,  or  invert  these  derivatives.  Instead,  the  computation 
of  the  current  corresponding  to  the  exit  time  t^^  must  be  carried  out 
through  a  summation  process,  adding  up  all  the  incremental  charges 
that  had  different  entrance  times  t^  but  that  have  the  same  exit  time  t|^  . 

Fortunately,  we  do  not  need  to  concern  ourselves  with  such  large 
signal  levels  here.  We  can  have  a  large  degree  of  velocity  change  and  of 
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current  bunching  within  the  gap  without  the  crossover  actually  taking 
place.  Indeed,  one  might  expect  the  point  of  the  divergence  of  the  series 
in  Equation  (25)  to  be  the  point  of  crossover.  To  investigate  this  point, 
we  use  a  graphic  technique. 

Consider  the  Kepler  equation 

4  =  u  -  c  sinu  (63) 

For  the  case  where  c  =  1  ,  the  series  diverges.  In  this  case,  we  have 

I 

the  behavior  illustrated  in  Figure  3.  Now  we  can  see  that  the  slope  zero 
is  indeed  the  reason  for  the  divergence,  for  when  the  equation  is  inverted, 
as  in  Figure  4,  an  infinite  slope  r  esults.  If  the  case  proceeded  for  larger 
values  of  c  ,  then  upon  inversion  there  would  be  no  unique  value  for  u 
for  a  given  (  .  Indeed,  there  could  be  three  values  of  u  for  the  properly 
chosen  value  of  C  • 

To  illustrate  that  this  point  of  series  divergence  is  not  necessarily 
the  point  of  maximum  bunching,  let  us  make  a  rough  graph  of  the  deri¬ 
vative  of  u  with  respect  to  .  This  is  given  in  Figure  5. 

This  derivative  goes  to  infinity  at  ^  =  0  ,  2ir  ,  4ir  ,  etc.  Note 
that  for  values  of  f  <  1  ,  the  derivative  does  not  go  to  infinity  at  these 
points  but  remains  finite  and  positive  for  all  ^  .  For  values  of  c  >  1  , 
the  slope  of  the  curve  u  versus  ^  becomes  negative  for  (  =  0  , 
indicating  that  the  derivative  function  in  Figure  5  would  become  negative. 
The  behavior  of  this  derivative  is  important,  because  it  is  related  to 
the  current.  Letting  u  =wt|^  and  6  =  0  for  an  input  gap,  we  have 


^  d; 

dtj^  ~  d(«tj^)  "  d;  d(«tj^) 


(64) 
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Figure  5.  ^  versus  ^  for 


The  current  as  derived  in  Equation  (56)  is  merely  the  inverse  of 
Equation  (64).  In  order  to  know  where  zeros  and  infinite  peaks  occur  in 
the  current,  one  nM;i|,<only  know  the  behavior  of  the  two  derivatives 
du/d4  and  d4/d(«*t  )  . 

3L 

Examining  the  expression  for  c  in  Equation  (32),  we  see  that 
c  approaches  one  most  rapidly  with  increasing  a  if  the  time  value  wt^ 
is  equal  to  ir  .  For  c  =  1  and  wt^  =  ir  ,  the  restriction  in  Equation  (39) 
can  be  derived.  We  would  therefore  like  to  examine  the  current  at  this 
first  point  of  series  divergence.  Knowing  the  approximate  behavior  of 
du/d(^  ,  we  calculate  d4/d(Mt^)  by  differentiating  Equation  (33): 


dt  -  _ 1 

'“'"‘a' 


coswt^  +  —  sincdt 
a  2  a 


(65) 


From  Equation  (32),  we  can  deduce  that  if  c  s  1  and  wt^  is  set 
equal  to  tr  ,  then  this  determines  that  a/2^^  s  l/2  .  Using  this  relation 
in  Equation  (65)  together  with  wt^^  =  ir  we  find 


(66) 


Therefore,  d(/d(«*t^)  is  finite  and  not  zero  for  these  values. 

We  know  that  du/d(  is  positive  for  all  values  of  ^  when  c  =  1  ,  and 
therefore  we  find  the  product  (du/d^)  Cd&/d(wt  )^  is  nonzero  for  these 

cl 

values.  Thus  there  is  no  crossover  at  c  =  1  ,  wt.  =  ir  ,  even  though  the 
series  diverges.  We  can  see  from  the  graph  of  du/d(  that  nothing  un¬ 
usual  happens  to  this  derivative  as  c  increases  from  0  to  1  .  It  is 
always  positive  and  finite.  Therefore,  we  expect  any  erratic  behavior 
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of  the  current  to  be  due  to  the  zeros  of  di^/d(«it  ). 

3l 

3 

This  erratic  behavior  is  well  illustrated  if  one  makes  6  10  , 

o 

a  10  .  Then,  in  Equation  (65),  terms  containing  the  ratio  a/2^^  will 
be  negligible,  and  the  equation  can  be  written 

=  1  +  y  sinwt^  (67) 

3i 

Here,  t  <<  1  for  all  time  so  that  du/d^  is  positive.  If  a  >  2  ,  zeros 
of  Equation  (67)  are  periodic.  Therefore,  we  must  have  peaks  of  in¬ 
finite  current  appearing  periodically  at  the  exit  plane  of  this  very  long 
gap.  This  means  that  for  a  very  long  gap,  this  ballistic  theory  will 
predict  velocities  and  currents  beyond  crossover  and  will  converge  fairly 
rapidly. 

We  have  therefore  found  a  region  in  which  the  series  diverges 
but  crossover  has  not  yet  occurred,  and  a  region  in  which  crossover 
has  occurred  but  the  series  still  converges  rapidly.  There  are  many 
intermediate  steps  between  these  two  extremes.  We  have  three  para¬ 
meters  that  determine  crossover  and  series  divergence:  a  ,  a/2^^  , 
and  wtj^  .  If  these  are  fixed,  then  one  can  determine  both  series  conver¬ 
gence  and  crossover  at  that  point.  Actually,  only  «st  and  a/2^  are 
necessary  to  fix  the  series  convergence.  If  the  relation. 


(68) 


is  stipulated,  then  the  series  will  converge,  since  c  will  be  less  than 
unity.  Then  the  a  chosen  will,  together  with  the  other  two  parameters, 
determine  crossover.  The  most  rapid  way  to  check  this  is  by  graphic 
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or  computer  eolutione  of  the  zeros  of  Equation  (65).  If  inequality  (68)  is 
satisfied,  then  the  series  expressions  will  hold  even  if  crossover  has 
occurred.  However,  in  klystron  gaps,  the  gap  length  is  usually  too  short 
to  allow  crossover,  and  if  crossover  did  occur,  the  space-charge  forces 
in  a  bunch  would  caase  a  deviation  from  this  ballistic  theory. 

.Perhaps  we  should  emphasize  that  there  is  a  whole  range  of 
values  of  a/2^^  >  1/2  for  which  the  series  converges  for  some  wt^  < 

The  behavior  of  the  beam  for  such  large  drive  levels  can  be  found  by 

/ 

this  convergent  series  only  for  ir  -  5  >  ti>t^  >  ir  -f  6  ,  where  6  is  a 
small  parameter,  which  is  related  to  how  much  a/2^^  exceeds  1/2. 

If  a/2^^  is  less  than  l/2  ,  however,  the  series  converges  for  all  time. 

We  can  conclude,  therefore,  that  this  series  converges  to  a 
definite  result  for  all  values  of  time  if  a  <  .  If  a  ^  ,  then  the 

series  will  converge  for  all  time  except  when  wt^  is  in  a  neighborhood 
of  ir  ,  the  size  of  the  neighborhood  being  determined  implicitly  by 
Equation  (68).  There  seems  to  be  no  connection  between  crossover  of 
electrons  and  divergence  of  the  series.  The  series  can,  in  many  cases, 
give  a  convergent  result  for  the  region  beyond  crossover.  Whether  or 
not  crossover  also  occurs  in  the  region  of  series  divergence  is  not 
known.  We  can  now  proceed  to  evaluate  these  equations  with  confidence 
within  their  region  of  convergence.  The  ultimate  criterion  for  this 
theory,  of  course,  is  its  actual  numerical  evaluation. 

C.  MEANING  OF  A  VARIATION  IN  ENTRANCE  VELOCITY 

As  we  have  indicated  already,  a  possible  application  of  this 
theory  is  to  the  output  gap  of  a  high-power  klystron  to  determine  the 
velocity  and  current  at  the  exit  plane  of  this  gap.  This  information  can 
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be  used  in  the  present  problem  of  eliminating  some  of  the  harmful  X-rays 
that  are  produced  when  high-velocity  electrons  strike  the  anode. 

To  apply  this  theory  to  an  output  gap,  we  must  keep  the  most 
general  initial  conditions  of  arbitrary  velocity  and  current  at  the  entrance 
to  the  gap.  In  the  ballistic  theory,  the  current  at  the  entrance  plane  enters 
into  the  calculation  of  the  exit  current  only  as  a  multiplying  factor.  The 
velocity  at  the  entrance  plane,  however,  significantly  affects  the  transit 
time,  exit  velocity,  and  exit  current.  The  time -varying  part  of  this 
initial  velocity,  normalised  against  the  d-c  beam  velocity,  appears  as 
the  quantity  g(t^)  in  the  preceding  analysis. 

Perhaps  a  detailed  explanation  of  the  meaning  of  this  term  g(t  ) 
should  be  made.  In  the  preceding  section,  we  considered  the  extension 
of  the  ballistic  treatment  for  an  electron  with  a  single  entrance  time  t^ 
and  a  single  entrance  velocity  We  extended  the  theory  to  different 

electrons,  each  with  its  own  different  value  of  t.  ,  but  each  was  assumed 

A 

to  have  the  same  velocity,  the  d-c  beam  velocity  of  .  There  is  no 
reason  why  each  electron  must  have  the  same  velocity,  since  the  analysis 
was  made  for  an  arbitrary  entrance  velocity.  If  each  electron  had  a 
different  value  of  v(tj^)  for  its  own  entrance  time  t^^  ,  the  analysis  would 
still  be  valid.  We  therefore  may  assume  a  velocity  function  v(t^)  that 
varies  slightly  about  the  d-c  beam  velocity.  The  term  g(t  )  is  defined 

A 

from  this  in  Equation  (30). 

One  may  perhaps  object  that  the  initial  conditions  at  the  entrance 
plane  of  the  gap  are  not  the  same  with  a  velocity  variation  as  well  as  an 
assumed  entrance  time  variation.  This  is  indeed  true,  in  one  respect. 

The  number  of  electrons  per  unit  time  that  enter  the  gap  may  be  different 
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for  each  value  of  .  Furthermore,  the  velocity  of  these  electrons  may 
be  different  for  each  value  of  t^  .  This  means  that  the  arbitrary  phase 
6  of  the  applied  gap  voltage,  which  we  could  choose  to  be  zero  in  the 
input  gap  case,  can  no  longer  be  chosen  as  zero.  Indeed,  it  makes  a 
significant  difference  whether  6  is  chosen  so  as  to  make  the  applied 
voltage  in  phase  or  out  of  phase  with  respect  to  the  entrance  current. 

But  the  phase  6  must  be  given  as  one  of  the  known  measurable  para¬ 
meters. 

Insofar  as  the  electrons  themselves  are  concerned,  it  makes 
absolutely  no  difference,  in  the  ballistic  theory,  whether  a  single 
electron  or  thouiaads  of  electrons  are  entering  the  output  gap  at  a  par¬ 
ticular  time  t^  .  The  ballistic  theory.as  we  have  formulated  it.does 
require,  however,  that  all  electrons  entering  at  one  particular  time, 
tj^  ,  have  the  same  velocity.  This  neglects  any  fluctuations  resulting 
from  noise  in  the  beam  current. 

With  the  extension  of  the  analysis  to  the  case  of  arbitrary  entrance 
velocity  and  current  one  must  still  keep  the  general  requirement  that  c 
be  less  than  unity.  If  this  condition  is  satisfied,  the  series  solutions  will 
be  convergent.  This  convergence  requirement  cannot  be  reduced  easily 
to  the  simple  form  given  for  the  input  gap,  but  can  be  easily  checked  in 
the  numerical  computations  of  a  particular  case.  In  the  next  section, 
the  explicit  method  of  calculation  will  be  shown  for  the  general  case  of 
an  output  gap. 
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IV.  APPLICATION  TO  X-RAY  SUPPRESSION 


A.  INFORMATION  REQUIRED  FOR  X-RAY  SUPPRESSION 

The  simple  small -signal  klystron  theory  predicts  that  the  elec¬ 
trons  that  are  not  in  a  bunch  at  the  output  gap  will  be  out  of  phase  with 
the  voltage  induced  on  the  last  gap  and  thus  will  be  accelerated  instead 
of  retarded  in  this  gap.  Since  the  output  gap  voltage  is  usually  of  the 
order  of  the  beam  voltage,  the  out-of -phase  electrons  will  have  veloc¬ 
ities  roughly  equivalent  to  twice  the  beam  voltage.  For  beam  voltages 
of  the  order  of  magnitude  of  tens  and  hundreds  of  kilovolts,  very  hard 
X-rays  are  produced,  which  are  extremely  dangerous  for  operating 
personnel.  Current  methods  of  solving  this  problem  by  shielding 
techniques  using  lead  and  other  metals  have  resulted  in  very  bulky 
equipment.  The  proposed  method  has  been  to  limit  the  X-ray  produc¬ 
tion  by  slowing  down  the  fast  electrons. 

To  apply  this  method,  information  is  needed  about  the  velocity 
and  current  distribution  of  the  beam  as  it  leaves  the  output  gap.  If 
the  behavior  of  the  beam  is  known  one  may  be  able  to  place  a  structure 
after  the  output  gap  to  extract  energy  from  the  fast  electrons.  The 
velocity  of  the  electrons  is  important,  for  if  the  electrons  are  of  the 
same  order  of  magnitude  in  velocity  throughout  a  cycle,  then  the 
simplest  solution  would  perhaps  be  a  d-c  retarding  field.  On  the  other 
hand,  if  the  velocity  and  current  distributions  vary  widely  with  time, 
then  perhaps  an  r-f  structure  should  be  used,  with  the  magnitude  and 
phase  of  the  applied  voltage  specified.  To  design  such  a  structure, 
the  current  as  well  as  the  velocity  must  be  known,  so  that  the  number 
of  electrons  in  a  particular  cycle  that  possess  the  same  velocity  may  be 
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calculated.  In  this  way,  one  can  obtain  information  about  the  intensity 
of  the  X-riays  produced,  as  well  as  their  wavelength. 

B.  LARGE  SIGNAL  SPACE -CHARGE  WAVE  THEORY  IN  DRIFT  REGION 
Very  few  large -signal  analytical  results  for  the  drift  region  have 
been  published  to  date.  Some  machine  computation  has  been  carried  out, 
but  the  most  prominent  large -signal  analytical  efforts  have  been  those 
of  Paschke.^'^'^  In  his  first  paper,  Paschke^  derives  the  nonlinear  space - 
charge  wave  equation.  He  solves  this  equation  by  the  method  of  successive 
approximations,  extending  his  results  to  the  third  harmonic  in  his  most 
recent  work.^  Rather  than  develop  a  similar  approach  to  the  problem  or 
attempt  a  closed -form  solution,  we  will,  for  purposes  of  computation, 
use  Paschke's  most  recent  results  directly.  His  derivation  of  the  non¬ 
linear  space -charge  wave  equation  assumes  one  very  important  simpli¬ 
fication;  that  the  d-c  beam  velocity  does  not  vary  with  distance.  If  Z 
is  the  distance  variable  in  the  drift  tube,  then  his  assumption  is 


8Z 


S  0 


This  assumption  implies,  from  the  equation  of  continuity,  that 


(63) 


8po 


8Z 


0 


(64) 


There  is  no  reason  to  suppose  that  this  is  the  case  for  extremely  large 
signals,  for  the  pbrtion  of  the  total  velocity  that  does  not  vary  with  time 
could  conceivably  change  in  magnitude  with  the  distance  Z  .  However, 
this  assumption  is  necessary  in  order  to  put  the  nonlinear  equation  into 
a  workable  form.  One  should,  however,  remember  that  the  solutions 
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are  subject  to  this  limitation,  and  therefore  cannot  give  correct  results 
if  carried  to  higher  orders  indefinitely.  Paschke's  third-order  result 
is  probably  the  best  nonlinear  analysis  yet  published  for  the  drift  region. 

For  the  general  case  of  a  finite  diameter  beam,  Paschke^  gives 
the  following  results  for  the  current  and  velocity: 


1.  Fundamental  Current 
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2.  Fundamental  Velocity 
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Second -Harmonic  Velocity 
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5.  Third -Ha rnaonic  Velocity 
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In  these  equations, 


a  =  p(«)  PpX 

(71) 

ft  =  P(«) 

^  pTZS] 

(72) 

o-  r  -EH 
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and  p(nw)  is  the  plasma  frequency  reduction  factor  of  the  fundamental 
mode  at  the  n^  harmonic  frequency.  Paschke's  initial  condition  at 
X  s  0  is 

V  ^^1 

—  sinwt  (74) 

o  o 

To  adapt  Paschke's  results  to  the  special  case  of  an  infinite  beam, 
we  set  I  s  <r  =  p(u)  =  1  .  Then  his  results  simplify  considerably.  All 
of  the  terms  containing  the  factor  (  1  -  )  in  the  numerator  drop  out. 

This  drops  the  last  three  terms  from  the  fundamental  current,  five 
terms  from  the  fundamental  velocity,  a  term  from  the  second  harmonic 
current  and  velocity,  three  terms  from  the  third  harmonic  current,  and 
four  terms  from  the  third  harmonic  velocity.  In  the  final  simplification 
we  have  the  current  as 
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and  the  velocity  as 


^  Pp*  -  -JT  [‘=°*  V  ■  ‘=°®  ^ 

+  ^  (J)^(~)(«in2  PpX)ein2(«t-PgX) 

+  ^(5)^3^^  (cosppX-cos3Ppx]sin3(s»t-P^x)  (76) 

These  results  will  be  used  to  provide  the  initial  conditions  for  the  bal¬ 
listic  treatment  of  the  output  gap  in  the  subsequent  calculations. 

C.  LARGE -SIGNAL  BALLISTIC  THEORY  IN  THE  OUTPUT  GAP  REGION 
To  apply  the  large -signal  ballistic  theory  to  the  output  gap,  it 
is  necessary  to  keep  the  most  general  form  for  the  velocity  and  current, 
as  is  given  in  Equation  (26)  and  (57).  This  results  in  a  somewhat  un¬ 
wieldy  expression  for  numerical  calculation;  therefore  the  results  are 
obtained  most  easily  from  machine  computation.  The  results,  however, 
may  be  carried  to  any  desired  degree  of  accuracy  in  this  way.  Subject 
to  the  restrictions  of  the  preceding  section  for  convergence  of  the  series, 
the  general  ballistic  theory  yields  the  current  and  velocity  at  any  plane 
in  the  output  gap. 

One  might  inquire  into  the  basic  assumptions  of  the  theory  at 
this  point;  that  is,  just  how  valid  is  the  ballistic  theory  for  an  output 
gap?  The  answer,  of  course,  is  that  it  is  an  approximation  to  the 
actual  behavior,  an  approximation  that  neglects  space -charge  variations 
across  the  gap.  This  theory  does  allow  for  a  variation  in  the  quantity 
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of  charge  per  unit  time  that  enters  the  gap.  It  merely  assumes  that  the 
potential  depression  at  any  instant  across  the  gap  is  negligible  compared 
to  the  potential  caused  by  the  alternating  voltage  applied  to  the  grids. 

For  short  gaps,  of  the  order  of  zero  to  two  radians  in  length,  there  is 
no  reason  why  these  assumptions  should  not  be  a  relatively  valid  approx¬ 
imation  of  the  actual  behavior. 

These  assumptions  are  further  validated  by  the  following  con¬ 
sideration.  Recall  that  in  the  mathematical  statement  of  the  problem, 
we  assumed  that  there  was  a  known  voltage  applied  to  the  grids  of  the 
arbitrary  gridded  gap.  All  that  was  actually  assumed  was  that  the 
voltage  existed  there  on  the  grids.  This  voltage,  in  the  case  of  the 
input  gap,  was  supplied  from  the  outside  and  was  assumed  to  be  known, 
but  for  the  ouQ>ut  gap,  the  cavity  is  usually  tuned  to  resonance  at  or 
near  the  fundamental  frequency  and  extracts  power  from  the  beam. 

That  is,  the  voltage  appearing  on  the  grids  of  the  gap  is  not  supplied 
from  an  external  source,  but  is  induced  on  the  grids  by  the  bunched 
beam.  The  beam,  in  so  doing,  tives  up  power  to  the  external  circuit 
and  changes  its  values  of  current  and  velocity  across  the  gap.  Now, 
because  of  the  induction  of  the  voltage  on  the  grids  by  the  beam  and 
because  of  the  voltage  producing  a  change  in  beam  current  and  velocity, 
wexan  consider  this  behavior  as  a  reciprocal  relationship,  and  know¬ 
ledge  of  the  effect  of  the  voltage  on  the  beam  is  equivalent  to  knowledge 
of  the  effect  of  the  changing  beam  upon  the  induced  voltage.  The  reason 
for  this  equivalence  is  that  there  can  be  only  one,  unique  voltage  on  the 
grids  for  a  given  beam  behavior,  and  conversely,  there  can  be  only  one, 
unique  change  in  the  beam  configuration  resulting  from  a  given  voltage 
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function.  This  change  in  the  beam  resulting  from  the  known  voltage 
function  is  determined  by  Newton's  second  law  together  with  the  boundary 
conditions.  Therefore,  if  the  magnitude  and  phase  of  the  output  voltage 
is  measured,  the  ballistic  theory  using  this  information  should  give 
fairly  reliable  results  when  applied  to  the  output  gap. 

D.  CALCULATION  PROCEDURE 

The  calculation  proceeded  in  eight  steps. 

1 .  Input  Current 

^c  =  ^  <1^)  [-inPpX  -y  (J)^  '  Pe*) 

+  j(|)^  -cos2ppx]cos2(trt-P^x) 

+  -^(y)  [3sinppX-sin3Ppx]cos3(«t-p^x)  (77) 


2.  Input  Velocity 


Vc  =  1  +  g  =  1  +  (|)|cos  PpX  -  ^  (|)^  [cos  PpX  -  cos  3Ppx]| 


/  *  sin  (ut  -  p^x) 


+  j  (|  )^  (j^)  [sin2  PpxJ  sin2  (wt  -  p^x) 

+  (2)^(2;^)^  [cosPpX-cos  3ppx]sin3(«t-p^x)  (78) 


These  two  basic  initial  conditions  were  calculated  from  Paschke's  results. 
Equations  (75)  and  (76).  The  velocity  is  used  in  subsequent  computations 
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of  transit  time. 

3.  Transit  Time 
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With  the  calculation  of  the  transit  time,  the  final  velocity  calculation  is 
simplified. 


4.  Final  Velocity 
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(80) 


In  order  to  calculate  the  final  current  at  the  exit  plane  of  the  output  gap, 
one  needs  the  derivative  of  Paschke's  velocity. 
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5.  Input  Velocity  Derivative 
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Now  one  can  BubBtitute  Bteps  2  and  5  into  the  expresBion  for  the  derivative 
of  the  traneit  time. 


6.  Traneit  Timr  Derivative 
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The  calculation  of  the  final  current  is  then  easily  performed. 


7.  Final  Current 
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In  order  to  find  the  average  power  added  to  the  beam  per  period  as  the 
phase  of  the  output  voltage  changes,  we  integrate  Equation  (60). 

8.  Average  Kinetic  Power  Gain  of  the  Beam  per  Period 

2ir 

Pav'Z?  I  “("'c'  ■ 

o 

The  computation  was  done  on  a  Burroughs  220  with  punched -card 
input,  and  the  program  was  designed  so  that  any  of  the  six  parameters 
a,  a^,  6,  could  be  given  any  values  desired.  Here,  the 

subscripts  1  refer  to  the  normalized  voltage  and  drift  length  of  the 
second  cavity.  Each  step  was  printed  out  for  16  values  of  wt^  between 
zero  and  2ir  . 

Several  comments  should  be  made  about  the  calculations.  In 
steps  1  and  2,  the  computation  of  the  input  current  and  input  velocity 
from  Paschke's  theory,  it  is  necessary  to  choose  values  of  and 

a  so  that  the  product  (a/2)  (w/w^)  is  always  less  then  unity.  If  this 
condition  is  violated,  then  the  magnitudes  of  the  second  and  third  har¬ 
monics  could  be  greater  than  the  fundamental  for  some  values  of  B  x  . 

P 

This  implies  that  the  assumed  series  solution  that  Paschke  uses  may 
diverge,  and  that  the  expression  for  current  and  velocity  up  to  the  third 
harmonic  is  an  extremely  poor  approximation  to  the  actual  beam  be¬ 
havior.  For  our  purposes,  we  used  the  values  of  a  =  .  3  and  w/w^  =  5  , 
These  values  gave  a  product  of  .  75  ,  which  was  significantly  less  than 

one.  We  wished  to  choose  the  large  value  of  and  the  value  of  a 

P 

to  correspond  more  closely  to  the  parameters  used  in  modern  high- 
power  two -cavity  klystron  amplifiers,  rather  than  to  the  parameters  in 
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the  domain  of  strict  small-signal  theory.  The  value  of  p  x  was  com- 
puted  from  the  relation, 

P 

Since  most  two-cavity  Jdystrons  are  designed  so  that  the  drift 
length  is  a  quarter  of  a  plasma  wavelength,  the  value  of  p^x  =  ir/Z  was 
chosen.  This  makes  the  current  have  the  maximum  amplitude,  and 
makes  the  velocity  unity  for  all  time;  that  is,  the  velocity  at  p^x  =  ir/Z 
is  merely  the  d-c  velocity.  If  the  value  of  (a/2)(u)/up)  is  larger  than 
unity,  then  the  resulting  current  will  be  negative  for  some  values  of  ut 
and  positive  for  other  Values  of  ut .  Since  the  velocity  is  always  a  constant 
positive  unity,  this  corresponds  to  the  physically  impossible  condition,  that 
the  charge  density  change  sign  as  a  function  of  time.  The  condition, 

.  (86) 

P 

must  always  hold. 

In  the  computation,  the  steps  2  and  5  were  included  for  generality, 
although  they  are  one  and  zero  respectively  for  all  values  of  ut^  for  the 

P  X  we  have  chosen. 

P 

In  the  evaluation  of  the  series  in  steps  3  and  6,  the  criterion  was 
included  that  the  argument  of  the  Bessel  function  must  be  less  than  the 
order.  This  is  equivalent  to  condition  (40),  and  insures  convergence  of 
the  series.  The  further  restriction  was  imposed  that  the  niimber  of  terms 
in  the  series  should  not  exceed  50.  Neither  of  these  restrictions  were 
violated  for  the  values  of  the  parameters  chosen. 
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The  Bessel  functions  were  computed  from  the  power  series 
expression  for  small  arguments  given  in  Equation  (47).  Accuracy  was 
obtained  to  at  least  ten  significant  figures  by  this  method.  Accuracy  of 
the  total  computation  was  compromised  to  five  figures,  since  a  graphic 
presentation  of  the  results  was  to  be  given. 

Paschke's  expressions  for  the  current  and  velocity  use  wt  , 
the  time  used  to  describe  the  voltage  at  the  input  gap,  as  an  independent 
variable.  There  is  a  certain  phase  delay  that  enters  into  the  expression 
for  the  current  and  voltage,  if  one  wishes  to  express  them  as  functions 
of  time  of  arrival  at  the  entrance  plane  of  the  output  gap.  Since  Paschke 
assumes  that  the  d>c  beam  velocity  is  constant  in  space  and  time,  we  can 
deduce  that  the  difference  between  expressing  the  output  gap  entrance 
current  and  entrance  velocity  as  a  function  of  wt  or  as  a  function  of 
(st^  is  merely  some  phase  constant.  In  the  calculations,  therefore,  it 
was  assumed  that  wt  s  ,  and  this  phase  constant  was  implicitly 
contained  in  the  phase  6  .  Strictly  speaking,  therefore, the  meaning  of 
the  constant  6  was  slightly  different  in  the  calculation  scheme  than  in 
its  original  derivation.  However,  0  was  varied  between  zero  and  2tr  , 
and  results  were  therefore  obtained  to  cover  all  possible  cases.  It 
still  remained  to  find  out  what  the  relationship  was  between  a  particular 
6  and  the  tuning  of  the  output  cavity.  Therefore,  the  change  in  kinetic 
power  was  computed  for  each  value  of  wt^  and  was  integrated  over  a 
period,  giving  the  average  power  added  to  the  beam.  The  phase  6 
which  gives  the  most  negative  value  for  this  average  power  was  assumed 
to  be  the  operating  point  of  the  klystron  amplifier  under  normal  condi¬ 
tions;  that  is,  usually  the  output  cavity  is  tuned  to  resonate  at  the 
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fundamd'ritail  frequency  ftnd  esktracit  the  maximum  amount  o£  from 

the  beam.  It  kihould'be  pointed  out  that  the  express  ion. 


'  'a  <''b  - 


(87) 


was  used  in  the  calculations,  but  this  was  merely-  a  convenient  Choice. 
We  could  have  used,  with  the  same  validity,  the  expression, 

Pb ' ‘b 


These  functions  of  course  are  different,  but  when  integrated  over  a 
period,  they  are  equivalent.  This  is  easily  shown  as  follows. 

From  Equation  (59), 

«  =  'a'“aS<''b  -'^1  •  <*” 


for  the  kinetic  energy  change  of  the  beam  across  ^o  planes,  a  and 
b  .  But 


V*a  -  V‘b 


(90) 


therefore 


(91) 


where  and  represent  the  total  current  at  planes  a  and  b  , 
respectively.  Then  we  can  define 


dE  T  m  .  2  2. 

Pb  =  3t^  =  ^b  ZF  ^'^b  ■  "^a^ 


as  a  function  of  ett. 


But 


Pa  'a  ‘^'b 


(92) 


(93) 
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(94) 


and,  changing  notation  to  the  output  planes  c  and  d  ,  we  have 

2ir  2ir 

2¥  J  =  Z¥  J"  ^d‘*<"‘d) 


(95) 


This  equation  is  the  result  desired.  Strictly  speaking,  the  average 

power  is  equal  to  the  value  computed  in  step  8,  multiplied  by  . 

The  currents  and  velocities  in  the  steps  of  calculation  are  all  total 

currents  and  total  velocities,  normalized  against  the  d-c  beam  current 

I  and  the  d-c  velocity  t  Therefore,  the  normalization  is  also 
o  o 

implicitly  included  in  the  power  calculation,  as  simple  algebra  will  show. 

In  the  actual  integration  to  find  average  power,  only  approximate 
values  were  necessary  for  the  integral,  so  the  trapezoidal  rule  was 
used  with  16  equally  spaced  intervals  per  period. 

The  following  values  of  the  six  parameters  and  the  time  variable 
were  used: 

o  =  .  3  , 


ppX  =  ir/2 
^2  =  1  •  5  , 

Cj  =  0.6,  0.7,  0.8,  0.9,  1.0,  1.  1 
0  =  0,  ir/8,  2ir/8,  3^/8,  ...,  15ir/8  , 

wt^  =  0,  ir/8,  2ir/8,  3ff/8,  .  .  ,  ,  15it/8 
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E.  RESULTS 


For  each  possible  combination  of  the  parameters  used,  all  eight 
steps  were  evaluated  and  printed  out.  The  most  important  results  are 
the  transit  time,  the  final  velocity  and  current,  and  the  average  power. 
Each  set  of  parameter  values  yielded  a  unique  value  of  arrival  time  at 
plane  d,  ut^  •  Since  we  are  interested  primarily  in  the  current  and  vel¬ 
ocity  at  plane  d  as  a  function  of  time,  we  must  express  most  of  the 
graphic  results  using  wt^  as  a  time  variable.  We  must  do  this  because 
we  are  eventually  interested  in  the  beam  behavior  beyond  the  output 
gap,  at  the  collector.  A  structure  will  eventually  be  placed  in  the  region 
between  the  output  gap  and  the  collector  which  will  reduce  the  velocities 
of  the  fast  electrons  in  the  beam.  In  order  to  provide  the  entrance  con¬ 
ditions  to  such  a  structure,  we  must  express  our  current  and  velocity  in 
terms  of  the  entrance  time  to  this  structure,  which  is  the  same  as  the 
exit  time  of  the  output  gap.  The  situation  here  is  essentially  analogous 
to  that  encountered  i providing  proper  entrance  conditions  to  the  output 
gap.  In  that  case,  we  expressed  Paschke's  current  and  velocity  in  terms 
of  entrance  time  wt^  and  included  the  constant  phase  difference  between 
(dt^  and  wt  in  the  parameter  9  .  However,  in  this  case,  the  phase 

difference  between  wt ,  and  wt  is  not  constant,  but  is  a  function  of 

o  c 

entrance  time  wt^  .  It  is  contained  in  the  transit  time  expression. 
Therefore  we  must  include  it  explicitly  by  evaluating  it  numerically  and 
then  plotting  current  and  velocity  as  functions  of  exit  time  wt^  . 

The  values  of  exit  current  versus  exit  time  for  four  different 
values  of  phase  are  plotted  in  Figure  6. 

ti  '  '  '  * 

From  this  point,  "time"  refers  to  normalized  time;  that  is  the  angle 

wt^  or  wt,  . 
c  a 
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Figure  7.  Average  Normalized  Kinetic  Power  per  Period  Added 
to  Beam  versus  Phase  of  Output  Voltage. 
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The  single  value  of  s  1.1  was  used  in  this  graph,  and  the  effect  of 
a  change  in  phase  upon  the  current  was  examined.  As  can  be  seen,  the 
phase  of  the  output  gap  voltage  has  a  marked  effect.  For  a  phase  of 
ir/2  ,  the  current  is  very  flat  over  a  large  portion  of  the  curve.  How¬ 
ever,  increasing  the  phase  to  ir  produces  a  very  large  current  peak 
at  (a)tj  s  3.8  radians  with  a  relatively  low  value  elsewhere.  The  other 
two  values  of  phase  given  are  intermediate  between  these  two  extremes. 
A  more  detailed  graph  for  the  phases  between  ir/2  and  ir  shows  the 
flattened  portion  of  the  curve  rising  continuously  from  «/2  to  the  large 
peak  observed  at  ir  .  One  would  suspect  that  these  two  extremes  cor- 
respond  to  maximum  power  extraction^  from  the  beam  for  the  phase  of 
ir/2  ,  and  to  maximum  addition  of  power  to  the  beam  to  produce  the 
highly  bunched  condition  for  the  phase  of  «  .  In  order  to  check  this, 
the  average  power  per  period  versus  the  phase  was  calculated  as  out¬ 
lined  in  the  previous  discussion.  This  is  given  in  Figure  7. 

The  original  curve  for  all  values  of  6  was  made  with  Oj  s  1.1  . 
As  was  suspected,  the  beam  loses  the  maximum  amount  of  power  in  the 
vicinity  of  ir/2  ;  however,  it  gains  the  most  power  at  0  s  7ir/4 .  This 
is  far  from  the  point  of  greatest  bunching,  ir  .  Indeed,  at  8  =  v  ,  the 
beam  is  still  losing  power  to  the  gap.  This  would  seem  to  indicate  that 
a  condition  of  maximum  harmonic  content  of  the  beam  current  could  be 
achieved  while  power  is  still  being  extracted  from  the  beam. 

The  point  of  maximum  power  extraction  is  closer  to  8  s  5ir/8  . 

A  check  of  this  point  for  a  different  value  of  (a^  =  .7)  was  made  and 
is  also  presented  on  Figure  7.  It  seems  that  for  a  large  change  in  , 
the  minimiun  shifts  slightly,  being  closer  to  ir/2  for  larger  values  of 
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.  However,  over  the  region  of  of  interest,  the  change  in  the  phase 
6  for  maximum  power  extraction  is  very  small.  The  maximum  kinetic 
power  extraction  takes  place  when  the  phase  0  is  in  the  interval  (ir/2,  3ir/4). 
We  therefore  examine  the  exit  current,  exit  velocity,  and  exit  charge  den¬ 
sity  as  functions  of  time  for  these  values  of  phase  in  Figures  8-1 6,  since 
klystron  amplifiers  are  normally  operated  under  these  conditions. 

For  a  particular  ,  all  three  velocity  graphs  are  the  same,  but 

for  each  different  0  they  are  displaced  from  each  qther.  This  is  to  be 

expected  from  the  form  of  Equations  (79)  and  (80),  where  the  quantity 

(<i)t  4-  0)  always  appears,  instead  of  wt  or  0  separately.  Therefore, 
c  c 

a  variation  in  phase  merely  causes  a  change  in  phase  of  the  output  velocity. 

The  current  graphs  show  remarkable  changes,  but  have  some 
similarities.  For  0  =  ir/2  ,  the  curve  for  =  1.1  is  the  same  as  is 
given  in  Figure  6.  But,  Figure  8  shows  an  increase  in  the  current  with 
decreasing  a^  ,  for  electrons  in  the  region  2.7  l«t^<5  .  Figure  9 
shows  that  these  electrons  have  the  lowest  velocity.  Therefore,  it  is 
very  probable  that  the  electrons  with  exit  times  in  this  region  gave 
power  to  the  external  load  and  were  slowed  down.  With  less  retarding 
field  evidenced  by  a  lower  a^  ,  the  electron  velocity  would  not  be  de¬ 
creased  as  much,  and  the  current  would  be  greater,  since  it  is  the  prod¬ 
uct  of  velocity  and  charge  density.  The  charge  density  function  does 
vary  with  time  towards  the  end  of  this  interval,  as  can  be  seen  from 
Figuire  10.  This  is  evidence  of  a  further  bunching  effect  coming  into 
play.  It  is  interesting  to  note  that  the  limiting  case  of  smaller  a^  must 
be  a  value  of  one  for  the  velocity  at  all  time,  and  a  current  equal  to  the 
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Figure  10.  Output  Gap  Exit  Charge  Density  versus  Exit  Time; 
0  =  ir/2  . 


- 


Figure  11.  Output  Gap  Exit  Current  versus  Exit  Time;  9  =  5ir/8  . 


Figure  12.  Output  Gap  Exit  Velocity  versus  Exit  Time;  0  s  5ir/8  . 
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e  14.  Output  Gap  Exit  Current  versus  Exit  Time;  6  =  3^/4. 


;ure  15.  Output  Gap  Exit  Velocity  versus  Exit  Time;  6  s  3ir/4 


Figure  16.  Output  Gap  Exit  Charge  Density  versus  Exit  Time; 
e  =  3Tr/4 . 
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entrance  current. 


The  mechaniam  just  mentioned  for  power  extraction  at  8  s  tr/Z 
can  also  be  observed  at  6  =  5ir/8  and  9  =  3it/4  .  In  these  cases,  the 
region  of  minimum  velocity  is  displaced,  but  the  current  variations 
are  in  general  displaced  in  the  same  direction  by  approximately  the 
same  magnitude,  so  that  the  power  extraction  occurs  during  the  in¬ 
terval  of  the  large  current  bunch.  It  is  worth  noting  that  the  charge - 
density  functions  show  a  huge  peak  during  that  interval,  but  since  the 
velocity  is  larger  outside  the  interval,  we  can  conclude  that  the  vast 
majority  of  the  electrons  give  up  their  energy  during  that  interval  and 
that  the  current  peaks  outside  that  interval  are  caused  by  small  numbers 
of  relatively  high-velocity  electrons.  The  magnitude  of  the  current  in 
the  region  5.2  <  eft^  <  8.5  seems  to  be  relatively  insensitive  to  changes 
in  for  a  given  0  .  However,  an  increase  in  0  gives  a  slight  de¬ 
crease  in  the  current  in  these  two  peaks,  since  the  current  in  the  largest 
peak  increases,  and  the  current  must  be  conserved  over  a  period  in 
time. 

The  two  smaller  peaks  in  the  current  functions  in  the  region 
5.2  <  wt^  51  8.5  are  caused  by  high-velocity  electrons,  as  can  be  seen 
from  the  velocity  curves.  They  have  a  relatively  low  charge  density, 
indicating  that  the  number  of  the  high-velocity  electrons  is  small.  This 
is  evident  from  the  charge -density  curves.  The  velocities  of  these  high- 
velocity  electrons  are  independent  of  phase  but  are  directly  proportional 
to  voltage  magnitude  .  The  number  of  particles  with  the  high  vel¬ 
ocities  is  not  a  function  of  voltage  magnitude,  but  does  decrease  some¬ 
what  with  increasing  phase.  The  position  in  time  of  these  high-velocity 
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bunches  does  not  vary  with  magnitude  or  phase  of  the  voltage. 

As  a  side  effect,  one  might  point  out  that  one  would  expect,  for 
a  value  of  s  1.1  ,  a  maximum  value  of  velocity,  which  corresponds 
to  approximately  times  the  d*c  beam  voltage,  for  the  electrons 
that  gained  the  most  energy  in  passing  through  the  gap.  This  is  indeed 
the  case,  as  can  be  seen  on  Figures  9,  12,  and  15. 

The  results  seem  to  indicate  that  the  X-rays  are  produced  by 
relatively  few  electrons  at  high  velocity,  which  appear  during  the  latter 
part  of  the  period.  Therefore  a  large  d-c  retarding  field  would  not  be 
practical,  for  if  it  were  effective  in  slowing  down  the  few  fast  electrons, 
it  would  cause  the  very  large  number  of  slow  electrons  to  be  repelled 
back  into  the  gap,  where  they  would  interact  with  the  later  electrons. 
There  is  even  the  possibility  of  reflex  oscillations  under  these  conditions. 

It  seems  very  possible,  however,  that  the  often  suggested  r-f 
structure  placed  ahead  of  the  collector  could  actually  be  used.  It 
appears  that  it  would  have  to  be  driven  by  an  outside  voltage  source, 
possibly  at  twice  the  fundamental  frequency  w  in  order  to  affect  both 
fast  bunches.  Its  effect  on  the  slow  large  bunch  would  probably  be  only 
to  speed  up  some  electrons  to  the  order  of  magnitude  of  the  beam  voltage, 
and  would  not  necessarily  cause  any  electrons  to  turn  back,  although 
there  would  be  a  tendency  for  some  electrons  entering  in  the  large  slow 
bunch  to  do  so.  Possibly  a  cavity  tuned  to  twice  the  fundamental  fre¬ 
quency  could  extract  power  from  the  fast  electrons,  but  only  a  more 
detailed  analysis  could  confirm  this.  At  any  rate,  the  coherent  charac¬ 
ter  of  the  beam  current  and  velocity  at  the  exit  plane  of  the  output  gap 
leads  one  to  believe  that  a  r-f  structure  placed  before  the  collector  with 
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the  proper  voltage  magnitude  and  phase  maintained  upon  it  could  be 
singularly  effective  in  slowing  the  two  small  bunches  of  fast  electrons 
without  seriously  increasing  the  velocity  of  the  rest  of  the  beam. 

V.  CONCLUSIONS  AND  RECOMMENDATIONS 

In  the  first  half  of  the  preceding  analysis,  the  ballistic  theory 
of  an  electron  beam  in  a  gridded  gap  was  developed  for  large  signals. 

The  limits  of  validity  of  the  expressions  for  velocity  and  current  were 
established.  With  the  neglect  of  space  charge,  this  analysis  provides 
a  rigorous  analysis  of  the  electron  beam  in  a  gridded  gap  region.  An 
experimental  verification  of  this  basic  gap  theory  could  be  undertaken, 
using  a  set  of  parameters  designed  to  realize  more  closely  the  assump> 
tion  of  an  infinite  beam  and  no  space -charge  effects. 

This  gap  analysis  was  applied  in  the  output  gap,  combined  with  the 
best  available  large -signal  theory  using  space -charge  waves  for  the  drift 
region.  The  exit  current  and  exit  velocity  of  the  beam  at  the  end  of  the 
output  gap  was  found  for  several  values  of  the  magnitude  and  phase  of  the 
output  voltage.  The  change  in  average  kinetic  power  of  the  beam  across 
the  gap  was  also  calculated,  and  it  was  concluded  that  under  conditions  of 
normal  operation  of  two -cavity  klystron  amplifiers,  there  are  two  small 
bunches  of  relatively  high-velocity  electrons  in  the  beam.  The  positions 
of  the  bunches  is  constant  for  several  voltages  and  phases,  and  the  am¬ 
plitude  of  the  current  bunches  decreases  with  decreasing  voltage. 

As  a  logical  continuation  of  this  theoretical  analysis,  some  ex¬ 
perimental  investigation  of  the  possibility  of  slowing  down  the  fast 
electrons  with  a  r-f  structure  should  be  undertaken.  It  seems  very 
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probable  that  such  a  structure  could  effectively  slow  the  fast  electrons 
and  prevent  production  of  hard  X-rays  at  the  collector. 


APPENDIX 


A.  DERIVATION  OF  POWER  SERIES  INVERSION 

The  derivation  of  Equation  (41)  has  been  quoted  by  many  authors, 
but  originated  with  Lagrange. We  are  following,  the  method  given  by 
Tisserand^^  (pp  262-269).  Consider  the  equation 


Z  -  a  -  of(Z)  =  0  ,  (A.  1) 

in  which  a,  a,  and  Z  designate  real  or  imaginary  quantities.  Con¬ 
sider  also  a  closed  contour  S  such  that  one  has,  on  all  its  points, 

mod  <  1  .  (A.  2) 


We  further  suppose  the  function  f(Z)  to  be  holomorphic  on  the  interior 
of  S. 


It  has  been  shown  that  the  Equation  (A.  1)  admits  one  and  only 
one  solution  Z  in  the  interior  of  S  and,  in  designating  by  ir(Z) 
any  holomorphic  function  of  that  root,  we  can  develop  ir(Z)  in  a  con¬ 
verging  series  in  powers  of  a  : 

n  s  oo 


n  +  1 


ir(Z)  =  ir(a)  + 


TTT. 


ill 


n  =  o 


(A.  3) 


We  can  take  for  Equation  (A.  1)  the  Kepler  equation, 


u  -  ^  €  sinu  =  0  ,  (A.  4) 

in  which  we  assume  ^  and  c  to  be  real.  According  to  the  preceding 
statement,  if  one  can  find  a  closed  contour  S  on  all  points  of  which 
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(A.  5) 


u  u  -  ; 


then  Equation  (A.  4)  admits  one  and  only  one  root  u  in  the  interior  of 
S  ,  and  one  can  develop  the  converging  series, 

n  =  00 


ff{u)  =  IT  (C)  + 


n+ 1  .n 
c  d 


(n  +  i).'  j^n 


ir'(4)  sin 


n  +  1^ 


n  =  o 


(A.  6) 


In  making 


and 


ir(u)  =  u  , 


Tr(u)  =  cos  u  . 


(A.  7) 


(A.  8) 


successively  in  Equation  (A.  6),  it  follows  that 


u  s  ^  4  c  sin|^  -f  . .  . 


and 


m  ,m  -  1 
at 


r  2, 

cos  u  =  cos  t  *  <  Bin  t  *  •  * 


(  A.  9) 


m -  1  ,m  -  2 
<  d  .  m . 

(m  -  iji  j^m  -  2  t  ■  •  • 


(A.  10) 


But,  one  also  has  these  known  formulas.  For  m  even, 
m 

sin™t  =  I  |cos  mt  -  ^  cos  (m-2)t  f  cos  (m  -  4)t  -  .  .  • 


(A.  11) 
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and  for  m  odd, 


m-1 


™  •in  (m-2)  4  ®i>^  {ni-4)4  -  .  .  .  ] 


(A.  12) 


From  these  formulas  one  can  infer  easily  that 

d4^  2' 


.  m,  1  r  m-1  ,  m  ,  ,vm-l  .  , 

— ^  am  4  =  J  [  sinm4  - -j- (m-2)  8in(m-2)4 


+  Hjl2^!l(m-4)”'‘l  8in(m-4)4  -  .  .  .  ]  (A.  13) 

This  formula  holds  for  the  two  cases  of  m  even  or  odd.  It  is  terminated 
at  the  term  in  sin2  4  if  ni  is  even,  and  at  the  term  in  sin4  if  m  is  odd. 

Using  this  relation  in  Equation  (A.  9)  and  (A.  10),  we  have  the  final 

rfssult: 

'  2  m  ,  ,  m  m  \ 

us4'*’<  •inC  •in24  +  -  •  •  - |^m"'"^sinm4  --j-(m-2)“'"*sin(m-2)4 

2  mi 

+  (m-4)"*‘^  •in(m-4)4  -  .  .  .  ]  ,+....  (A.  14) 

B.  CONVERGENCE  OF  POWER  SERIES  INVERSION 
1 8 

Laplace  was  the  first  to  show  the  limit  on  the  convergence  of 
the  power  series.  His  method  was  very  complicated,  and  so  we  shall 
give  a  simpler,  shorter  method  ascribed  to  Rouche^  that  leads  to  the 
same  result. 

Let  A  be  the  point  of  the  x  axis  where  the  abcissa  is  4  • 
m  be  the  point  where  the  abcissa  is  4  u  .  We  take  for  the  contour  S 
a  circle  of  radius  p  drawn  from  A  as  center.  Let  us  make  the  point 
M  move  on  this  circle,  and  designate  the  angle  between  the  radius  AM 
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and  the  x  axis  by  <j>  .  Let  us  adopt  Tisserand's  notation  of  i  = 
and  £  =  Naperian  base  =  2.71828.  .  .  .  Then  we  have,  from  these  defin¬ 
itions,  , 

u  =  C  +  .  (A.  15) 

Condition  (A.  5)  becomes 

^  mod  sin  (4  +  p  <  1  .  (A.  16) 

Let  F(p)  be  the  maximum  value  of  the  modulus  of  sin  (4  p  S*^^) 
when  ^  varies  from  zero  to  2ir  and  4  takes  on  all  possible  real  values. 
If  the  value  for  c  is  such  that 

‘  ^  ^ 

the  condition  given  in  Equation  (A.  16)  will  be  satisfied  and  the  series  will 
converge  for  all  values  of  4  •  one  could  determine  the  value  of  p^ 
of  the  radius  p  so  that  the  expression  p/F(p)  is  as  large  as  possible, 
and  then  one  makes  <2  ~  series  would  converge  for 

€  <  €  j 

One  must  first  find  the  form  of  F(p).  In  order  to  do  this,  we 
compute  the  square  of  the  modulus  of  sin  (4  pE^^)  •  It  i*  equal  to 

sin  (4  p  E^^)  sin  (4 -t  p  E~^^)  =  sin  (4 'f  p  cos  ip  sin  ^)  • 

2  2 

'  sin  (4 p  cos  <^  -  i  p  sin  ^)  =  cos  (ip  sin  ^) -cos  (4  p  cos  ^)  .  (A.  19) 


-66- 


One  therefore  has 


mod  sin 


(l*pE'*)=^j  l£ - - j 


cos^  (C  +  p  cos  <^)  .  (A.  20) 


The  miniirium  of  p/F(p)  occurs  when  Equation  (A.  20)  is  at  its 
maximum,  for  a  given  value  of  p  .  This  takes  place  when 

l^pein^^  ^-psin,^j 

is  as  large  as  possible,  at  sin^  =  1  .  Then,  one  will  have  at  the  same  time 


cos  (C  -t  p  cos  ^)  s  0 
cos  C  =  0  , 

t  . 


(A.  21) 
(A.  22) 
(A.  23) 


For  this  case,  it  will  happen  that 

P,p, .  £4^ 


(A.  24) 


In  order  to  find  the  maximum  of  the  expression. 


IJP  +  E'P 


(A.  25) 


we  set  its  derivative  equal  to  zero: 


eP(p  -  1)  -  E"P<p  +  1)  =  0 


(A.  26) 


The  first  term  in  this  equation  takes  values  of  opposite  sign  for 
p  =  1  and  p  =  2  ;  its  derivative  p  (E^  E~^)  is  always  positive.  There¬ 

fore  this  equation  has  one  and  only  one  root  p^  ,  which  is  positive.  One 
finds  from  numerical  computations  that 
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(A.  27) 


Pj  =  1.9967. . . 


so  that  the  limit  of  convergence  is 

2 


‘l  =-p 


’1 


=  0.6627  . . . 


+  E 


'1 


(A.  28) 


This  is  the  limit  of  convergence  of  the  power  series. 
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ABSTRACT 


In  this  report,  the  design  of  the  coaxial  cavities  used  by  high- 
power  klystrons  is  specified  so  as  to  minimize  the  spurious  output 
power.  A  detailed  theoretical  analysis  of  a  general  doubly  re-entrant 
grid-less  gap  coaxial  cavity  is  made  which  enables  one  to  calculate 
the  fundamental  and  first  three  higher-order  resonances.  It  is  shown 
that  the  fiuidamental  resonance  may  be  held  fixed  while  the  first  few 
higher -order  resonances  are  adjusted  in  frequency  so  as  to  avoid  the 
second  and  third  harmonics  of  the  signal  frequency.  In  this  way,  the 
output  cavity  of  a  klystron  may  be  optimized  with  respect  to  minimum 
spurious  output  power,  which  occurs  largely  at  the  harmonics.  Experi¬ 
mental  measurements  made  on  a  symmetrical  cavity,  i.  e. ,  no  large 
output  port  or  loop,  indicate  accuracies  of  a  few  per  cent  for  the  cal¬ 
culated  resonances.  Limited  design  curves  are  presented  which  have 
a  universal  application  and  equations  are  given  so  that  additional  curves 
may  be  calculated. 

The  second  section  of  this  report  treats  the  general  theory  of 
the  perturbation  of  resonant  cavities  from  a  physical  viewpoint.  The 
application  of  this  theory  can  aid  the  large  signal  solution  of  the  beam 
dynamics  in  a  klystron,  which  would  lead  to  the  effects  of  the  interaction 
phenomena  on  the  spurious  output.  In  particular,  formulas  are  given 
for  the  calculation  of  absolute  values  of  the  fields  from  perturbation 
measurements.  It  is  shown  that  the  perturbing  object  must  be  on  the 
order  of  two  or  three  hundredths  of  a  wavelength  in  size  for  these 
formulas  to  be  accurate. 
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I.  INTRODUCTION 


Within  the  past  few  years,  engineers  have  become  increasingly 
aware  of  the  serious  problems  presented  by  the  spurious  output  of  high- 
power  microwave  systems,  and  a  recent  paper  by  R.  D.  Campbell^  has 
shown  the  severity  of  this  problem  as  related  to  interference  caused  by 
radar.  The  paper  cites  experiments  in  which  measurements  were  made 
on  several  radars  to  determine  their  power  output  spectrum.  All  of  the 
radars  showed  "very  substantial  magnitudes"  of  power  radiated  at 
frequencies  other  than  the  signal  frequency.  The  radars  tested  employed 
both  magnetrons  and  klystrons  as  the  microwave  power  source. 

There  are  two  broad  methods  of  approaching  this  spurious  ou^ut 
problem.  The  first  is  to  suppress  the  unwanted  output  by  some  means 
external  to  the  high>power  microwave  tubes,  and  the  second  is  to  improve 
the  basic  design  of  these  tubes  to  prevent  the  generation  of  spurious 
signals.  The  first  method  can  be  used  to  improve  existing  installations 
while  the  second  is  desirable  for  future  installations.  This  thesis  is 
devoted  to  the  second  of  these  methods  and  deals  with  high*power  klystrons. 

As  related  to  spurious  outputs  there  are  two  approaches  to 
klystron  design: 

(1)  active  -  the  electron  beam  and  interaction  phenomena, 

(2)  passive  -  the  microwave  circuitry  of  the  tube. 

These  two  subdivisions  are  by  no  means  unrelated  to  each  other.  The 
large  signal  operation  of  klystrons  produces  nonlinearities  which  are 
the  basis  for  the  development  of  power  at  frequencies  other  than  the 
signal  frequency.  Then  for  this  power  to  become  spurious  output^ 
the  microwave  circuit  must  couple  to  it.  Part  of  the  design  of  klystrons 
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is  therefore  associated  with  the  optimum  design  of  the  output  cavity, 
optimum  in  the  sense  that  only  the  fundamental  signal  frequency  is 
coupled  out.  This  means  that  the  pass  bands,  or  essentially  the  resonant 
modes  and  frequencies,  must  be  delineated,  for  it  is  only  at  or  near 
these  frequencies  that  the  cavity  can  have  any  output. 

The  microwave  circuitry  of  klystron  design  is  undertaken  in 

the  first  part  of  this  thesis.  The  analysis  is  directed  at  calculating  the 

resonances  of  the  usual  type  of  klystron  cavity  so  that  the  results  may 

be  put  into  a  form  usable  in  design  work.  The  general  case  of  the  doubly 

re-entrant,  gridless -gap,  coaxial -type  cavity  is  treated.  To  date,  no 

known  solution  of  this  type  of  cavity  exists,  although  the  singly  re-entrant, 

gridded-gap  cavity  has  been  treated  numerically  at  the  fundamental 

2.  3 

resonance  with  very  good  results. 

After  the  analysis  is  completed  and  verified  experimentally, 
design  criteria  are  outlined  that  will  permit  the  design  of  a  cavity  whose 
higher  order  resonances  are  removed  from  the  harmonics  of  the  signal 
frequency.  The  design  has  been  developed  to  eliminate  harmonics  up 
to  and  including  the  fourth.  This  particular  criterion  based  on  the 
elimination  of  resonances  at  the  signal  harmonics  has  been  used  because 
a  large  part  of  the  spurious  output  encountered  in  practice  occurs  at 
these  harmonics. 

The  second  part  of  this  thesis  is  concerned  with  the  general 
theory  of  the  perturbation  of  resonant  cavities.  This  theory  is  useful 
in  that  its  application  indirectly  aids  the  large  signal  solution  of  the 
beam  dynamics  in  a  klystron,  a  solution  that  is  a  prerequisite  to  the 
study  of  the  interaction  phenomena  as  related  to  spurious  outputs. 

In  particular  perturbation  techniques  can  be  used  to  determine  relative 
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field  distributions.  This  experimental  fact  is  well  known,  but  perhaps 
not  so  well  known  is  precisely  how  various  parameters  affect  the  amount 
and  type  of  perturbation.  Some  effort  has  been  made  here  to  clarify 
this  and  to  show  the  limitations  on  absolute  measurements. 
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II.  RESONANT  CAVITY  THEORY 


It  is  always  a  difficult  problem  to  treat  a  high-frequency  circuit 
by  means  of  the  lumped -circuit  type  of  analog.  Invariably  several  ways 
of  defining  the  lumped -circuit  parameters  arise,  the  three  usual  ones 
being  the  following: 

(1)  definitions  based  on  geometry, 

(2)  definitions  based  on  stored  energy, 

(3)  definitions  based  on  "equivalent"  admittances. 

The  model  used  here  to  represent  the  microwave  cavity  is  the 
one  based  on  "equivalent"  admittances.  Essentially,  this  will  involve 
finding  the  electromagnetic  fields  as  a  series  summation  for  different 
regions  of  the  cavity  and  then  matching  these  solutions  at  the  boundaries 
between  the  regions.  Then  as  an  aid  in  using  these  field  expansions, 
admittance  functions  are  defined  so  that  the  absolute  value  of  the  total 
admittance  function  is  a  minimum  at  a  cavity  resonance.  This  treat¬ 
ment  idealizes  the  cavity  structure  to  be  lossless  so  that  these  minima 
correspond  to  true  zeros  of  the  total  admittance  function.  This  means 
that  the  model  consists  of  a  sequence  of  ideal,  parallel,  resonant,  lumped 
circuits.  In  fact  some  of  the  properties  of  these  two  resonant  devices 
are  very  similar. 

The  cavity  geometry  of  particular  interest  is  the  doubly  re¬ 
entrant  gridless -gap  type,  which  is  the  one  usually  associated  with 
high -power  klystrons.  A  typical  cavity  with  its  co-ordinate  system 
is  shown  in  Figure  1 . 
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Figure  1.  Co-ordinate  System  for  Cavities. 
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For  finding  the  field  expansions,  the  cavity  may  be  divided  into 
two  regions: 

1)  The  region  external  to  the  drift  tubes,  henceforth  called 
region  1.  The  range  of  r  is  r  i  r  i  r,  . 

2)  The  region  inside  the  outer  radius  of  the  drift  tubes,  hence¬ 
forth  called  region  2.  The  range  of  r  is  0  ^  r  ^  r  . 

The  physical  structure  of  the  cavity  lends  itself  well  to  a  radial-line 
treatment,  which  is  the  reason  for  the  two-part  division  above. 

Region  1  corresponds  to  a  radial  line  which  is  short-circuited  at  the 
outer  radius,  while  region  2  involves  the  "capacitive"  loading  placed 
on  this  line. 

The  usual  fundamental  mode  of  oscillation  of  this  cavity  is  the 
coaxial  line  mode,  where  the  subscripts  indicate  the  spatial 

variations  (of  the  electromagnetic  fields)  in  the  ^ ,  r  ,  and  z  directions 
respectively  and  the  TM  refers  to  the  transverse -magnetic  type  of 
mode,  that  is,  one  having  no  component  of  the  magnetic  field  in  the  z 
direction.  When  used  as  an  output  device  this  type  of  cavity  is  excited 
by  means  of  a  z-directed  electric  field  located  in  the  gap  region.  This 
type  of  excitation  can  only  couple  to  modes  having  a  z  component  of 
electric  field,  so  that  the  higher  order  resonant  modes  will  also  be  of 
the  TM  type.  Since  the  first  region  of  this  cavity  is  treated  as  a 
radial  line,  the  above  nomenclature  may  be  confusing  in  that  it  applies 
to  the  transverse  area  dimensions  r  ,  ^  while  the  radial  line  has 
transverse  area  dimensions  z.  It  would  therefore  be  convenient 
to  develop  a  nomenclature  that  would  apply  to  the  radial  line.  Since 
a  nomenclature  for  the  radial  line  has  not  been  standarized  yet,  the 
following  somewhat  arbitrary  definitions  are  adopted: 
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TEM  mode  —  the  mode  which  has  electromagnetic  field 

components  only  in  the  area  transverse  to 

the  direction  of  propagation,  i.  e.  ,  only 

E  and  H.  exist, 
z  o 

Hmnp™°^^  —  H  =0,  subscripts  refer  to  the  spatial  variations 
of  the  electromagnetic  field  in  the  z  ,  6  i  and  r 
directions  respectively. 

E  mode  —  E  =0,  these  correspond  to  the  coaxial  line  TE 
mnp  z 

modes  and  will  not  be  used. 

Thus  it  is  apparent  that  the  TEM  mode  corresponds  to  the  coaxial  line 
TM  mode,  while  the  H  modes  correspond  to  the  coaxial  line 

010  *nnP 

TM  modes.  From  this  point  on  all  mode  designations  are  in  terms 
npm 

of  the  radial -line  nomenclature  unless  otherwise  specified. 

A.  THEORETICAL  ANALYSIS 

Maxwell's  equations  may  be  solved  in  the  cavity  subject  to  the 
boundary  conditions  to  give  the  electropnagnetic  field  expansions  in  the 
form  of  series  summations.  From  these  field  expansions,  voltages 
and  currents  may  be  defined  and  the  admittance  function  associated  with 
the  cavity  resonances  can  then  be  obtained.  One  must,  of  course  be 
consistent  in  these  definitions;  i.  e. ,  all  voltages  and  currents  must  be 
defined  at  one  particular  reference  "plane"  or,  as  it  turns  out  in  this 
case,  at  one  value  of  radius.  The  reference  "plane"  used  here  is  the 
cylinder  located  at  the  outer  radius  of  the  drift  tubes.  For  the  principal 
part  of  this  analysis  symmetry  in  the  ^  direction  is  assumed,  since 
this  is  how  these  cavities  are  theoretically  designed  to  operats.  How¬ 
ever,  it  would  also  be  of  interest  to  include  a  treatment  of  nonsymmetrical 
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operation,  eince  perturbations  such  as  struvjtural  inmper factions,  and 
the  tuners  and  antenna  one  invariably  inserts  in  these  cavities  can 
easily  excite  nonsymmetrical  modes.  A  first  approximation  to  this 
nonsymmetrical  operation  is  made  in  a  later  section,  where  a  rather 
artificial  technique  is  used  to  introduce  the  nonsymmetrical  modes. 

The  assumptions  made  in  the  following  derivation  are  quite 
usual,  and  are  good  approximations  for  the  practical  cavity  which 
is  either  lightly  loaded  or  symmetrically  loaded.  The  conductivity 
of  the  metal  walls  of  the  cavity  is  assumed  infinite  and,  most 
important,  the  gap  edges  are  assumed  to  be  square.  The  one 
further  assumption  (as  mentioned  above)  is  that  perfect  mechanical 
symmetry  exists  in  the  direction. 

1 .  Cavity  Admittance  of  Region  1 

Since  this  region  is  to  be  treated  as  a  short-circuited  radial 
line,  one  need  not  begin  with  Maixwell's  basic  equations;  instead  a 
modal  analysis  can  be  used.  These  radial -line  modes  themselves 
have  been  derived  from  Maxwell's  equations,  so  there  is  no  loss 
of  generality  in  applying  them  directly. 

The  general  expressions  for  the  field  components  in  a  radial - 
line  cavity  for  the  H  type  modes  are  as  follows: 
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and  p^  is  the  radial  phase  constant  (wave  number), 
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Actually,  in  Equation  (1)  there  is  a  choice  of  cosn^  or  sinn4>  for  E  . 

z 

The  only  time  these  occur  simultaneously  (termed  mode  splitting)  is 
when  the  circular  cavity  is  slightly  out  of  round  or  some  other  perturbation 
in  the  ^  direction  is  present. 

For  symmetrical  modes,  when  n  =  0,  Equation  (1)  becomes 

_  (1)  jwt 

®  coB^  z  e*' 

z  mot  mOt  z 


H  =  0 
z 

(P  r)  +  B  (P  r)D  cosp  z 

9  r^Lmit  ml  t  z 

H  =  0 
r 

These  equations  now  specify  modes  of  the  H  .  type.  The  particular 

mOp 

mode  with  m  =  0  has  E  =0  and  is  the  TEM  mode.  At  this  point 

r 

the  meaning  of  "P  variations  in  r"  may  not  appear  clear.  This  corres 
ponds  to  having  P  variations  in  the  standing  wave  pattern  in  the 
longitudinal  direction  of  a  rectangular  cavity  which  is  resonating. 

Equation  (4)  can  now  be  applied  to  the  cavity.  Figure  2  shows 
a  section  of  the  cavity  and  its  dimensions.  To  find  the  coefficients 
A^  and  in  Equation  (4),  a  field  may  be  assumed  at  the  gap. 
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Figure  2.  Cross-Section  of  Cylindrical  Cavity  Showing  Dimensions. 


From  this  same  assumed  field  a  voltage  at  the  gap  can  be  defined. 

Then  knowing  H  ,  a  current  may  be  defined  which  can  be  used  to 
give  us  the  admittance  associated  with  the  short-circuited  radial  line. 
For  simplicity  the  actual  field  that  will  exist  at  the  gap  is  approximated 
by  a  constant  as  follows: 
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Figure  3  shows  the  relationship  between  this  assumed  field  and  the  true 
field. 


For  region  1  (in  Figure  2,  the  region  between  r  and  r,  ),  the  field 
coefficients  can  be  found  by  knowing  £  and  H  at  a  particular  value 
of  radius,  or  by  knowing  one  of  these  at  two  values  of  radius.  Here 
tl^e  latter  is  chosen  since  £  at  r  is  known  by  £quation  (5),  and 

Z  ft 
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since  E  at  r  is  zero  to  match  the  boundary  conditions.  These 
z  b 

conditions  and  Equation  (4)  give: 


=  E  e'^"*  =\  CA  J  (P.i*  )  +B  H  r  G  cos 

0  /  mo^a  mO  ta  / 

(6) 


ut 


m=l 


and 


(1) 


/miTz\ 


jut 


=  0  .  >  CA^J„(P,r^)  .  (p^r^)3co.  . 

m=^ 

«# 

(7) 

which  now  include  a  sum  over  all  the  possible  modes  allowed.  Because 
of  the  orthogonal  properties  of  cos  ^  >  msO,  1,2...,  over 
the  transverse  area  both  of  these  equations  can  be  multiplied  by  a 
particular  cos  ^  ,  integrated  over  the  transverse  area,  and  have 
one  term  left  on  the  right  hand  side.  Thus  Equation  (7)  gives 


m 


(m  =  0  ,  1  ,  2  ,  .  .  .  ) 


and  Equation  (6)  becomes 
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which  hold*  for  m  0  .  The  result  is  two  equations  with  two  unknowns 

A  and  B  Solving  for  B  gives  the  following: 
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Equations  (8),  (9)«  and  (10)  give  the  coefficients  of  Equation  (4) 
so  that  the  complete  fields  in  the  cavity  are  now  known.  The  voltage 
at  the  gap  follows  immediately  from  Equation  (5)  as  follows: 
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The  current  induced  at  the  gap  may  be  found  from  an  integration  of  H 


4> 


over  the  gap  area  as  follows: 
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Substituting  for  from  Equation  (4)  and  inserting  the  series  coefficients 

just  calculated,  gives 
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From  Equation  (3)  p^  becomes 


However,  the  cutoff  frequency  for  the  mth  mode  is 


w_ 


mrr 


m 


L  V  uc 


(14) 


so  that  P  becomes 
t 


(15) 


In  calculating  the  admittance,  it  is  convenient  to  normalize  with 
respect  to 


Ztr  r. 


(16) 


The  admittance  is  found  by  using  Equations  (11),  (13),  (15)  and  (l6) 
as  follows: 


-15- 


In  the  form  given  this  equation  applies  strictly  to  the  "propagating" 
modes,  where  a  "propagating"  mode  is  defined  as  having  u  >  a»^  .  How¬ 
ever,  there  will  certainly  be  modes  for  which  u  <  is  true  and  these 
modes  are  called  "cutoff.  "  Equation  (17)  can  be  written  to  apply  directly 
to  all  modes  by  using  the  proper  value  of  the  radial  phase  constant 
for  the  cutoff  modes.  This  value  is 


P  cutoff  modes 


(18) 


Using  the  modified  Bessel  functions 


I  (x)  =  j  J  (jx) 
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n 


and  making  the  substitution 


(20) 


where  the  negative  sign  of  Equation  (18)  is  now  taken  into  account  by 
realizing  that 


J^(-x)  =  -J^(x)  , 


(21) 


one  can  find  the  admittance  function  which  accounts  for  all  the  modes 
in  region  1 : 
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In  this  equation  the  first  K  modes  are  "propagating"  while  all  the 
higher  modes  are  cutoff. 

Before  Equation  (22)  is  actually  used,  it  must  converge.  The 
only  term  that  requires  checking  is  the  third,  since  this  is  the  only 
infinite  summation.  From  Equation  (20), 


lim  6  '  =  lim 
m - ►00  m - ►  00 


mu. 


u 


(23) 


where  u^  is  the  first  cutoff  frequency  and 


u  =  mu  , 
m  1 

which  follows  from  Equation  (14).  It  is  also  known  that  as  the  argument 
becomes  large,  the  modified  Bessel  functions  become 


I„(X) 


X 


e 


K  (x) - ^7 

n  ^ 


-x 

e 


x - ►»oo 


(24) 


With  the  aid  of  Equations  (23)  and  (24)  it  follows  that: 


lim 

m— ♦-oo 


.2p 

mu| 

Ir 

-r  ) 

a 

1  +  e 

U) 

'^b 

1  -  e 

.2p 

mu^ 

Ui 

-'a> 

(25) 


-19- 


The  term  — -  in  the  denominator  is  conr  non  to  all  the  terms  from 
u 

here  on  and  may  be  dropped.  Then,  since 


and  is  known  to  converge,  the  series  in  Equation  (22)  also  con- 

m 

verges,  provided  w  /  w  .  Unfortunately  Equation  (22)  is  in  a  form 

m 

such  that  the  usual  "convergence  techniques"  cannot  be  applied  con¬ 
veniently;  that  is,  the  method  of  separating  a  low  frequency  admittance 
from  the  infinite  summation  rather  than  adding  contributions  to  this 
admittance  from  the  various  modes  is  not  applicable  in  this  case.  This 
means  that  a  calculation  of  the  effect  of  the  cutoff  modes  must  certainly 
be  delegated  to  a  digital  computer  since  experience  has  shown  that 
twenty  or  thirty  terms  of  this  series  are  required  for  any  reasonable 
accuracy. 


2.  Cavity  Admittance  of  Region  2 

Region  2,  as  Figure  2  shows,  can  be  divided  into  two  subregions. 

The  first  is  the  region  directly  between  the  ends  of  the  drift  tubes,  for 

the  drift  tubes  do  have  finite  thickness  ind  the  second  is  inside  the 

drift  tubes.  To  be  strictly  correct,  one  should  find  the  admittance  of 

the  inner  region  at  the  inside  of  the  gap  edge  and,  by  the  use  of  the 

general  radial  line  equations  transform  this  admittance  from  r.  to 

r  Also,  the  distributed  admittance  between  the  ends  of  the  drift 
a 

tubes  should  be  accounted  for  by  use  of  these  same  equations.  However, 
for  the  cavity  usually  met  in  practice,  these  effects  can  be  approximated 
by  finding  the  admittance  of  the  inner  region  at  the  inner  edge  of  the  gan 
and  adding  this  directly  to  the  radial  line  admittance  in  Equation  (22). 
This  admittance  is  called  the  inner  storage  admittance.  The  second 
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region,  between  the  drift  tubes,  can  be  considered  to  have  a  lumped 

admittance,  thus  treating  it  essentially  like  a  parallel -plate  capacitor. 

This  admittance  is  called  the  annular  admittance. 

The  inner  storage  admittance  has  previously  been  calculated 
4 . 

by  C.  C.  Wang  in  a  manner  consistent  with  the  definition  of  admittance 

used  here.  He  also  makes  the  assumption  of  square  gap  edges  and 

infinite  wall  conductivity,  and  makes  the  same  approximation  to  the 

field  between  the  drift  tube  edges  that  is  used  here  (see  Figure  3). 

However,  Wang  does  put  the  further  restriction  on  his  solution  that 

Pr,  <  2.  405  in  order  to  eliminate  modes  that  can  propagate  down  the 

drift  tubes  in  the  z  direction.  This  is  not  a  serious  limitation  since 

no  klystron  would  be  operated  with  its  fundamental  frequency  in  this 

region.  Wang  also  points  out  that,  to  be  perfectly  safe,  one  ought  to 

have  pr  <  1.  84  so  that  the  circular  wave  guide  TE, ,  mode  will  not 
i  11 

propagate,  since  a  slight  dissymmetry  may  excite  this.  This  again 
gives  no  real  difficulty.  Wang's  expression,  normalized  by  the  use 
of  Equation  (16)  and  written  in  terms  of  the  quantities  used  here,  is 
as  follows: 


where  6  =  R  -  (3r.)  ,  and  R  is  the  nth  root  of  J  . 

Generally  then,  for  consideration  of  modes  having  symmetry  in 

the  ^  direction,  the  expression  in  Equation  (26)  will  hold  for  pr  <  2.  405  . 

1 

To  calculate  resonances  above  this  point  some  other  means  of  approximating 
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the  fields  inside  the  drift  tubes  is  needed,  or  else  one  must  be  content 
with  the  approximation  the  rest  of  the  admittance  function  gives  by  itself. 
In  most  cases,  resonances  up  to  the  fourth  harmonic  of  the  signal 
frequency  can  be  included  and  still  be  within  the  range  stated  above; 
resonances  beyond  this  point  are  not  considered  here. 

As  a  first  approximation  for  the  annular  admittance,  the  ex¬ 
pression  for  the  parallel  plate  capacitor  can  be  used.  This  expression 
yields  the  following: 
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i  2  2, 

ir  (r  -  r.  ) 
_ a  X 


(27) 


At  this  point,  a  better  approximation  for  this  particular  admittance  is 
possible  if  one  admits  that  the  gap  edges  of  the  drift  tubes  are  not  square. 
In  the  usual  cavity  the  gap  edge  is  rounded  so  that,  in  a  cross-sectional 
view,  the  edge  will  approximate  a  semicircle.  The  varying  distance 
between  the  gap  edges  can  be  accounted  for  by  a  straightforward  averaging 
process.  Let  the  corrected  gap  width  be  d* .  From  Figure  4, 

d*  *  *  d  +  t(2-J)  .  (28) 

where  2t  =  r  -  r.  .  Now  simply  replace  d  in  Equation  (27)  by  d*  as 
&  1 

given  above. 

The  annular  admittance  just  calculated  applies  strictly  only 
when  the  cavity  is  resonating  in  the  TEM  mode  and  does  not  apply 
directly  for  the  higher -order  resonances.  However,  this  expression 
is  used  here  when  calculating  the  higher -order  resonances  since  it 
does  provide  a  reasonable  approximation. 
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Figure  4.  Average  Gap  Width  for  Rounded  Edges. 


3.  The  Total  Admittance  Function 

All  of  the  cavity  admittances  associated  with  the  resonances  have 
now  been  calculated.  Each  of  the  e}q>ressions  can  be  combined  to  form 
a  total  admittance  function  or  a  "resonance  equation",  so  called  since 
the  fundamental  cavity  resonance  occurs  at  the  first  zero  of  this  function. 
Before  combining  the  admittance  terms,  shorthand  techniques  can  be 
employed  to  simplify  some  of  the  expressions.  In  particular,  the  cavity 
admittance  of  region  1  is  especially  cumbersome  to  write,  to  say 
nothing  of  calculating  this  in  the  form  given.  Examining  the  expression 
for  the  TEM  mode  admittance  in  Equation  (22),  one  can  substitute  for 
the  Hankel  functions  as  follows: 
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and  obtain 
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If  a  plot  of  this  quantity  is  made  as  a  function  of  P(r  •  r  ),  it  is  very 

b  a 

similar  to  the  circular  cotangent  function  with  P(r^  -  r^)  as  an  argument. 
In  fact  at  large  values  of  this  argument,  the  expression  in  Equation  (30) 
approaches  the  value  of  the  cotangent  as  an  "asymptote.  "  Because  of 
the  similarity  of  these  two  quantities,  we  shall  venture  to  call  the  ex¬ 
pression  in  Equation  (30)  the  Bessel  cotangent  function  of  seroth  order, 
and  define  it  as  follows: 
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In  addition,  for  aid  in  computation  one  may  write  this  as 


b  cotj,P(rj^  -  r^l 


b  cot  Br  (a  -  1 ) 
0  a 


(32) 


where  a  = 


Figure  5  shows  a  comparison  between  the  Bessel  cotangent  function 
and  the  circular  cotangent  function  for  a  =  2 .  The  similarity  mentioned 
earlier  is  very  apparent  in  this  plot.  Figure  6  shows  the  absolute  value 
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Figures.  Comparison  of  b  cot  (06-6)  and  cotG  for 


«»3.0 


26. 


of  b  cot  Pr  (a  >  1)  plotted  as  a  function  of  pr  for  several  values 
0  a  a 

of  a  .  This  also  represents  the  absolute  value  of  the  TEM  mode 
admittance  for  a  short-circuited  radial  lines,  hence  Figure  6  is  so 
labeled. 

As  it  turns  out,  there  are  more  combinations  of  Bessel  functions 
which  are  similar  to  other  trigonometric  functions,  such  as  the  Bessel 
sine  and  cosine  functions.  These  two  sets  of  functions  also  have  many- 
similar  properties;  for  example  the  following  is  true  for  the  Bessel 
cotangent  function: 


b  cot (b  -  a) 


b  cos  (b  -  a) 
b  sin  (b  -  a)  ' 


In  addition,  Bessel  hyperbolic  trigonometric  functions  can  be  derived 
from  the  Bessel  circular  trigonometric  functions  as  follows: 


b  coth  (b  -  a)  =  j  =  j  b  cot  j  (b  -  a) 

b  sin  j  (b  -  a) 


where  b  coth  is  the  Bessel  hyperbolic  cotangent  function.  Substituting 
the  quantities  in  the  equation  above  in  terms  of  modified  Bessel  functions 
gives 


b  coth  (b  -  a) 


K^(a)  Ip(b)  +  Kp(b)  Ij(a) 
Kgta)  l^{h)  -  KQ(b)  1^(3) 


(33) 


This  is  just  the  form  of  the  part  of  Equation  (2Z)  dealing  with  the  cutoff 
modes.  At  large  values  of  argument  this  function  asympotically  approaches 
the  value  of  one  just  as  the  hyperbolic  cotangent  function  does.  A  com¬ 
parison  between  the  b  coth  and  coth  functions  is  shown  in  Figure  7. 

Figure  8  shows  the  b  coth  function  for  several  values  of  a. 
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Figure  7.  Comparison  of  booth  (a  6  -  0)  and  coth  6  for 


29 


Now  by  combining  all  the  admittance  terms  and  making  the 
substitutions  for  the  Bessel  functions  just  discussed,  the  so-called 
"resonance  equation"  is  obtained: 


K 


y  =  -jbcotQ  Cpr^(a-l)H  -  j 


.  8L 

TT  d 


2  \  bcotQ!IPjr^(a-l)Ilcos' 


f  mirh  1  •  2  n 


mird 

“2xr, 


m=l 


m 


.  8L 

IT  d 


b  coth 


p;  r^(a-l)Dcos^(i^)  sin^ 


m=K+l 


-  1 


.  Finally,  the  resonances  may  be  found  by  determining  the  values  of 
frequency  and  the  conditions  on  the  parameters  which  give  the  zeros 
of  this  total  admittance  function.  The  simplest  way  to  find  the  values 
of  the  parameters  that  meet  this  condition  is  to  solve  the  resonance 
equation  graphically. 

The  procedure  adopted  for  this  graphic  solution  is  to  plot  the 

absolute  value  of  each  of  the  admittance  terms  as  a  function  of  6r^  . 

a 

Then  the  intersection  of  the  total  inductive  and  capacitive  curves 
indicates  a  resonance.  This  procedure  is  used  in  lieu  of  a  chart  or 
nomograph  type  of  presentation  because  of  the  large  number  of  para¬ 
meters  involved  and  consequently  the  extensive  calculation  that  would 
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be  necessary.  This  graphic  procedure  does  not  seriously  limit  the 
accuracy  in  finding  resonances,  for  once  one  is  familiar  with  the 
type  of  curves,  all  that  need  be  plotted  is  a  small  section  of  each 
curve  in  the  region  of  resonance.  This  of  course  may  be  on  an 
expanded  scale  so  that  graphic  errors  can  be  on  the  same  order  of 
magnitude  as  any  error  inherent  in  the  resonance  equation. 

As  an  aid  in  using  Equation  (34),  several  parts  of  it  have  been 
calculated  and  are  included  here.  The  TEM  admittance  function  is 
shown  in  Figure  6.  In  addition,  the  Appendix  includes  a  short  table 

for  noninteger  values  of  a  ,  and  integer  values  of  a  are  given  by 

5* 

Marcuvitz.  Figure  6  shows  that  the  short-circuited  radial  line 
behaves  just  like  any  other  short-circuited  wave  guide;  that  is,  if 
one  looks  at  the  short  circuit  in  the  plane  of  the  short  circuit,  the 
admittance  is  infinite.  Then  as  one  proceeds  away  from  the  short 
circuit,  first  the  short  circuit  appears  inductive  and  finally  the 
admittance  goes  to  zero;  then  the  admittance  becomes  capacitive 
and  goes  to  infinity  again.  Then  the  cycle  repeats  itself.  The  same 
behavior  can  be  observed  at  a  fixed  position  in  the  line  by  increasing 
the  frequency  from  zero.  The  important  difference  between  the 
radial  line  and  an  ordinary  guide  (rectangular  or  coaxial)  is  that 
the  admittance  discontinuities  do  not  appear  at  an  integer  number 
of  half-guide  wavelengths  from  the  short  circuit.  The  expression 
guide  wavelength  is  meaningless  in  a  radial  line  due  to  the  ever- 
changing  transverse  area.  Also  in  Figure  6  it  should  be  pointed 
out  that  since  the  absolute  value  of  the  admittance  function  is  plotted, 
the  inductive  parts  have  a  negative  slope  while  the  capacitive  parts 

Marcuvitz  calls  the  "b  cot  "  function  the  "ct"  function. 

0 
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have apositive  slope.  In  the  actual  use  of  these  curves,  it  is  convenient 
to  remember  this. 

The  second  term  of  Equation  (34)  deals  with  the  "propagating" 
higher -order  mode  admittance.  These  modes  behave  essentially  as 
the  TEM  admittance  does,  except  that  the  modes  do  not  appear  until 
the  frequency  is  above  the  cutoff  frequency  for  any  particular  mode. 

s 

These  modes  add  another  set  of  discontinuities  to  the  admittance  plot 
at  us  fa)  for  m  =  1,  2,  3,  . . . .  Note  that  if  the  gap  is  exactly 

m  h  j 

centered,  Tu  ~  "Z  •  only  the  even-numbered  modes  are  excited. 

However  it  has  been  standard  practice  in  such  a  case  to  include  the 
discontinuities  of  the  odd  modes  on  the  admittance  plot,  for  the  physical 
cavity  will  never  have  the  gap  exactly  centered  and  a  resonance  is 
possible  no  matter  how  weak.  This  admittance  term  should  be  calculated 
for  each  particular  cavity;  although  an  attempt  has  been  made  to 
produce  a  universal  type  of  curve  for  this  which  may  be  used  as  an 
initial  approximation.  This  curve  is  shown  in  Figure  9>  This  has 
been  calculated  for  a  centrally  located  gap  in  a  doubly  re-entrant  cavity 
and  is  the  inductive  part  of  the  second  higher -order  mode  for  which 
m  =  2  .  The  type  of  nomenclature  involving  R  and  x  will  be  discussed 
in  the  section  on  optimum  cavity  design  (II-D).  The  value 


for  the  m  =  2  cutoff.  From  Equation  (14), 


2irr. 


2irR  = 


is  merely  the  value  of  pr. 


2ir 


u 


L  'sf^ 
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c  SCALE  FOR  2  RADIAN  GAP  AT 
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Figure  9>  Universal  Curve  Second^Order  Higher  Mode  Inductive 
Susceptance. 
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or 

Zirr 

«  's/ilir  =  pr  =  .  (35) 

Since  the  graphs  presented  here  cannot  be  used  with  any  accuracy,  a 
table  of  this  higher -mode  admittance  is  included  in  the  Appendix  as 
Table  A-2.  The  table  has  been  calculated  on  the  basis  of  a  two-radian 
gap  at  u^/c  =  0.  5  ( u^  is  the  d-c  beam  velocity).  A  method  for 
adjusting  the  values  in  this  table  for  other  gap  lengths  is  also  shown 
in  the  Appendix. 

The  third  term  of  Equation  (34)  involves  the  admittance  due  to 
the  energy  storage  in  the  fields  of  the  cutoff  modes  about  the  gap.  This 
term  has  discontinuities  in  it,  for  as  the  frequency  increases,  each 
successive  mode  will  "pop"  into  "propagation"  and  will  no  longer  con¬ 
tribute  to  this  term.  Little  can  be  done  in  presenting  curves  of  this 
admittance  term  because  of  the  large  numbers  of  parameters  and  because 
it  is  particularly  sensitive  to  the  value  of  cutoff  frequency,  or  to  the 
ratio  which  is  related  to  this  frequency.  Because  of  the  large 

niunber  of  terms  needed  in  this  series  a  computer  solution  is  almost 
a  necessity.  No  known  tables  of  the  Bessel  hyperbolic  cotangent  function 
exist,  so  this  part  is  probably  best  calculated  separately  then  approximated 
by  curve -fitting  techniques  for  use  in  a  digital  computer.  The  data  for 
the  different  values  of  a  used  so  far,  as  well  as  the  range  of  errors 
involved  in  the  curve  fittings, are  presented  in  Table  A-3  in  the  Appendix. 

As  previously  mentioned,  twenty  to  thirty  terms  of  the  series  are 
required  to  keep  the  error  in  the  series  summation  on  the  same  order 
of  magnitude  as  that  in  the  approximations  for  .the  b  coth  fimction,  which 
is  less  than  one  per  cent. 
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The  fourth  term  in  Equation  (34)  ii  the  inner  storage  admittance 


discussed  previously.  This  has  been  calculated  for  all  cases  and  is 
shown  in  Figure  10.  The  admittance  multiplied  by  r  /L  is  plotted 
as  a  function  of  pr  with  as  a  parameter.  For  actual  use  a 

table  with  this  same  information.  Table  A-4,  is  given  in  the  Appendix. 
When  using  this  admittance  on  a  plot  of  the  "resonance  equation"  the 
abscissa  must  be  changed  from  pr.  to  pr^  ,  which  is  a  trivial  con¬ 
version. 

The  final  term  in  Eqviation  (34)  is  the  annular  admittance  term. 

As  is  easily  seen,  this  is  a  straight  line  and  no  difficulties  are  presented 
in  using  it. 

Besides  the  value  of  frequency  at  the  different  resonances; 

Equation  (34)  can  give  the  value  of  Rsh/Q^  at  the  fundamental  resonance. 
At  the  resonant  frequency  <a 


(36) 


Then  the  point  of  intersection  indicating  the  fundamental  resonance 
on  the  plot  of  the  absolute  value  of  the  "resonance  equation"  is  b^  , 


b 

0 


1 


'TirEY 

0 


Combining  Equation  (36)  and  (37)  gives 


(37) 


(38) 


where  is  the  normalizing  admittance  given  in  Equation  (l6). 

4.  Nonsymmetrical  Higher -Order  Modes 

The  assumption  of  symmetry  in  the  ^  direction  of  a  klystron 
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Figure  10.  Inner  Storage  Admittance  at  the  Gap  of  Two  Hollow  Cylinders. 


cavity  is  rarely  adhered  to,  especially  in  the  output  cavity  of  a  high 
power  klystron.  In  present-day  klystrons  some  sort  of  iris  or  antenna 
that  is  large  compared  with  the  cavity  is  usually  used  as  a  means  of 
removing  the  output  power  from  the  tube.  The  heavy,  nonsymmetrical 
loading  caused  by  this  coupling  will  certainly  cause  nonimiformities 
either  in  the  gap  electric  field  or  in  the  cavity  fields  or  both.  A  first 
approximation  of  the  effects  of  any  of  the  resulting  modes  that  have 
variations  in  the  ^  direction  is  made  by  an  analysis  similar  to  the 
one  just  performed.  In  a  distortion  of  the  electric  field  at  the  gap 
caused  by  a  large  coupling  loop,  the  resulting  gap  field  can  be  approximated 
as  shown  in  Figure  11. 
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One  can  now  proceed  to  find  the  fields  given  in  Equation  (1) 
and  the  respective  series  coefficients  in  the  same  manner  as  for 


no  ^  variations.  However  if  one  carried  this  out,  it  would  be 
difficult  to  define  a  current  at  the  gap  since  H  now  depends  on 
^ .  This  difficulty  can  be  overcome  by  defining  an  admittance  Y 
at  the  input  gap  as  follows: 


Y*V^ 

2 


1 

2 


(E  X  H  *)  dA 
*  ^ 


(39) 


where  the  integral  is  over  the  transverse  area  of  the  gap.  Then  to  find 
the  voltage, 


V  s 


_1_ 

2ir 


E  dA 
z 


(40) 


rsr 

a 

which  essentially  is  an  approximation  of  the  voltage  based  on  the  "average" 
of  the  electric  field  at  the  gap. 

The  detailed  steps  are  not  shown  here  since  they  are  involved  and 
lengthy.  If  the  definition  of  the  Bessel  cotangent  function  as  given  in 
Equation  (31)  is  extended  to  include  higher -order  Bessel  functions. 


bcot  (b-a) 
n 


N^(a)  J„(b)  -  N^(b)  J^(a) 


(41) 


A  similar  extension  can  be  made  of  the  Bessel  hyperbolic  functions. 
Using  Equations  (39)  and  (40)  and  these  new  combinations  of  Bessel 
functions,  the  normalized  admittance  for  the  short-circuited  radial 
line  is  as  follows: 
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m=K+l 


In  this  equation,  the  last  three  terms  are  introduced  by  the  modes  with 
^  variations.  Of  course,  these  terms  go  to  zero  if  either  v  or  9  is 
zero,  which  is  as  it  should  be.  since  this  gives  the  symmetric  case 
originally  calculated. 

As  it  turns  out  the  magnitude  of  these  last  terms  is  small 

compared  with  the  admittance  terms  which  do  not  depend  on  variations 

except  at  the  discontinuities.  In  the  fifth  and  sixth  terms  of  Equation  (42), 

the  discontinuities  are  caused  by  the  denominator  and  they  occur  at  the 

same  values  of  6r  as  do  the  discontinuities  of  the  admittance  function 
a 

for  the  symmetric  case.  The  fourth  term  has  discontinuities  which 
depend  on  the  roots  of  terms  like 

*  N„(par  )  J  (Pr  )  =  0  .  (43) 

nan  a  n  ana 

However  for  a  fixed  value  of  a .  this  equation  is  satisfied  at  values  of 
pr  ,  which  are  nearly  independent  of  the  order  n . 

A 

The  next  few  sections  show  that  the  higher -order  resonances 
of  the  cavity  occur  largely  in  the  region  of  some  of  the  discontinuities. 

Thus  the  fact  that  the  nonsymmetrical  modes  only  contribute  a  significant 
admittance  near  these  discontinuities  means  that,  to  the  first  approximation, 
these  modes  and  their  admittance  terms  can  be  neglected  when  compared 
with  the  whole  admittance  function.  At  the  best,  they  can  shift  only 
slightly  the  higher-order  resonances  which  occur  near  the  discontinuities. 

It  should  again  be  emphasized  that  this  treatment  of  the  non¬ 
symmetrical  modes  of  a  klystron  cavity  is  only  a  first  approximation. 

It  is  not  meant  to  brush  aside  these  modes  with  only  a  cursory  glance 
as  it  might  at  first  appear.  Certainly  a  careful  treatment  of  them  would 
prove  worthwhile  since  any  practical  output  cavity  of  a  present-day 
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high-power  klystron  will  have  these  nonsymmetric  modes  present. 

It  is  just  that  the  first  consideration  of  this  whole  problem  has  been 
to  treat  those  first  few  higher -order  resonances  which  are  known 
to  exist  experimentally  and  which  involve  primarily  the  symmetrical 
type  of  mode. 

B.  CONTROL  EXPERIMENTS  AT  THE  FUNDAMENTAL  RESONANCE 
Several  experiments  were  performed  as  a  check  on  the  results 
that  Equation  (34)  predicts.  The  first  of  these  was  the  comparison  of 
the  fundamental  resoance  of  a  gridded-gap,  singly  re-entrant  cavity 
calculated  by  curves*  similar  to  the  ones  in  Moreno^  and  calculated 
by  Equation  (34).  This  equation  may  be  revised  to  apply  to  this  type 

of  cavity  by  letting  r. - ^0.  This  essentially  approximates  the  total 

capacity  inside  the  gap  as  the  d-c  parallel  plate  capacity.  Table  1 
shows  the  results  of  these  calculations.  Agreement  is  very  good  if 

Table  1.  Calculated  Fundamental  Resonance  for  a  Sinj^y  Re-entrant  Cavity. 


Cavity  Parameters 

Calculated  Resonant  Frequencies 

1 

B 

^a 

(cm) 

Hansen's 

Curves 

Eq.  34 

Percent  of 
Difference 

1.  0 

0.  61 

2.0 

1.59 

3000  Mc/s 

3060  Mc/s 

1.9 

1.0 

0.  122 

m 

2.50 

955  Mc/s 

959  Mc/s 

0.  4 

one  considers  that  for  the  first  example^  the  d-c  approximation  to  the  gap 
capacity  is  applied  to  a  very  large  gap.  Hansen's  curves  have  been  used 


An  enlarged  set  of  the  curves  appearing  in  Moreno  were  used. 

Although  not  so  named  in  this  reference,  these  curves  are  called  Hansen's. 
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here  for  a  comparison  since  these  are  known  to  be  an  accurate  means 
of  computing  resonances  for  this  type  of  cavity.  In  addition  Figures 
12  and  13  show  the  plot  of  each  term  of  Equation  (34)  for  these  cavities. 
About  seventeen  terms  of  the  infinite  series  for  the  admittance  of  the 
"nonpropagating"  modes  were  necessary  for  these  cavities.  In 
Figures  12  and  13  the  quantity  value  of  pr^  at  resonance. 

The  second  experiment  performed  was  the  measurement  of  the 
fundamental  resonant  frequency  and  Rsh/Q^  of  an  existing  cavity. 

This  is  a  doubly  re-entrant  cavity  with  solid  posts,  which  represent 
a  gridded  gap.  One  of  the  posts  was  threaded  so  that  it  could  be  moved 
to  adjust  the  width  of  the  gap.  Since  only  one  of  these  posts  moved,  this 
meant  that  the  gap  was  off  center  except  at  one  position,  which  occurred 
when  the  gap  was  as  wide  as  possible  (corresponded  to  the  highest 
resonant  frequency).  Equation  (34)  was  used  to  predict  resonances 
for  several  gap  widths  and  the  value  of  Rsh/Q^  at  one  of  these. 
Hansen's  curves  were  used  to  predict  two  resonances,  one  with  the 
gap  centered  and  the  other  with  the  gap  slightly  off-center.  Hansen's 
curves  apply  strictly  to  a  singly  re-entrant  gridded -gap  cavity,  but  a 
symmetrical  doubly  re-entrant  cavity  may  be  treated  by  considering 
one -half  of  this  cavity.  The  nonsymmetrical  cavity  was  divided  in 
the  middle  of  the  gap  and  the  predicted  resonance  for  each  half  was 
averaged. 

Two  other  means  of  calculating  Rsh/Q^  were  used  to  see 
how  well  these  compared  with  the  experimental  value.  The  first  of 
these  was  from  the  same  data  used  for  Hansen's  curves,  ^  while  the 
second  was  from  an  expression  often  used  for  this  type  of  cavity; 
that  is; 
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Figure  12.  Sample  Problem  -  Singly  Re-entrant  Gridded -Gap  Cavity. 
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SAMPLE  PROBLEM 


Doubly  Re-entrant  Gridded  Gap  Cavity 
Gap  Distance  is  Mechanically  Variable 


Table  2.  Resonances  of  S-Band  Cavity 


Measured 

Hansen's  Curves 

Predicted  Values 

_ _ 

(cm) 

a  0 

(cm) 

%  Diff.  to 
Pleasured 

»'a)o 

\ 

(cm) 

%  Diff.  to 
Measured 

0.  3125 

0.  3182 

9.  43 

1.  77 

0.  3170 

9.46 

1.46 

0. 2971 

'  ■ 

— 

■  — 

0. 3000 

10.  00 

1. 00 

0. 2857 

♦O. 2906 

10.  32 

1.  71 

0.2895 

10.  36 

1.  33 

loil 

0.  2703 

— 

'  ■ 

0.  2730 

10.  99 

0.99 

Approximate  Due  to  Asymmetrical  Cavity 


Table  3.  R  . /Q-  at  X.  =  10.  50  cm 
sh'  0  0 


Method  Used 

^sh'^o  * 

%  Djff.  to 
Measured 

Measured 

65.  3 

Hansen 

56.  8 

13.  0 

Simple  Formula 

54.  1 

17.  2 

Predicted 

60.  5 

7.  4 

(44) 


The  experimental  and  calculated  data  for  this  cavity  are  listed  » 

in  Tables  2  and  3.  The  values  marked  "predicted"  are  from  Equation  (34). 

The  cavity  dimensions  are  shown  in  Figure  14. 

These  simple  control  experiments  indicate  that  good  accuracy 
may  be  obtained  from  calculations  using  Equation  (34).  The  experiment 
including  Rsh/Q^  in  particular  indicates  that  this  equation  gives  better 
accuracy  for  this  quantity  than  the  other  methods  used  here. 

C.  SPECIFIC  APPLICATION  TO  HIGHER -ORDER  RESONANCES 

The  utility  of  the  admittance  function  as  applied  to  fundamental 
resonances  has  now  been  demonstrated  by  experiment.  It  remains  to 
be  shown  that  this  function  predicts  the  higher -order  resonant  frequencies 
with  a  reasonable  accuracy.  For  this  purpose  a  model  of  an  L-band 
klystron  cavity  was  constructed.  This  cavity  was  doubly  re-entrant  and 
had  a  gridless  gap  with  rounded  gap  edges.  The  drift  tubes  of  this  model 
were  extended  far  enough  in  the  z  direction  so  that  the  discontinuity  at 
the  ends  of  the  tubes  would  not  be  reflected  back  into  the  cavity.  Figure  15 
shows  the  approximate  sizes  involved  in  this  model.  The  gap  was  very 
nearly  centrally  located  in  the  z  direction. 

Two  rectangular  loops  were  made  so  that  the  cavity  could  be 
tested  as  a  transmission -type  cavity.  The  loops  were  approximately 
3/4"  X  1  1/4",  although  for  the  resultant  tests  these  were  only  projected 
1/4  in.  into  the  cavity  since  this  position  gave  the  best  match.  When 
the  cavity  was  tested  for  resonances  through  the  third  harmonic  of  the 
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Figure  15.  Model  of  L>Band  Klystron  Cavity. 


fundamental,  six  resonances  other  than  the  ones  predicted  were 
observed.  It  turned  out  that  these  were  resonances  of  ^-directed 
(nonsymmetrical)  modes  excited  by  the  presence  of  the  two  coupling 
loops  in  the  cavity.  Because  of  the  method  of  exciting  these  modes, 
a  simple  approximation  to  the  expected  resonant  frequencies  can  be 
made. 

For  this  consider  Equation  (1)  with  m  =  0.  An  admittance 
for  region  1  of  the  cavity  may  be  defined  as 
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(45) 


y  '  "TT 

0  z 


The  series  coefficients  of  the  field  expansions  can  be  found  using  the 
boundary  conditions  at  the  outer  radius  of  the  cavity,  since  only  the 
relationship  between  and  in  Equation  (45)  is  needed.  Making 
this  substitution  and  using  the  definitions  of  the  Bessel  trigonometric 
functions,  the  admittance  for  the  nth  mode  follows: 


n 

'n  =  j  I 


+  bcot  Ppr  (a-l)J 


(46) 


This  eiepression  may  be  combined  with  the  inner  storage  admittance  given 

in  Equation  (26)  and  the  annular  admittance  given  in  Equation  (27)  to  give 

the  resultant  admittance  function  for  the  H  modes.  Of  course  this 

Onl 

is  only  approximate  since  the  other  susceptances  will  not  apply  precisely 
to  these  ^-dii^ccted  modes. 

The  computation  of  these  resonances  was  performed  and  showed 

that  the  first  eight  4-‘di*'ected  modes  had  a  resonance  within  the  frequency 

test  range  of  the  cavity.  A  comparison  of  the  model  lobes  in  ^  and  the 

positions  of  the  two  coupling  loops  showed  that  extremely  v  eak  coupling 

was  present  for  just  two  modes,  precisely  the  two  that  had  not  been 

measured.  Agreement  between  the  measured  and  calculated  values  of 

resonance  for  the  other  six  ^-directed  modes  was  reasonable.  Differences 

ranged  from  less  than  3  per  cent  to  about  19  per  cent  (this  large 

difference  was  for  the  H  mode,  which  occurred  close  to  the  cutoff 

041 

frequency  of  the  H  mode). 

'  201 

The  resonances  of  the  symmetrical  modes  of  this  cavity  were 
calculated  by  Equation  (34).  For  this  case  approximately  25  terms  of 
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the  infinite  series  representing  the  cutoff  modes  were  required.  To 
give  an  idea  of  what  these  terms  look  like  graphically,  a  complete 
plot  for  this  cavity  is  given  in  Figure  l6.  The  experimental  results 
and  the  calculated  predictions  for  the  symmetrical  modes  are  shown 
in  Table  4. 


Table  4.  Resonances  of  lj>Band  Cavity. 


Radial 

Line 

Mode 

Coaxial 

Mode 

Calculated 
f  Resonance 
(Mc/s) 

Measured 

Resonance 

(Mc/s) 

Error  in 
Calculation 

TEM 

1 

TM 

010 

819.  7 

823.  0 

-0.  40% 

101 

TM^ 

Oil 

1004.  2 

1048.  0 

-4.  18% 

H, 

201 

012 

2301.  8 

2538.  9 

-9.  34% 

As  can  be  seen  from  the  results,  the  admittance  function  is  excellent 
for  the  fundamental  resonance  and  is  a  good  approximation  at  the  first 
few  higher 'Order  resonances.  The  error  is  expected  to  increase  as 
one  considers  these  higher -order  resonances  since  the  admittance 
terms  connected  with  the  "gap  capacity"  are  not  precise,  as  explained 
in  tlje  section  on  theoretical  analysis. 

D.  OPTIMUM  CAVITY  DESIGN 

All  of  the  data  and  equations  given  up  to  this  point  are  quite 
convenient  for  analyzing  an  existing  cavity.  However,  to  design  a 
cavity  by  using  this  theor)^  the  problem  of  where  to  begin  presents 
itself.  For  this  reason  it  is  desirable  to  present  a  design  technique 
and  the  general  limits  for  certain  of  the  cavity  parameters.  In  the 
interest  of  serving  the  most  usual  cavity,  the  following  type  of  doubly 
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Figure  16.  Total  Admittance  Function  of  the  L>-Band  Cavity  Model  (r  =  1.  750  in.  ,  a  =  1. 85). 


re-entrant  gridless -gap  cavity  is  considered: 


(1) 

(2) 

(3) 


(i  •  ?)  • 


the  gap  is  centrally  located 
the  gap  is  one  or  two  radians  long, 

the  d-c  electron  beam  velocity  is  one -half  the  speed  of  light, 


The  first  of  these  restrictions  means  that  the  modes  having  an 
odd  number  of  variations  in  the  z  direction  are  only  weakly  excited. 
The  second  and  third  restriction  give  a  measure  of  the  gap  length  or 
strength  of  excitation  of  the  higher -order  modes.  The  cavity  is  to  be 
designed  so  that  harmonics  of  the  signal  frequency  up  to  and  including 
the  fourth  wil'  not  occur  at  or  near  any  higher -order  resonances  of  the 
cavity. 


If  one  examines  the  curves  of  the  TEM  inductive  susceptance, 

one  can  see  that  on  the  basis  of  the  discontinuity  points  of  these  curves 

and  the  type  of  cavity  cross-sections  represented,  the  usual  range 

of  a  where  o  =  —  would  be  from  1.  5  to  3.  0.  A  lower  value  of  a 
*’a 

would  give  an  extremely  "thin"  cavity  which  would  be  unsuitable  for 

high  -power  applications.  A  larger  value  of  a  would  necessitate  a 

very  low  value  of  (pr  )  in  order  to  avoid  trouble  with  the  harmonics 

a  0 

at  the  discontinuities  of  the  inductive  curves.  Before  qualifying  the 
selection  of  a  particular  value  of  a ,  one  should  know  something  about 
the  effects  of  the  other  parameters.  The  length  of  the  cavity 
controls  the  higher-order  cutoff  frequencies,  as  shown  in  Equation  (14). 
From  this  equation,  one  can  find  the  value  of  the  first  higher  mode 
cutoff  by 

irr 

pr^  =-T~  =  ifR  •  (47) 
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This  can  be  rewritten  to  apply  to  the  value  of  pr  at  the  fundamental 

a 

resonance  (pr  )  by  introducing  the  parameter  x  as  follows: 

^  0 


This  equation  now  makes  it  possible  to  find  a  value  for 
length  to  cavity  length  ratio.  For  a  one -radian  gap, 


(48) 

or  the  gap 


p  d  =  ~  d  =  1 
e  u 

0 

Substituting  from  Equation  (48)  gives 


which  is  for  a  one -radian  gap.  For  a  two-radian  gap  it  is  obvious  that 


Now  returning  to  the  T£M  admittance  curves,  one  can  see  that 
the  likely  range  of  ^or  1. 5  <  a  <  3.  0  is  from  0.  3  to  1.0  (perhaps 

even  a  bit  higher  for  a  =  1.5).  pacifically,  one  can  now  recommend  the 
range  of  ^or  a  particular  a .  Figure  1 7  shows  these  recommended 

values  for  a  =  1.  5,  2.  0  and  3.  0. 

From  Equation  (48),  since  the  first  higher -mode  cutoff  frequency 
is  normally  between  the  fundamental  and  second  harmonic  of  the  signal 
frequency,  a  condition  on  the  parameter  x  is  obtained  as  follows: 

1.0  <  X  <  2.0  (51) 

One  is  now  in  a  position  to  make  calculations  for  some  specific 
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cases,  then  interpolate  between  these  so  as  to  cover  all  the  values  of 
the  parameters  mentioned  here.  Cavities  with  the  following  properties 
are  considered: 

a  :  1.5,  2. 0,  3. 0 

X  :  1.25,  1.  50,  1.  75 

(pr  )  :  0.2,  0.6,  0.9 

»0 

P  d  :  1,  2 
e 

Figure  18  shows  the  cross-sectional  views  of  the  cavities  these  values 
represent  so  that  the  physical  effect  of  varying  one  of  these  parameters 
can  be  observed. 

The  effects  of  these  parameters  on  the  higher -order  resonances 

of  a  cavity  can  be  observed  by  marking  the  values  of  on  a  scale 

of  pr  ,  and  indicating  the  harmonics  of  these  values  of  (Pr  )  .  Then 
a  ^0 

the  cutoff  values  of  pr^  for  the  particular  values  of  x  chosen  can  be 
marked.  Once  this  is  done,  the  inductive  susceptance  of  the  second 
higher -order  mode  m  =  2  in  z  is  calculated  and  plotted  as  a  function 
of  pr^  .  This  susceptance  is  added  to  the  inductive  susceptance 
produced  by  the  TEM  mode  which  gives  the  total  inductive  susceptance 
of  the  cavity.  From  past  experience  one  knows  that  the  first  few 
higher -order  resonances  occur  in  the  region  of  the  large  values  of 
the  higher-order  mode  inductive  susceptances  and  can  stipulate  the  "safe" 
range  of  the  parameter  x.  The  conceptual  process  is  not  obstruse 
but  the  actual  solution  is  involved  since  one  must  also  be  careful  to 
keep  the  cavity  resonances  away  from  the  harmonics  at  values  of 
«(Prj^)^  other  than  the  three  chosen  here  for  the  sample  calculations. 

The  results  of  these  calculations  are  shown  in  Figures  19,  20,  and  21. 


-54- 


(FOR  a  =  3.0.  2.0,  1.5  AND  WITH  ONE  AND  TWO  RADIAN  GAPS  AT  =  0.5) 


O 

<g  «  fe  “ 
60  6^ 
“  ^  “  01 
(r^«o  p 
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RADIAN  GAP  AT  l^/C  •  0.5  I  RADIAN  GAP  AT  L^/c*0.5 


56 


Figure  19.  Recommended  Ranges  of  Parameter  x  as  Indicated  by  Shaded  Areas 


a » 2.0  _ ,  a  *  2.0 
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Each  of  these  figures  is  for  a  specific  value  of  a  (1*  2.  0,  and  3,  0) 

and  shows  the  recommended  range  of  x  for  any  (Pr  )  .  Each  of 

a  0 

these  figures  has  been  subdivided  into  curves  for  a  one-  and  a  two- 
radian  gap  so  that  a  rough  interpolation  may  be  made  for  any  other 
expected  gap  length. 

The  cavity  has  not  been  completely  specified  yet,  for  there 

is  the  ratio  to  be  specified.  This  can  actually  be  chosen  last, 

since  it  affords  a  minor  adjustment  on  the  resonant  frequencies  of  the 

cavity.  Now  the  ratios  of  all  the  physical  measurements  of  the  cavity 

are  completely  specified,  although  a  “scaling  factor"  must  still  be 

determined.  This  may  actually  be  found  from  the  value  of  (^r  )  . 

a  0 

The  operating  frequency,  or  value  of  p ,  will  be  known  when  the  tube 
design  is  begun  and  from  the  resulting  electron  beam  characteristics 
a  value  of  r^  may  be  estimated.  Once  this  particular  physical  para¬ 
meter  is  known,  the  rest  of  the  cavity  parameters  can  be  found  from 
the  information  given  above,  and  in  the  figures  following. 

Perhaps  it  should  be  pointed  out  that  this  design  does  not  consider 
the  effect  of  only  a  few  higher -order  modes;  rather  it  accounts  for  all 
the  symmetric  higher -order  modes  which  are  of  the  type. 

Table  5  lists  the  important  ones  that  were  considered  in  specifying 
the  design  technique  outlined  here.  Each  mode  is  identified  in  both 
the  radial  and  coaxial  nomenclature.  In  addition.  Figure  22  shows 
the  electric  field  configurations  on  cavity  cross  sections  of  several 
typical  modes.  This  field  structure  for  the  other  modes  may  be  easily 
obtained  by  extension. 
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Table  5.  Modes  Considered  in  Design 


Mode 

Description 

Radial  Line 
Designation 

Coaxial  Line 
Designation 

Fundamental 

TEM 

1 

TM 

010 

1*^  Cutoff  Mode 

”l01 

^“0.. 

2"**  Cutoff  Mode 

^201 

™0.3 

3*^**  Cutoff  Mode 

“30. 

^“0.3 

2^*^  Order  Fundamental 

TEM 

TM--^ 

020 

2“**  Order  of  1*^  Cutoff 

“.03 

TM.,, 

021 

2"^*  Order  of  2“**  Cutoff 

”202 

TM^„ 

022 
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III.  PERTURBATION  THEORY  IN  RESONANT  CAVITIES 


The  general  case  of  perturbation  of  cavities  by  the  insertion  of 
some  object  has  not  been  treated.  Slater^  treated  the  general  case  of 
boundary  perturbation  but  did  not  consider  an  object  apart  from  the 
boundary.  Over  the  past  few  years  several  articles^’  have 

treated  cavity  perturbation,  but  largely  from  the  experimental  point 
of  view.  A  few  articles  calculate  the  calibration  constant  mentioned 
later  in  this  section  for  several  types  of  beads,  but  most  use  an  assumed 
form  of  perturbation  equation.  Casimer  does  give  a  brief  derivation 
for  certain  types  of  perturbation  in  terms  of  total  magnetic  and  electric 
moments  which  will  depend  on  field  strength  and  the  sise  of  the 
perturbing  object.  Here  perturbation  theory  is  developed  in  general, 
with  the  approximations  and  assumptions  clearly  shown.  Rather  than 
deal  with  moments,  the  equations  are  developed  so  that  they  are  true 
in  general,  with  only  a  calibration  constant  to  differentiate  between 
types  of  beads.  This  calibration  constant  is  independent  of  field 
strength  and  siae  of  the  perturbing  bead.  It  depends  only  on  the 
geometry  and  orientations  of  the  bead. 

A  convenient  way  to  handle  perturbation  is  to  use  normalised 
orthogonal  functions  for  the  fields  in  a  cavity.  A  short  development 
of  these  is  shown  in  the  next  sections.  For  those  who  are  unfamiliar 
with  this  type  of  treatment,  an  illustrative  problem  which  shows  how 
these  functions  might  be  used  has  been  included. 

A.  ORTHOGONAL  ELECTROMAGNETIC  FIELD  FUNCTIONS 

In  dealing  with  cavities,  it  is  often  convenient  to  expand  the 
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fields  in  terms  of  a  series  of  orthogonal  functions.  This  reduces  the 
solutions  of  Maxwell's  equations  in  the  cavity  to  the  problem  of  finding 
the  coefficients  of  expansion.  The  form  of  this  expansion  is  immediately 

g 

suggested  by  Helmholtz's  theorem,  which  states  that  any  vector  field 
F,  that  is  finite,  uniform,  and  continuous,  maybe  expressed  as  the  sum 
of  a  gradient  of  a  scalar  4>  ^  curl  of  a  zero -divergence  vector  A,  as 

F  =  grad  ^  +  curl  A  (51) 

The  vectors  grad  ^  and  curl  A  are  often  termed  "ir rotational"  and 
"solenoidal"  respectively,  since  one  has  zero  curl  while  the  other  has 
zero  divergence. 

6 

This  expansion  for  the  cavity  fields  was  performed  by  Slater 

9 

several  years  ago.  However,  Kurokawa  recently  pointed  out  a 
fundamental  error  in  Slater's  derivation.  Since  Slater's  expansion 
has  not  been  entirely  corrected,  it  is  worthwhile  to  indicate  this 
expansion  here.  For  purposes  of  reference  the  nomenclature  used 
here  will  in  large  measure  agree  with  that  used  by  Slater. 

In  general  a  cavity  is  considered  of  volume  V  bounded  by 
short-circuited  surfaces  s  and  open-circuited  surfaces  s^  which 
are  so  named  because  of  the  type  of  boundary  conditions  at  the  surfaces. 
For  the  expansions  of  the  scalar  fields  in  the  volume  there  are  two 
sets  of  orthogonal  functions  which  satisfy  the  following  relationships 
and  boundary  conditions: 

2  2 
V  4  +  k  =  0 

8^1 

=  0  (on  s)  ,  -5r^=  O(ons')  (52) 

a  on 
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0 


2  X  X  1,  2  . 

V 


^  =0  (on  s') 

X 


8 


=  0  (on  s) 


(53) 


Similarly  for  the  expansions  of  the  vector  fields  in  the  volume,  two 
sets  of  orthogonal  functions  may  be  used: 

2  2 
7  tp  +  k  «li  =0 

"IP  P  ^ 


n  X  =  0 
“  “P 


V  •  4,  =  0 

“  -P 


/  (on  s) 


n  X  V  X  <(i  =0 

“  “  ”P 

n  •  4*  =0 

“  “P 


>  (on  s') 


(54) 


2  2 

V  #  +  k  f  =0 

-q  q-<J 


n  X  i  s  0 

”  q 


7  •  i  *  0 
”  -q 


>  (on  s') 


n  X  V  X  f  s  0 

-  -  q 


n  •  £  =  0 

-  ■*-q 


(on  s) 


(55) 


Some  of  the  functions  of  Equations  (54)  and  (55)  have  the  common 
eigenvalues  k  (a  =  1,  2,  3,  . . . )  and  are  related  by  the  expressions 

A 


k  tV  =7x4 
a  -p  —  “q 


k  £  s  V  x  iL 

a  q  “  ^ 


(56) 


These  particular  «)i  and  £  will  be  denoted  as  E  and  H  respectively, 

— p  q  —a 

so  that  Equation  (56)  becomes 


k  E,  =  V  X  H 
a  —a  —  —a 


k  =  V  X  E, 

a  —a  —  —a 


(57) 


E  and  H  will  now  satisfy  the  boundary  conditions: 
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n  X  E 
—  —a 


n  • 


H 

—a 


0 

>  {on  s) 

0 

i 


n  *  E  = 
—a 

n  X  H  = 
a 


(on  s')  (58) 


All  the  other  t)'  which  do  not  satisfy  Equation  (56) 

are  denoted  as  F  (a  =  0,  1,  2,  . .  .  )  and  G  (X  =  0,  1,  2,  . .  ,  ).  Since 
“o  “X 

these  functions  have  zero  curl,  they  satisfy  the  relations: 


la  =1  ^ 

^X  ^  (X  /  0) 


(59) 


where  the  4*  6  are  from  Equations  (52)  and  (53).  The  particular 

solutions  F  (o  =  0)  and  G.  (X  =  0)  correspond  to  the  electrostatic 
—a  —X 

and  magnetostatic  fields  under  certain  conditions;  however  these  are 
neglected  here  since  they  are  not  important  in  the  development  of 
perturbation  theory. 

By  examining  Equation  (57)  and  (59),  it  is  apparent  that  grad  ^ 
in  Equation  (51)  corresponds  to  the  F  and  G.  functions  (the 
"ir rotational"  functions)  while  curl  A  corresponds  to  the  E  and  H 
functions  (the  "solenoidal"  functions).  Thus  one  can  now  expand  the 
the  electromagnetic  fields  by  using  Equation  (51)  as  a  guide.  The 
general  form  can  be  seen  by  considering  the  electric  field  which  can 
be  expressed  in  terms  of  F^  and  E^  (since  these  have  boundary 
conditions  similar  to  the  actual  field  E  in  a  cavity): 


Here 


E  =  yeE  +  yfF 

^  ^  &  ^3-  Cl  *-*01 

3  Cl 

e  and  f  are  the  series  coefficients  of  the  expansion, 
a  a 


(60) 

These 
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coefficients  can  in  general  be  found  by  using  the  orthogonal  properties 
of  the  functions  ,  and  and  by  considering  them 

normalized  over  the  cavity  volume  so  that 


6  =  0  ,  n  /  m 

nm 

8  =  1  ,  n  =  m 

nm 


(61) 


Of  course,  the  "cross  products"  of  different  functions  integrated  over 

the  volume  are  always  zero,  that  is,  equations  of  the  form 

\ 


(62) 


Thus  to  find  the  series  coefficient  e_  in  Equation  (60),  multiply  through 

by  E  and  integrate  over  the  volume  to  obtain 
•“b 


By  virtue  of  Equations  (61)  and  (62)  the  only  term  remaining  on  the  right 
side  occurs  when  b  =  a  which  gives 


I 


(E  .  E^)  dv  =  e^ 


(63) 


The  coefficients  f^  may  be  found  in  a  similar  way  and  are  as  follows: 


/, 


(E  •  F  )  dv  =  f 
y  —  o  a 


(64) 


5 

nm 
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It  is  easily  seen  that  all  further  series  coefficients  of  expansion  will 

have  the  form  of  Equations  (63)  and  (64). 

Now  in  considering  the  expansion  of  the  magnetic  field,  since 

the  divergence  of  this  field  is  always  zero,  it  is  a  temptation  to  expand 

this  in  terms  of  a  solenoidal  function  only,  for  this  also  has  zero 

divergence.  One  must  remember,  however,  that  the  magnetic  field 

being  considered  does  not  exist  in  infinite  space,  but  in  a  bounded 

space,  that  is,  ins  de  a  cavity.  Because  of  this  an  irrotational  function 

is  also  required  to  expand  the  magnetic  field.  Mathematically  this 

may  be  seen  by  examining  the  general  solution  for  ^  and  A  in 

o 

Equation  (51)  which  's  as  follows: 


where  R  =  r  -  r'  and  2^'  is  a  differential  operator  on  r'.  Even 
though  V  *  F  =  0  there  is  no  assurance  that  F  =  v  x  A  only.  Hence, 
the  magnetic  field  must  be  expanded  in  terms  of  the  (solenoidal) 
and  G  (irrotational)  fimctions. 

Also  because  of  similar  boundary  conditions,  J  (electron  current) 
is  expanded  in  terms  of  E  and  F  ,  and  p  (electron  charge  density)  is 
expanded  in  terms  of  t|i  .  This  gives  the  following  field  expansions: 
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H 


J 


(H  •  H  )dv  + 

a 


I 


(66) 


In  order  to  use  Maxwell's  equation^  expansions  are  needed  for  the 
curl  and  divergence  of  the  electric  and  magnetic  fields.  Since  2  ^  — 
behaves  like  H  this  can  be  expanded  in  terms  of  and  as  follows: 


This  can  be  put  into  a  more  tractable  form  by  making  use  of  the  vector 
identity 

=  *  (68) 

so  that 

Y  •  CE  X  (2  X  E^)l  =  2xE  -  Ea  '  E  ’  Z  ^  Z  ^ 

or 

7-(ExkH)  =  7xE-kH,  -  E-  k^E 
Integrating  this  over  the  cavity  volume  and  using  Gauss'  theorem  gives 
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where  n  is  the  outer  normal  to  the  surface.  The  order  of  the  vectors 


can  be  permuted  in  the  surface  integral  and  with  n  x  H  =0  on  s', 


f  (7  xE  ■  H 


,, -  =  k  J  {E-  E  )dv  , 

y  —  —  -a  ay  —  -a  ^-a—  - 


+  H  •  (n  X  E) 


dA  .  (69) 


An  expression  for  the  second  term  of  Equation  (67)  can  be  obtained  by 
using 

7  •  (Ex  Gj^)  =  •  2  ^  ^  ~  f  * 

=  •  2  X  E 

since  the  curl  of  an  irrotational  function  is  zero.  Again,  integrating 
over  the  cavity  volume,  and  using  Gauss'  theorem  gives 


/ 


(n  •  E  X  G. )  dA  =  J  {7  x  E  >  G.)  dv 

s  s'  -  -  -\  ^  j 

One  can  permute  the  order  of  the  vectors  in  the  surface  integral  and 
with  n  X  G^  s  0  on  s',  obtain 


/ 


(^  X  E  •  dv 


=  /  •  (n  X  E) 


dA 


(70) 


V  s 

Equations  (69)  and  (70)  may  be  inserted  into  Equation  (67)  to  give  for 
the  expansion  of  2  ^  ^ 

(£•  £Jdv  +  u  (n  X  £)  dA 


+  Z  G.  /  G^-  (nxE 


£)  dA 

X  8 - -  " 

In  a  corresponding  manner  ^x  H  behaves  like  E  so  that 


(71) 


-69- 


S‘Zs, 


[  f 

i  (H  •  H  )  dv  +  J  E  •  (n  X  H)  dA 

mm  “"A.  “"a 


-a  -a 


+  Z  la 


dA  (72) 

^  aw 

a  “s’ 

Since  is  related  to  H,  expand  7  ■  B  in  terms  of  the 
as  follows: 


f  (V  •  B 


7  •  B  =  Z  4» 

"  ”  \ 

However,  by  using  the  vector  relation 


^  B^Jdv  . 

\  a  V  “  “ 


(73) 


7  •  (ij^B)  =  B  •  7  +  ^,^7  •  B 

=  ®  *  -X  *X  5  *  - 

integrating  over  the  cavity  volume,  and  using  Gauss'  theorem,  one  obtains 

f  7  •  B)dv  =  -k  j' (B  •  G  )dv  +  (A  B  •  n)dA  . 

V  X  -  “  V  y  -  -X  8,  s'  “ 

Since  =  0  on  S'  the  surface  integral  over  this  part  of  the  area  of 
zero.  Now  Equation  (73)  becomes 


.  .-I.,  k  /..-a..-.  /».  B  *  n)  dA 


.  (74) 


In  a  corresponding  manner  one  can  obtain' for  ^  B 
y  .  D  =  Y.K 


-k 


i  (D  .  F^)dv  +  / 


(4>jj^  D  •  n^)  dA 


(75) 
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Equations  (66),  (71),  (72),  (74)  and  (75)  give  the  expansions  for 
the  various  electromagnetic  field  quantities  which  can  now  be  used  in 
Maxwell's  equations. 

B.  ILLUSTRATIVE  EXPANSION  OF  CAVITY  FIELDS 

If  one  lacks  experience  in  handling  cavity  fields  with  the  type 
of  expansions  made  here,  it  might  prove  helpful  to  make  a  specific 
expansion  at  this  time.  Consider  a  rectangular  cavity  with  ideal 
metal  walls  as  shown  in  Figure  23,  which  is  operating  in  the  TE^^^mode. 


AY 


Figure  23.  Rectangular  Cavity. 
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The  fields  in  this  cavity  may  be  expressed  as  follows: 


H  =  H-  cos  6  X  cosS  y  sinp_z  e 
.  z  0  X  y  * 


jut 


P_Px 

H  =  -H-  -'■■■■w-  sin  a  X  cosp  y  cos  P  z  e 
X  0  a  2  X  y  z 


H  =  -H. 


fcos  p  X  sinp  y  cos  p  z  e 

V  V  Z 


P 

=  j  .  ■■  y  *  COS  p  X  sinP  y  sinp  z  e 


jut 


jut 


(76) 


where 


P 

E  =  -j-- - i -  sinp  X  COS  py  sinp  z  e 

y  Q  i  *  y  * 

^c 


-  nit  ^  mir  pit 

P*  =  a“  •  Pv  =  "IT  i  P,  =  “ar 


Pc"  =  Px'  " 


(77) 


The  particular  cavity  chosen  needs  only  the  solenoidal  parts  of  the 
expansion  for  £  and  H  to  describe  these  fields.  One  can  find 

/(£  •  EJdv  by  taking  /  (E  *  E’^)  dv  and  realizing  that  the  residue 


E  )dv  by  taking 
V  “  "  V 


’  r 

"  r 

is 

J  (E  •  E  )  dv 
Lv  "  J 

J  (E  •  E  )  dv 
Lv  -  r*  J 

.  Then  E  may  be  found  by 
& 


comparing  the  coefficient 


f  ^ 

j  (E  •  E 


E  )  dv  with  E .  These  operations  give 
a. 


^0  /abd  jut 
(E  •  EJdv  =  j  — s /—  e 


5a 


I 

V 

nr  r  P  P 

O-  i  jX  cosBx  sinp  y  sinp  z  -  j  sinp  x  cos  p  y  sinp  z 

</  abd  I”  Pc  ^  z-P^  X  y 
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A  similar  operation  on  H  gives 


C.  MAXWELL'S  EQUATIONS  WITH  ORTHOCX)NAL  FUNCTIONS 


It  is  now  possible  to  substitute  the  series  obtained  in  section  A 

into  Maxwell's  equations  (assuming  pi  and  c  are  constant  throughout 

the  cavxty  volume).  Of  course,  the  resulting  equations  will  be  a 

summation  of  series  involving  various  combinations  of  the  basic 

functions  E  ,  H  .  G  ,  and  F  .  Actually  the  primary  interest 
—  a  —a  — \  —a 

of  these  equations  are  the  coefficients  of  these  functions,  as  expressed 
by  the  various  volume  and  surface  integrations  [see  for  example 
Equations  (63)  and  (64;.  A  particular  "set"  of  coefficients  may  be 
examined  by  using  the  orthogonal  properties  of  the  basic  functions. 

This  process  is  illustrated  in  the  next  few  paragraphs. 

From 
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Multiply  through  this  equation  by  H  ,  integrate  over  the  volume,  and 
obtain  an  equation  in  terms  of  the  coefficients  for  this  particular  H 
The  same  technique  can  be  used  with  so  that  for  the  coefficients, 


f  (E-  E  )  dv  +  (h  •  (n  X  E 


E)dA=  J  {H  HJdv 


(80) 


and  for  the  coefficients. 


r 


G  •  (n  X  E)  dA  =  -  ^  g|- 


/(H  .  G  ) 


dv 


(81) 


dD 


Similarly  2.x  H  -  =  J  becomes 


.  /(H  HJdv-.  -  /  (E  •  E J  dv  =  J" (J  .  E  )  dv  -  f  •  (n  x  H) 

a  Y  -  -a  dt  Y  -  Y  “  8 '  ~  ” 


€  i  f  (E  .  F Jdv  s  /(J. 


F  *  (n  x  H)  d  A  -  «  __  ^  -  *  » 

— o  —  —  dt  Y  “■* 


The  equation  v  *  D  =  p  gives 


F  )dv 

Y  -  “o 


dA  , 

(82) 

(83) 


f  (D  •  F  )dv  +  /"  (il»  D  •  n)i 


-k  J  (D  •  F  )dv  +  J  (i)(  D  •  n)dA  =  J  (pi);  )dv 
o  Y  “  s '  *  ”  ~  V 


(84) 


The  equation  £  •  B  =  0  gives 


k  f  (B  •  G Jdv  =  /(♦^  B  .  n) 
X  Y  ""  “"X 


dA 


(85) 


Now  from  Equations  (80)  to  (85)  the  expansion  coefficients 


(E  .  EJdv  ,  j  (H  •  H  )dv  ,  j  (E  •  F 


F^)  dv  ,  and  (H  •  G^)  dv 


-74- 


may  be  found  and  then  substituted  into  Equation  (66)  to  give  the  electric 
and  magnetic  fields  in  the  cavity.  Slater^  has  shown  that  Equations  (83) 
and  (84)  lead  to  the  continuity  of  electric  charge, 

Equations  (81)  and  (85)  may  be  combined  by  taking  the  time  derivative 
of  (85)  and  substituting  into  (81)  to  obtain 


-If/**,.? 


J  G  •  (n  X  E)  dA  =  J  V 


•  (n  X  E)  dA 


(87) 


n)  dA  = 

To  simplify  the  second  term  introduce  the  differential  operator  on 

the  surface,  which  is  equivalent  to  ^  -  n  ^  .  By  vector  identity, 

on 

^A  “  ®A  \  '  <11  ^  ‘  (nxE) 

By  integrating  this  over  the  surface  and  realizing  that  =  Z  4>x. 

a 

(since  -r—  A.  =  0  on  s  and  =  0  on  s'), 
on  \  X 


n  X  E)  dA  s:  •  (n  X  E)  dA  +  ^ 

3 1  8  8  8 


^  E  '  (n  X  E)  dA 
^  A 


The  operator  E.  treats  the  surface  as  an  equivalent  "volume"  so 


that 


f 


2  6  ■  (n  X  E)  dA  =  - 

sTs'  A  X 


j  E) 


dl 


where  dl  is  on  the  perimeter  of  the  cavity  surface.  However,  the 
line  integral  about  the  surface  surrounding  the  cavity  volume  is  zero 
since  its  perimeter  vanishes;  therefore 


/ 


■  (fl  *  E)  dA  ~  ~  J  •  (n  X  E)  dA 


and  Equation  (87)  becomes 
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or 


S  8 


n)  =  2^^  •  (n  X  E) 


(n  X  E)  dA 


(88) 


Equation  (88)  expresses  the  conservation  of  the  fictitious  surface  magnetic 
charge  just  as  Equation  (86)  expresses  the  conservation  of  electric  charge. 
This  leaves  Equations  (70)  and  (82),  which  may  be  combined  and 


/  (E  • 


/ 


solved  for  J  (E  •  E  )  dv  and  J  (H  •  H  )  dv  .  This  solution  is 
V  “  -a  V  -  -a 

straightforward  since  these  equations  are  of  the  form 


ax  +  by  =  c 
dx  +  cy  =  f 


Their  solution  gives 
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These  two  equations,  (89)  and  (90),  form  the  basis  for  the  rest 
of  this  development.  In  fact,  these  two  equations  will  usually  be  the 
most  important  in  any  type  of  analysis  using  this  orthogonal  expansion 


(£  •  E  )dv  and 
—  —a 


for  the  cavity  fields.  The  series  coefficients 
determine  the  solenoidal  parts  of  E  and  H ,  which  is  the  part  showing 
properties  of  wave  propagation.  This  part  of  the  field  is  usually  regarded 


as  the  "radiation"  field.  The  other  coefficients 


/(H-  G 


J  (E'  F  )dv  and 


G  )dv  which  determine  the  irrotational  parts  of  E  and  H  may  be 


found  from  Equations  (81),  (83),  (84),  and  (85).  This  part  of  E  and  H 
of  course  has  zero  curl,  and  hence  is  derivable  from  a  scalar  potential 
as  shown  ly  Equation  (51).  Except  for  the  a  =  0  and  X.  =  0  functions, 
the  problem  of  finding  these  parts  of  E  and  H  is  similar  to  the  electro¬ 
static  or  magnetostatic  problem  but  with  time -varying  charge  distributions. 
Note  that  no  account  has  been  taken  here  of  time  retardation  or  of  finite 
velocities  of  propagation  of  a  disturbance. 


D.  ILLUSTRATIVE  PROBLEM 


It  would  be  of  value  at  this  point  to  see  just  how  the  material 
derived  may  be  used.  Consider  a  cavity  with  ideal  metal  walls  (no  s' 
surfaces)  filled  with  a  lossy  dielectric.  Thus  in  the  cavity  v  ^  0, 
therefore  J  =  (tE  ;  also  (n  x  E^)  =  0,  since  the  walls  are  ideal.  With 
this  information  Equation  (89)  becomes 


d  ,  2 


j(E. 


Ea)dv 


=  0 


(91) 


Assume  a  solution  of  the  form 


where  w  is  the  complex  frequency; 


u  =  w 

1 


Now  Equation  (91)  may  be  solved  by  realizing  that 


(ideal  resonant  frequency)  , 


and 


(dielectric  Q) 


This  solution  yields 

2  2^2 
«a  =  Wj  +  <-2 

and 


so  that  the  oscillation  frequency  becomes 


E.  GENERAL  PERTURBATION  OF  CAVITIES 

The  general  case  of  perturbation  by  the  insertion  of  any  substance 
in  a  cavity  can  be  divided  into  these  three  classes: 

(1)  fields  do  not  penetrate; 

(2)  fields  totally  penetrate  the  substance;  and 

(3)  fields  partially  penetrate  the  substance. 

These  correspond  approximately  to  the  respective  conditions 

^  »  (  —  «  t  ,  and  —  S  t  .  Since  the  first  and  second  classes 

jw  ju  jw 

most  nearly  match  the  usual  physical  situation  when  a  substance  is 
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intentionally  inserted,  only  these  will  be  treated  here. 

The  general  method  of  attack  in  both  these  cases  will  be  to  assume 
a  given  ideal  cavity  in  which  the  fields  E  and  H_  are  known.  Then  if 
the  object  inserted  in  the  cavity  is  small  so  that  the  perturbation  is  also 
small,  the  original  fields  will  be  only  slightly  altered.  This  slight  change 
can  be  accounted  for  by  applying  a  correction  term  to  the  original  fields. 

In  this  manner  the  frequency  shift  caused  by  the  perturbing  object  can  be 
calculated. 

Perhaps  some  insight  may  be  gained  into  the  problem  of  correcting 
the  original  fields  by  considering  Poynting's  power -flow  theorem. 


(E  X  H)  dA 


dE 


E)  dv 


(92) 


where  the  left  side  is  the  closed  integral  on  the  area  surrounding  the 
volume  and  represents  the  power  flow  into  the  volume.  Suppose  only 
the  a^^  mode  is  present  so  that  the  following  may  be  substituted  for 
E  and  H : 


H 


H 

—a 


H  )dv 
—  a 


(93) 


where  only  the  solenoidal  functions  is  now  considered  since  this  is  the 
important  one  in  determining  the  "radiation"  field  mentioned  earlier. 
Equation  (92)  becomes 
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Equation  (92)  has  now  been  applied  to  an  ideal  cavity,  hence  cr  =  0  inside 

2  * 

the  cavity.  In  Equation  (94)  is  understood  to  mean  and 

E  ^  means  E.  ‘  E  The  same  is  true  for  the  square  of  the  series 
a  —a  —a  ^ 

coefficients, 


(E  •  E2)dv  and 


(H*  H^)dv 


These  two  coefficients  may  be  related  by  using  Equation  (80)  and  the 


fact  that 


for  an  ideal  cavity.  Equation  (80)  becomes 


k 

a 


f(E 


E  )dv 
—a 


-JliU 


L 


H  )dv 
—a' 


/(E 


EJdv 


(95) 


Of  course,  this  90°  time  relationship  between  the  electric  and 
magnetic  fields  is  just  what  is  expected  in  an  ideal  cavity.  Substituting 
Equation  (95)  into  Equation  (94)  gives 


* 
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(96) 


)  dv 


This  is  the  result  sought.  It  shows  the  form  of  the  expression  for  the 
energy  flow  into  and  out  of  a  given  volume.  If  one  approaches  the  problem 
of  perturbation  from  the  energy  viewpoint;  that  is,  by  considering  the 
change  in  energy  over  the  volume  of  the  perturbing  object  caused  by  the 
different  electrical  properties  of  this  object,  then  one  should  be  using 
terms  of  the  form  of  Equation  (96). 

In  the  derivations  following,  the  approach  suggested  above  is  used 
directly  for  the  first  part,  and  indirectly  for  the  second.  In  reality,  the 
small  perturbations,  on  the  whole,  affect  the  cavity  fields  only  slightly; 
however,  immediately  in  the  vicinity  of  the  perturbing  object  the  fields 
are  often  grossly  changed.  In  the  first  analysis  this  effect  on  the  "near 
fields"  is  neglected  and  is  included  later  in  a  term  called  the  calibration 
constant.  This  method  is  used  to  give  a  clear  physical  picture  at  first, 
before  giving  a  more  precise  final  result. 


1.  Class  1  Perturbation  ( —  »  €  ,  or  metal  beads) 

jw 

Assume  a  cavity  with  perfect  metal  walls  and  operating  in  the 
a^^  mode  for  which  the  fields  and  and  their  respective  co¬ 
efficients  are  known.  Now  suppose  a  small  metal  bead  of  volume  r 
is  placed  in  this  cavity  such  that  the  perturbation  is  small.  Figure  24 
shows  the  original  and  perturbed  states  with  the  associated  volumes 
and  surface  normals  involved. 

First  examine  Equation  (89)  with  the  current  integral  dropped 
(since  J  =  0)  and  the  surface  integral  over  s'  dropped  (no  s'  surfaces). 
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ORIGINAL  CAVITY  CAVITY  WITH  BEAD 


Figure  24.  Perturbed  Cavity. 


This  becomes: 

2 


+  k^  /(E  E)dv=-k  H  (nxE) 
a  J  -  -a  aJ-a-~ 


dA  . 


dt  V 

Since  the  cavity  is  assumed  ideal,  the  right  side  should  be  zero 
(no  tangential  component  of  E  at  the  metal  walls).  Permute  the 
order  of  the  vectors,  and  substitute  Equation  (93),  since  only 
the  a^^  mode  is  present. 


The  right  side  of  this  equation  is  now  recognized  from  Equation  (96). 

This  terms  acts  as  a  "driving  force"  and  represents  energy  being  fed 
into  the  cavity.  In  an  ideal  cavity,  as  was  just  pointed  out,  this  term 
is  zero.  If  the  cavity  has  wall  losses,  then  the  term  is  not  zero  (there 
will  now  exist  a  small  tangential  component  of  E  at  the  walls)  since 
energy  must  be  fed  into  the  cavity  if  it  is  to  continue  to  operate. 

Now  to  return  to  the  metal  bead  inside  the  cavity  shown  in 
Figure  24,  it  is  apparent  that  the  bead  removes  a  small  portion  of  the 
cavity  fields  or  of  the  internal  energy  since  no  fields  can  exist  inside 
the  metal  bead.  Of  course  this  alters  the  original  fields  slightly.  However, 


4 
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the  original  fields  may  still  be  considered  to  be  present  if  the  region  of 
the  bead  is  considered  as  an  energy  "sink.  "  That  is  in  Equation  (97) 
the  right  side  of  the  equation  can  be  adjusted  to  represent  energy  removed 
from  the  cavity  by  the  bead.  The  surface  integral  is  then  over  the  surface 
of  the  bead  and  the  algebraic  sign  is  changed  to  indicate  energy  fed  out 
of  the  cavity  with  respect  to  the  normal  n.  Equation  (97)  becomes 


In  terms  of  n*  , 


n  =  -  n' 


Substituting  this  and  Equation  (96)  into  the  above  equation  gives 


(E  -  E  )dv  s  0  .  (98) 

V  -  -a 


Now  the  integral  over  r  represents  the  energy  removed  by  the  bead. 
Assuming  a  solution  of  the  form 


/, 


(E .  E  )  dv  =  C  e 

V  “  -a 


jut 


gives 


or 


^  4  f 

h)  =  «  U  +  i  i 

^  L  T 


(H  ^  -  E,^)dv 
'a  a  ' 


=  0 


(H^^  -  E^^)  dv 


(99) 


Equation  (99)  now  relates  the  new  resonant  frequency  u  to  the  original 
resonant  frequency  .  If  6  is  defined  as 
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0)  _  OJ 


6  = 


0) 


(100) 


and  a  small  perturbation  is  assumed,  Equation  (99)  becomes 


(H  ^  -  E,^) 

,  'a  a  ' 

6  =  - - -  T 


(101) 


The  calibration  constants  mentioned  earlier  can  now  be  included  so  that 
Equation  (101)  becomes: 


6  =  (K,^h/-K,^e/)  1 


(102) 


where  K  refers  to  class  1  electric  field  perturbation  and  K,.  to 
le  1“ 

class  1  magnetic  field  perturbation.  These  quantities  are  calculated 
in  a  later  section. 


O' 

2.  Class  2  Perturbation  ■: —  «  c  or  dielectrics  and  ferrites) 

ju 

This  case  is  quite  different  from  the  one  just  calculated  since 
the  fields  in  the  cavity  will  also  exist  inside  this  type  of  bead.  Figure  24 
still  applies,  only  now  the  cavity  has  the  properties  and  while 
the  bead  has  the  properties  (x  and  c  .  Equation  (92)  suggests  the  use 
of  the  displacement  current  9D/9t  and  the  fictitious  magnetic  current 
9B/9t  in  accounting  for  the  bead  effects.  In  this  way  it  may  be 
possible  to  calculate  an  equivalent  p'  and  t  ’  of  the  combination 
of  bead  and  cavity.  With  the  bead  inserted  in  the  cavity,  the  total 
displacement  current  is  given  by 


over  V 


over  T 


(103) 
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Equation  (93)  can  be  substituted  for  E,  and  the  resulting  expression 
used  in  Maxwell's  7  x  H  equation  to  integrate  over  the  volumes  con¬ 
cerned  to  find  the  relative  contributions  of  displacement  current: 


Now  one  can  differentiate  with  respect  to  the  volume  V , 


E  J  (H  •  H  )  dv  =  c  ^  /(E  •  E  )dv  +  (e  -  c  ~  /  E  (E  •  E 

-a  V  -  -a'  0  dt  -  -a'  '  0'  dt  1  -  - 


E. ) dv  . 


V  ~  -  —  y  -  —  T 

Substituting  for  E  in  the  last  term,  dividing  through  by  E  ,  and  collecting 
terms  gives 


f  (H‘  HJdv  =  €, 


ay'-  -a' 


where 


1  +  (K^-l)  J  E,  dv 


^  J  (E.  E^)dv  ,  (104) 

dt  V  "" 


A  similar  solution  of  Maxwell's  2.^  ^  equation  can  be  obtained  using  the 
fictitious  magnetic  current  which  gives 


-a  ^ 


1  +  (K„.l)  j 

^  ^  1 
a 

T 

dt 


(105) 


where 


K_ 


m 


=  JL. 

^^0 
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Combining  Equations  (104)  and  (105)  to  solve  for 
/  (H  •  H  )dv  •< 


)dv  and 


-  -a' 


leparately  results  in 


Vo 


lt(K.-l)/E^<lv]  [l+K  -1)  /h 


m  ■  ^  a 


dv 


/(E.E  Jdvl 


2 


—  —a 


/(H-H, 


>  =0 


.--)dv 
V  -  y 


(106) 

In  comparing  this  equation  to  Equations  (89)  and  (90),  it  is  seen  that  the 
equivalent  p'  and  « '  are  as  follows: 


c  '  =  c 


t(K  .  1)  JeJ 


dv 


*  K, 


1  +  (K  -  1)  J  h/ 

m  ^ 


dv 


(107) 


As  is  customary  a  solution  of  form  e  is  assumed  for  Equation  (106); 
then 

2 

(108) 


.2  . 


w. 


«  =  r 


1  t  (K^  .  1)  /e/ dv  1  +  (K^  -  1)  / dv 


m  ’  a 

T 


Since  the  perturbation  is  assumed  small,  second  order  terms  can  be 
neglected  and  Equation  (100)  can  be  used  to  obtain 

E.^t  H 

6  =  -(K^-l)-^^ -  .  (109) 

Upon  insertion  of  the  calibration  constants  Equation  (109)  becomes 
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-  •'Zh  '•'m  -  *> 


(110) 


Ha^ 


where  refers  to  class  2  electric  field  perturbation  and  K2j^  to 

class  2  magnetic  field  perturbation. 


3.  Practical  Form  of  Perturbation  Equations 

If  one  is  unaccustomed  to  using  the  normalized  fields  and 
H.  .  then  Equations  (102)  and  (110)  may  be  written  in  terms  of  the 
conventional  E  and  H  fields  as  follows: 


6  = 


6  = 


(Kj^^  H  -K,,.„E)t 


4U 


(Class  1)  , 


-  *>  •  0  ’  •'Zh  '•'m  ■  •)  ‘‘o  ^ 


W 


(111) 

(Class  2)  , 

(112) 


where  U  equals  the  total  stored  energy  in  the  cavity. 


F.  CALIBRATION  CONSTANTS 

It  was  mentioned  earlier  that  the  reason  for  introducing  the 
calibration  constants  is  the  effect  the  bead  has  on  its  near  fields  in  the 
cavity,  even  though,  on  the  whole,  the  fields  are  only  slightly  altered. 
This  calibration  constant  may  be  calculated  if  we  impose  the  conditions 
that  the  bead  is  small  compared  to  a  wavelength  both  inside  and  outside 
the  bead  so  that  the  bead  in  the  cavity  is  in  a  region  of  nearly  uniform 
fields,  and  also  that  the  bead  is  kept  several  diameters  away  from  the 
cavity  boundaries  so  as  to  avoid  image  effects.  This  calculation  may 
then  be  made  by  considering  the  effects  of  placing  the  particular  bead 
in  question  in  a  uniform  electrostatic  or  magnetostatic  field.  Thus 
the  problem  may  be  reduced  to  approximately  the  static  solution  of 
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Laplace's  equation. 


As  it  turns  out  the  calibration  constant  depends  on  four  factors: 

1)  material  of  the  bead, 

2)  shape  of  the  bead, 

3)  orientation  of  the  bead  in  the  field, 

4)  type  of  field  the  bead  is  perturbing  (£  or  H). 

Once  this  "constant"  is  determined  for  a  particular  set  of  these  conditions 

it  will  remain  constant,  independent  of  field  strength,  provided  the  assumptions 

made  concerning  size  and  nearness  to  boundaries  are  adhered  to. 

Since  the  static  solution  to  Laplace's  equation  involves  matching 

boundary  conditions,  only  beads  whose  boundaries  may  be  expressed 

simply  are  considered  here.  This  includes  spheres  by  the  use  of  spherical 

co-ordinates  and  ellipsoids  of  revolution  by  the  use  of  elliptical  co-ordinates. 

The  general  ellipsoid  with  three  unequal  axes  is  not  considered  here  since 

elliptic  integrals  of  the  second  kind  result.  As  an  example  the  calibration 

constant  for  a  metal  sphere  perturbing  an  electric  field  is  calculated. 

The  insertion  of  this  bead  into  a  uniform  electric  field  may  be 

shown  easily  by  the  use  of  Legendre  polynomials  to  cause  a  field 

£  -  3E  cos  6  ,  E-  s  0  ,  E  =  0  to  exist  at  the  surface  of  the  metal 
r  u  w 

sphere.  Here  r,  9,  and  4^  are  the  spherical  co-ordinates.  When 

this  metal  bead  is  placed  in  the  cavity  it  acts  as  though  it  is  placed 

in  a  region  where  the  electric  field  has  the  same  value  as  that  at  the 

2 

surface  of  the  bead.  Thus  E  over  the  volume  of  the  bead  is  taken  as 

2  2  2 
E^  =  9Eq  cos  0 

and  K,  is  defined  as  follows: 
le 
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(113) 


I 


dv 


K. 


le 


/ 


2 

E  dv 
0 


Substituting  into  Equation  (113)  gives 
a  ir/Z 


[ 


2  2  2 

9  E.  cos  9  2ir  r  sin0  d^  d 
0  0  r 


K 


le 


X.  2  4  3 

E  IT  a 

0  J 


=  3 


(114) 


where  a  is  the  radius  of  the  metal  sphere. 

Figure  25  shows  the  nomenclature  used  for  the  ellipsoids,  while 
Tables  6  ,  7,  anc  8  show  the  calibration  constants  for  different 

conditions.  Figure  26  gives  a  pictorial  representation  for  metal  prolate 
spheroids  as  a  function  of  the  slenderness  ratio.  In  the  tables  E//  and 
H//  mean  that  E  and  H  are  parallel  to  the  semimajor  axis  while  EX 
and  HX  mean  that  E  and  H  are  prependicular  to  the  semimajor  axis. 


G.  SPECIFIC  EXPERIMENTS 
1.  Calibration  Constant 

A  small  number  of  experimental  measurements  were  made  using 
metal  spheres  to  check  the  calculated  calibration  constant  for  this  type 
of  bead  in  an  electric  field.  A  rectangular  cavity  was  made  from  L-band 
wave  guide  so  that  the  internal  fields  could  be  calculated  accurately. 

The  cavity  size  was  12.  4  cm  x  24.  8  cm  x  105.  5  cm  and  the  resonant 
frequency  was  approximately  624  megacycles.  A  simple  arrangement 
was  used  to  measure  the  frequency  shift  directly.  Two  oscillators  were 
used,  one  to  supply  the  cavity  and  the  other  as  a  standard.  The  output 
signals  of  these  oscillators  were  mixed  and  then  fed  to  a  receiver.  The 
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Prolate  Spheroid 
a  >  b  =  c 


Figure  25.  Nomenclature  Used  for  Ellipsoids. 

« 


Table  6.  Calibration  Constants  for  ^heres. 


Type  of  Sphere 

H-Field 

E -Field 

Metal 

3 

2 

3 

3 

Dielectric  (K  =1) 
rti 

No  perturbation 

K^  +  2 

_ 

Ferrite  (K,  =  1) 
d 

K  +2 
m 

No  perturbation 

* 
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Table  7.  Calibration  Constants  for  Ellipsoids  of  Revolution 


Type  of 
Bead 

Material 

El 

HI 

Prolate 

^heroids 

Metals 

2e^ 

2e^ 

e  -  T  ^ 

(1  -  2p")  e  .  £l  in 

Dielectrics 

and 

Ferrites 

2e^ 

2e^ 

Oblate 

Spheroids 

Metals 

3 

e 

e^ 

Table  8,  Calibration  Constants  for  Ellipsoids  of  Revolution 


Type  of 
Bead 

— 

Material 

E// 

H// 

Prolate 

Spheroids 

Metals 

3 

e 

■ 

3 

e 

BQ 

OBI 

e  -  ^  /n 

Dielectrics 

and 

Ferrites 

e^ 

e^ 

(Kj-UP^ 

1^ 

m 

Oblate 

Spheroids 

Metals 

2e^ 

1 

2e^ 

„2  ^  -1  P  ^  2  -6^  IT 

p  [tan  £  +  -p-E_  e  - 
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•Pigure  26. 
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cavity  oscillator  was  adjusted  for  resonance,  then  the  bead  was  inserted 
and  the  oscillator  readjusted  for  resonance.  The  change  in  the  receiver 
reading  was  the  frequency  shift  caused  by  the  bead.  The  receiver  had 
been  previously  calibrated  against  a  100  kc  crystal  frequency  standard. 
For  measurements  of  this  type  the  suspension  used  for  the  bead  is  always 
a  problem,  since  this  can  also  perturb  the  cavity.  Silk  and  nylon  threads 
were  finally  used.  The  beads  were  inserted  through  the  side  of  the  cavity 
so  that  the  volume  of  the  thread  was  a  least  an  order  or  two  of  magnitude 
smaller  than  the  volume  of  the  bead.  The  results  are  shown  in  Table  9. 


Table  9.  Calibration  Constant  Measurements  for 
Electric  Field  Perturbation  by  Metal  Spheres 


Bead 

Diameter 

Diameter  in 
Wavelengths 

Measured 

Kle 

Percent 
of  Error 

Suspension 

Thread 

0.  177  in. 

.  0094 

3.  2404 

8.0 

.  002  in.  silk 

0.  500  in. 

.  0266 

3.  0148 

0.  5 

.  003  in.  nylon 

1 . 00  in. 

.  0532 

5.2760 

75.  9 

.  005  in.  nylon 

Since  the  frequency  shift  for  the  smallest  bead  was  on  the  order  of  six 

kilocycles  and  the  order  of  error  in  setting  the  cavity  on  resonance  was 

±Z  kilocycles.,  this  measurement  is  not  too  meaningful.  Of  qourse  this 

data  is  inadequate  for  making  any  predictions  but  the  results  do  indicate 

that  the  largest  bead  is  too  large  compared  with  a  wavelength  for  the 

theoretical  value  of  the  calibration  constant  to  be  valid.  The  measured 

frequency  shift  of  this  bead  was  1.  75  megacycles. 

Other  measurements  of  this  type  have  been  described  in  the 

10,  11 

literature. 
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2.  Electric  Field  Meaiurements* 

The  perturbation  technique  has  been  used  to  measure  the  electric 
field  along  the  axis  of  a  gridless -gap  coaxial  type  cavity.  A  sensitive 
scheme  for  measuring  the  frequency  shift  was  used  that  consisted 
essentially  of  a  low-frequency  sweep  generator  which  modulated  the 
signal  generator.  The  cavity  was  used  as  a  transmission  cavity  so 
that  the  dynamic  resonance  curve  could  be  viewed  on  an  oscilloscope. 

A  calibrated  signal  generator  used  in  conjunction  with  the  sweep 
generator  provided  a  frequency  marker  which  was  placed  on  the  "a" 
axis  of  the  oscilloscope.  This  marker  was  moved  to  line  up  with  the 
resonance  curve  and  was  used  to  measure  the  frequency  shift  directly. 
The  accuracy  of  the  equipment  alone  was  estimated  as  being  less  than 
one  per  cent.  Over- all  accuracy  for  the  final  measurements  was 
estimated  as  being  a  few  per  cent. 


W.  E.  Blair,  graduate  student,  now  studying  at  Cornell,  set  up  the 
eqqipment  and  made  these  measurements.  A  detailed  account  of  this 
may  be  found  in  Research  Report  ££  377,  Cornell  University,  May 
1958,  p.  35. 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  design  of  the  general  doubly  re-entrant  gridless -gap  coaxial 
cavity  has  been  completely  specified.  This  type  of  cavity  can  now  be 
designed  in  such  a  manner  that  at  least  the  first  three  higher -order 
resonances  do  not  occur  at  harmonics  of  the  signal  frequency.  The 
accuracy  of  this  method  at  the  fundamental  resonance  is  within  several 
percent  and  at  the  higher-order  resonances  the  error  is  somewhat 
greater.  The  design  techniques  given  should  cover  most  applications, 
but  if  not,  one  need  only  use  the  formulas  given  to  set  the  design  criteria 
for  any  special  case.  The  general  formulas  developed  will  of  course 
also  apply  to  any  cavity  that  is  derived  from  the  type  considered  here; 
i.  e. ,  singly  re-entrant  cavities  and  gridded-gap  cavities.  A  by-product 
of  this  theory  is  the  ability  to  obtain  the  value  of 

fundamental  resonance  with  an  accuracy  that  appears  equal  to  and  perhaps 
even  better  than  existing  methods.  No  general  calculations  of  either 
shunt  resistance  (Rg}^)  or  cavity  are  attempted  here  because  of 
the  strong  dependence  of  these  separate  quantities  on  surface  con¬ 
ductivity  and  incalculable  losses.  Accuracies  better  than  twenty  to 
thirty  per  cent  would  be  difficult  to  realize. 

Nonsymmetrical  modes  have  been  considered,  but  only  on  a 
small  signal  basis;  i.  e.  ,  these  were  only  considered  to  be  weakly 
excited  (see  section  II-A).  However,  any  present-day  high-power 
klystron  cavity  is  bound  to  excite  these  modes  strongly  because  of 
the  large  nonsymmetrical  devices  for  coupling  out  the  power.  The 
relatively  small  coupling  loops  used  in  the  cavity  model  described  in 
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section  II-C  excited  the  nonsymmetrical  modes  strongly  enough  so 
that  these  resonances  were  easily  measured.  One  can  imagine,  then, 
what  must  occur  in  an  output  cavity  with  a  large  coupling  loop  or 
even  an  iris.  Any  meaningful  analysis  of  a  cavity  under  these  con¬ 
ditions  would  be  exceedingly  difficult  to  make.  One  answer  to  the 
analysis  difficulties  and  to  the  practical  difficulties  presented  by 
present-day  output  coupling  systems  would  be  to  consider  a  symmetrical 
output  system  for  the  cavity.  It  would  certainly  be  worth  while  both 
analytically  and  experimentally  to  study  such  a  possibility. 

The  perturbation  of  resonant  cavities  has  been  treated  for 
the  types  of  materials  one  is  likely  to  use.  In  addition,  methods  for 
obtaining  the  formulas  for  the  calculation  of  absolute  values  from 
perturbation  measurements  have  been  developed.  It  has  been  shown 
that  the  perturbing  object  or  bead  must  be  on  the  order  of  two  or 
three  hundredths  of  a  wavelength  in  size  for  these  formulas  to  be 
accurate.  For  larger  beads,  one  can  usually  obtain  only  relative 
values  of  the  fields.  In  many  cases  however,  this  is  all  that  is  required. 
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APPENDIX 


Table  A-1 


b  cot  Sr  (a-1) 

0  a  ' 

Absolute  Value  of  the  Short-Circuited  Radial  Line  Admittance 


Pr 
^  a 

a  =  1.  5 

Q  =  1. 8 

a  =  2.  2 

0.  10 

24.  6920 

17.  0042 

12. 6626 

0.  15 

-  - 

11.  2736 

8.  3753 

0.  20 

12. 2905 

8.  4640 

6.  2644 

0.25 

6. 7221 

4.  9668 

0.  30 

8.  1653 

5.  5994 

4.  1058 

0.  35 

4.  7601 

3.  4818 

0.  40 

6.  1104 

4.  1463 

3.  0042 

0.  50 

4.  8361 

3.  2616 

2. 3237 

0.  60 

4. 0021 

2.  6676 

1. 8567 

0.  70 

.... 

2.  2355 

1. 5091 

0.  80 

2.  9466 

1.9039 

1.2344 

0.90 

1.  6380 

1. 0068 

1.00 

2.2961 

1.  4188 

0. ?114 

1.20 

1. 8493 

1.0714 

0.  4774 

1.40 

1. 5185 

0. 8002 

0.  1784 

1.  60 

1. 2609 

0. 5816 

0.  1234 

1.  70 

0.  4715 

0.  2883 

1.  80 

1.0507 

0. 3734 

0.  4721 

1. 90 

0.  2787 

0. 6861 

2.00 

0.  8728 

0.  1856 

0. 9469 

2.  20 

0.  7183 

0.  0050 

1. 7706 

2.40 

0. 5801 

0.  1926 

4.  0304 

2.60 

0.  4537 

0.4056 

44I9. 00 

2.  80 

0. 3350 

0. 6581 

4.  3027 

3.00 

0.  2219 

0.  9804 

2. 1101 
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Table  A-Z*.  UNIVERSAL  CURVE 

Second  Order  Higher  Mode  Inductive  Susceptance 
(Absolute  Value  Tabulated) 


a  =  1.  5 

a  =  Z.  0 

a  =  3.  0 

•^.ZirR)*^' 

x=l.Z5 

x=1.50 

x=1.75 

x=l.Z5 

x=1.50 

x=1.75 

x=l.Z5  x=:1.50  x=1.75 

.  05 

Z8.  75 

Z3.  05 

15.40 

16.  3Z 

13.05 

8.  65 

9.  38  7.  35  4.  77 

.  10 

14.  30 

11.07 

7.  84 

7.90 

6.01 

4.  ZZ 

3.98  Z.86  1.83 

.  15 

9.6Z 

-  — 

5.  11 

- 

2. 07  —  - 

.ZO 

7.  ZZ 

5.40 

3.  84 

3.  7Z 

Z.  73 

1. 90 

1.08  0.47  0.09 

.  30 

4.  83 

3.  63 

Z.  53 

Z.  30 

1.  64 

1. 07 

<0  <0  <0 

.40 

3.61 

Z.  64 

1.85 

1.5Z 

1.00 

0.  6Z 

.  50 

Z.  85 

Z.  17 

1.44 

1.01 

0.66 

0.  33 

.60 

Z.  36 

1. 77 

1.  16 

0.63 

0.  36 

0.  11 

.  70 

1.  96 

1.41 

0.96 

0.  35 

0.09 

<  0 

.  80 

1.69 

1.  19 

0.  79 

0.09 

<  0 

.90 

1.45 

1.01 

0.66 

<  0 

1.  00 

1.Z5 

0.56 

1.  10 

1.07 

— 

l.ZO 

0.94 

0.64 

0.41 

1.  30 

0.  8Z 

1.40 

0.  74 

This  table  has  been  calculated  by  assuming  a  two -radian  gap  at 
Uq/c  =  0.  5  .  The  scale  may  be  changed  by  multiplying  the  values  by 


t  “o^1 

sin 

r-/ 

"0 

sin  (x) 

where  p  is  the  gap  length  in  radians.  The  table  is  approximate 
and  is  meant  to  be  used  only  as  a  guide. 

Zirra 

♦♦  Note:  ZirR  =  —r —  is  the  second  higher  mode  cutoff  . 


-98- 


99 


REFERENCES 


1.  R.  D.  Campbell,  "Radar  Interference  to  Microwave  Communication 
Services,"  Electrical  Engineering,  77,  (October  1958),  p.  916. 

2.  L.  Denichel,  "Net  Point  Calculations  for  Fields  in  Re-entrant 
Cavities,"  Master's  Thesis,  Stanford  University  (1948). 

3.  T.  Moreno,  Microwave  Transmission  Design  Data,  New  York: 
Doyer  Publications  Inc. ,  1^53. 

4.  C.  C.  Wang,  "Electromagnetic  Field  Inside  a  Cylinder  with  a  Gap,  " 
J.  Ap.  Physics,  16,  (June  1945),  p.  351. 

5.  N.  Marcuvitz,  Waveguide  Handbook,  New  York:  McGraw-Hill, 
1951,  p.  40. 

6.  J.  C.  Slater,  Microwave  Electronics,  New  York:  D.  Van  Nostrand 
Co. ,  Inc.  ,  1950. 


7.  H.  B.  G.  Casimer,  "On  the  Theory  of  EM  Waves  in  Resonant 
Cavities,"  Philips  Res.  Rep.  No.  6,  (1951). 

8.  P.  M.  Morse  and  H.  Feshback,  Methods  of  Theoretical  Physics, 
New  York:  McGraw-Hill,  1953,  p.  52. 

9.  K.  Kurokawa,  "The  Expansions  of  Electromagnetic  Fields  in 
Cavities,"  Proc.  1.  R.  E.  PGMTT,  6,  (April  1958),  p.  178. 

10.  L.  C.  Maire,  Jr.  and  J.  C.  Slater,  "Field  Strength  Measurements 
in  Resonant  Cavities,  "  J.  Ap.  Physics,  23,  (Jan.  1952),  p.  68. 

11.  S.  W.  Kitchen  and  A.  D.  Schelberg,  "Resonant  Cavity  Field 
Measurements,  "  J.  Ap.  Physics,  26,  (May  1955),  p.  618. 


-100- 


NOTE  ON  THE  PRELIMINARY  RESULTS 
OF  SPURIOUS  OUTPUT  POWER  MEASUREMENTS 


A.  Ray  Howland,  Jr. 


NOTE  ON  THE  PRELIMINARY  RESULTS  OF  SPURIOUS 
OUTPUT  POWER  MEASUREMENTS 


A.  Ray  Howland,  Jr. 

The  output  spectrum  of  a  large  klystron  was  measured  up  to  the 
fifth  harmonic.  Specific  measurements  were  made  on  the  fundamental, 
second  harmonic,  third  harmonic,  fourth  harmonic,  and  fifth  harmonic. 
The  fundamental  power  output  of  the  klystron  was  measured  in  the  output 
wave  guide.  The  harmonic  measurements  were  made  in  the  output  cavity 
of  the  klystron  by  means  of  a  small  coupling  loop  (sampling  loop)  placed 
in  the  output  cavity.  The  sampling  loop  was  not  calibrated  so  that  all  the 
harmonic  power  measurements  show  the  relative  change  in  a  particular 
harmonic  as  a  function  of  r-f  drive  power  and  as  a  function  of  the  beam 
voltage.  All  of  the  spurious  output  measurements  made  to  date  have  used 
the  SAL- 36  klystron. 

Two  earlier  reports^’  ^  have  described  the  measurement  technique 
and  the  reasons  for  using  the  sampling  loop  rather  than  attempting  to  maki 
model  power  measurements  in  the  output  wave  guide  of  the  tube  under  test 

The  measurement  circuit  used  to  measure  the  fundamental  output 
power  as  a  function  of  r-f  drive  power  and  as  a  function  of  beam  voltage 
is  shown  in  Figure  1.  The  measurement  circuit  used  to  measure  the 
second,  third,  fourth,  and  fifth  harmonic  outputs  is  shown  in  Figures  1 
and  2.  Briefly  the  fundamental  power  was  measured  by  means  of  a  pair 
of  cross -guide  couplers  in  the  output  wave  guide.  The  harmonic  measure¬ 
ments  were  made  by  sampling  the  harmonic  power  of  the  output  cavity  by 
means  of  the  small  coupling  loop.  The  sampling  loop  output  was  fed  throu( 
the  appropriate  filters  as  shown  in  Figure  2.  The  harmonic  output  was 
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then  fed  to  a  mixer  where  it  was  combined  with  the  local-oscillator  signal 
to  produce  a  30-Mc/s  signal  that  was  first  amplified  by  an  i-f  amplifier 
and  then  further  amplified  and  detected  by  a  receiver  for  display  on  an 
oscilloscope.  In  making  the  measurements  the  vertical  displacement  of 
the  signal  on  the  oscilloscope  was  kept  constant.  The  relative  change  of 
a  signal  was  measured  by  a  calibrated  attenuator  in  the  harmonic  signal 
output  line  and  by  a  calibrated  attenuator  in  the  receiver.  The  relative 
change  of  a  harmonic  signal  in  decibels  is  the  amount  that  the  attenuators 
had  to  be  changed  to  keep  the  display  on  the  oscilloscope  constant.  The 
r-f  input  drive  power  was  measured  by  means  of  a  pair  of  coaxial  direc¬ 
tional  couplers. 

In  making  the  harmonic  power  measurements,  it  was  necessary  to 
use  a  high-pass  filter  with  a  large  attenuation  in  the  stop  band.  The  only 
type  of  high-pass  filter  presently  available  that  is  capable  of  removing 
the  fundamental  and  all  harmonics  of  lower  frequency  than  the  harmonic 
being  measured  sufficiently  well  is  a  wave  guide  operated  beyond  cutoff. 
The  high-pass  filters  for  the  second,  fourth,  and  fifth  harmonics  were 
constructed  for  these  measurements.  The  third  harmonic  filter  was  made 
from  standard  S-band  wave  guide. 

Measurements  of  the  relative  change  in  second  through  fifth  har¬ 
monics  of  the  SAIj-36  are  shown  in  Figures  3  and  4.  Figure  3  shows  the 
changes  in  the  harmonics  for  a  beam  voltage  of  approximately  65  kv; 
Figure  4  shows  the  changes  for  a  beam  voltage  of  approximately  80  kv. 
The  nominal  operating  beam  voltage  of  the  SAL-36  is  120.  It  is  seen 
from  both  Figures  3  and  4  that  in  the  small-signal  case  the  rate  of  change 
in  the  fundamental  or  in  the  harmonics  with  respect  to  r-f  drive  power 
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3 

is  at  the  rate  predicted  by  small -signal  theory.  The  drive  level  for  which 

this  small-signal  analysis  is  valid  is  a  function  of  the  beam  voltage  of  the 
tube  under  test.  This  is  demonstrated  by  making  a  comparison  of  Figures 
3  and  4,  where  the  beam  voltage  is  the  only  parameter  changed  for  the 
two  sets  of  data.  It  is  noted  that  as  the  beam  voltage  is  increased,  the 
harmonics  saturate  at  a  lower  drive  level. 
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INPUT  POVVER  (NORMALIZED  FOR  MAXIMUM  VALUES) 


Figure  3.  Variations  of  Fundamental  and  Harmonic  Powers  as  a 
Function  of  Input  Power  for  a  Beam  Voltage  of  65  kv. 


Figure  4.  Variations  of  Fundamental  and  Harmonic  Powers  as  a 
Function  of  Input  Power  for  a  Beam  Voltage  of  80  kv. 
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ABSTRACT 


The  perturbation  technique  described  is  a  method  by  which  the 
electric  and  magnetic  field  intensities  can  be  measured  at  any  point 
inside  a  resonant  microwave  structure.  Small  magnetic  or  dielectric 
objects,  generally  spheres,  are  used  to  distort  the  electric  or  xnag- 
netic  fields.  In  the  electric  case,  this  field  distortion  causes  a  slight 
decrease  in  the  resonant  frequency  of  the  cavity,  sfhich  is  propor¬ 
tional  to  the  field  strength.  Data  can  be  obtained  very  precisely  to 
within  0.  3  per  cent  and  the  accuracy  within  1  per  cent.  A  finite  bead 
causes  slight  errors  in  resolution,  but  a  decrease  in  bead  size  causes 
too  great  a  loss  in  sensitivity. 


A.  DESCRIPTION 


This  technique  provides  a  method  of  accurately  measuring 
electri<^  and  magnetic  field  intensities  at  any  position  in  a  resonant 
microwave  structure.  The  method  by  which  this  is  done,  is  to  use 
a  perturbing  object  which  disturbs  the  natural  field  pattern  of  the 
structure  and  changes  its  resonant  frequency. 

The  perturbing  object  can  be  made  of  any  substance  and  have 
any  geometry;  but  a  small,  spherical  bead,  made  of  metal  or 
dielectric,  is  the  simplest  to  analyze.  The  apparatus  is  somewhat 
complex  to  set  up,  but  once  initial  calibrating  and  balancing  pro¬ 
cedures  are  complete*  data  can  be  taken  quickly  and  accurately  with 
less  than  one  per  cent  er  or. 


B.  PERTURBATION  THEORY  AS  APPLIED  TO  THIS  TECHNIQUE 
The  general  perturbation  equations^  are  used  as  the  basis  of 
the  theory;  for  metals, 


K„  p  K  €  E^ 

,  df  Ih  o  -  le  o 

6  =  j-  =  - TO -  ’ 

a 


(1) 


and  for  dielectrics  and  ferrites. 


6  = 


df 


-K,^-  (k  -l)p  H^  -  K,  fK  -  l\  E' 
2h  \  m  _ 2e  V  d  ./  o  . 


■W 


T  ;(2-r 


where 


E  =  electric  field  intensity; 


^L.  A.  MacKenzie ,  tCly stron  Cavities  for  Minimum  Spurious  Output 
Power"  Research  Report  EE  418,  Equations  (111)  and  (112). 
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H  =  magnetic  field  intensity; 

K  =  relative  permeability  for  dielectric; 

m  o 

=  dielectric  constant  </c  for  dielectric, 
d  o 

T  =  volume  of  the  perturbing  bead; 

U  =  total  stored  energy  in  cavity; 

correction  constants  that  are  functions  of  bead  shape, 

=  size,  and  material;  of  type  of  field;  and  to  some  extent 
of  the  position  of  bead  in  the  cavity: 


df  =  incremental  change  in  the  resonant  frequency  caused  by 
bead,  frequency  deviation,  and 
f^  =  natural  resonant  frequency  in  a^  mode. 

Because  of  the  simplicity  of  analysis  and  results,  metallic  and 
dielectric  spherical  beads  will  be  the  only  types  considered  in  calcu¬ 
lating  the  correction  constants. 

Consequently  if  the  E  and  H  field  strengths  are  desired  at 
a  particular  point  in  a  microwave  cavity,  both  a  metallic  bead  and  a 
dielectric  bead  can  be  placed  at  this  point  and  the  frequency  deviation, 
df,  can  be  measured  for  both  beads.  Substituting  the  values  of  df  into 
Equation  (1)  and  (2),  and  knowing  the  other  variables,  E  and  H  can  be 
obtained  by  simultaneous  solutions  of  these  equations. 

For  convenience  of  experimentation,  these  equations  may  be 
simplified  by  making  the  proper  choice  oft 

(1)  the  geometry  of  the  microwave  structure, 

(2)  the  position  at  which  the  bead  is  placed. 
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In  a  cylindrical  microwave  cavity  of  a  general  geometry  of  Figure  1, 
aisume  that  the  bead  is  moved  along  the  z-axis  of  cavity.  This 
simply  means  that  with  small  enough  beads,  the  H  field  and  the  radial 


E  field  are  theoretically  zero.  Then  E  becoxnes  E  ,  the  electric 

cza 

field  along  the  z-axis; 


K,  «  E  2 
le  o  cza 

4U 


(3) 


and 

....  K-  (k,-l)*  e2 

r_df  2e\d/o  cza 

6=7-=- - TO - •  < 

a 

For  Equation(2}the  H-term  is  also  zero  because  K  for  most  dielec- 

m 

tries  is  approximately  one. 

Equations  (3)  and  (4)  can  be  written  in  the  form 


Figure  1.  Cylindrical  Microwave  Klystron  Cavity. 
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or 


This  form  is  maintained  for  a  more  general  presentation  of  results,  as 
shown  later.  Essentially,  varies  as  VdT  for  either  a  metal  or 

dielectric  bead. 

For  cavities  of  this  geometry,  obtained  from  an 

electrolytic  tank.  Its  general  form  is  shown  in  Figure  2. 


Figure  2.  E^^a  ^*'*'**  ^  Klystron  Cylindrical  Cavity 

Let  the  maximiim  value  of  E  =  E '  so  E '  oe  df  maximum.  The 

cza 

vertical  axis  can  be  normalized  to  E'  and  df  maximum.  This  allows 
Equation  (6)to  become 

^cJ'^  _  5cj;  .  IS - 

=77^  E' 
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The  horizontal  axis  can  be  normalized  by  d  ,  the  separation  of  the  ends 
of  the  drift  spaces.  This  normalized  is  shown  in  Figure  3. 


Figure  3.  E  ,  Normalized  by  E',  versus  z/d 

C  Za 

The  most  general  plot  E  can  be  obtained  by  normalizing 

cza 

E  by  Egjjgi  ®czo  uniform  equivalent  field  of  two  parallel  plates 

as  in  the  equivalent  gridded  klystron  gap  with  a  voltage  V,  separated 
by  the  distance  d  ,  such  that  E^^^  =  V/d  .  However,  for  normalized 
E^^^  in  Figure  3,V  =  area  under  the  curve  of  Figure  2  =  ^  E  •  dl. 
Since  E  has  a  tendency  to  spread  out,  it  is  logical  to  asstime  E'  to 

C  ZSl 

be  less  than  E,.__, 

E  E 

CZA  C  Zdr 

Letting  the  ratio  - -  =  « — ^  y  -^and  normalizing  Equation  6 

'^'czo  ■'^czo/'^U 

by  E  Equation  8  is  obtained: 

'  czo  ’ 

^cza/*^  _  E^ /  df 

E  /«fu  E  A/df  ' 

czo/  ^  ^  czo  V  max 

Plotting  Equation (8)  with  the  gridded  gap  produces  a  curve  similar  to 
that  in  Figure  4,  where  the  two  areas  are  equal. 
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Fisure  4.  E  Normalized  by  E  versus  z  Normalized  by  d 
®  cza  czo 


The  equatdoas  are  oi^en  written  in  the  form  of  E  /•/!] 

cza 

rather  than  merely  E  *  This  is  used  because  the  value  of  U  tn 
''  cza 

Equation  3  and  4  is  often  not  known.  Since  all  data  are  normalized, 
it  is  not  needed  unless  the  values  for  E^^^  in  volts/m&ter  aro.  desired. 
This  can  be  obtained  by  knowing  U  and  merely  scaling  Figure 
U  can  be  obtained  in  various  ways,  for  example,  by  measuring 
power  input  and  Q  from  the  basic  definition  cf  Q: 

w  (stored  energy)  u  U 

Q  -  — - 5 - .  (9) 

(power  input)  W 

C.  METHOD  OF  PERTURBATION  MEASUREMENT 

Figure  5  is  a  general  block  diagram  of  the  circuit  for  making 
perturbation  measurements.  The  nucleus  of  the  setup  is  the  balanced 
modulator.  It  mixes  a  main  signal  generator  output  f^  with  a 
Marka  Sweep  output  Af .  The  main  signal  f^  is  for  example  1  Mc/s 
above  the  natural  resonance  of  the  cavity  f^  .  The  Marka  Sweep 
generates  a  constant  voltage,  but  a  varying  frequency  output  that 
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Figure  5.  Diagram  for  Perturbation  Measurement 


modulates  the  carrier  .  This  Af  varies  in  frequency  linearly  and 
periodically  in  time  from  50  kc  to  5  Mc/s;  that  is,  it  is  a  sawtooth- 
varying  frequency  with  sweep  rate  of  60  cps. 

The  balanced  modulator,  which  suppresses  f^  ,  passes  the 
sidebands  f ^  ±  Af  into  the  cavity.  The  lower  sideband  will  pass 
through  the  cavity  resonance  since  f  >  f  .  Consequently,  energy 
of  the  upper  sideband  and  the  lower 'Sideband  (ether  than  at  resonance) 
will  be  reflected:  so  the  cavity  output  will  be  s  superposition  of  the 
two  sidebands.  (Figure  6a}  The  sharpness  of  the  dip  depends  on 
the  Q  of  the  cavity. 
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Figure  6.  Cavity  Output  versus  Frequency 
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The  output  for  high-powered  signals  may  be  taken  off  at  the 
directional  coupler;  low-powered  signals  may  be  coupled  directly 
out  of  the  cavity,  with  little  attenuation.  This  cavity  output  is 
detected,  differentiated,  and  finally  becomes  the  Y-axis  input  to 
the  Tektronix  oscilloscope  (Figure  6b).  The  wave  form  is  differ¬ 
entiated  because  the  interaction  of  the  curve  and  the  reference  line 
is  more  easily  distinguishable  in  Figure  6b  than  in  Figure  6a. 

In  order  that  the  wave  form  of  Figure  6b  be  viewed  by  an 
oscilloscope,  it  is  necessary  that  the  X-axis  be  swept  with  a 
60  cps  sawtooth.  In  addition,  it  is  necessary  that  the  X-axis 
voltage -varying  sawtooth  be  exactly  proportional  to  the  frequency- 
varying  sawtooth.  The  Marka  Sweep  provides  such  a  sweep; 
this  is  the  X-axis  input. 

Finally,  a  calibrated  signal  generator  is  used  as  an  external 
source  to  trigger  a  marker  pulse  in  the  Marka  Sweep.  This  pulse 
is  triggered  at  the  exact  instant  the  Marka  Sweep  modulating  output 
passes  through  the  frequency  at  which  the  calibrated  signal  generator 
is  set.  This  pulse  is  used  to  Z-axis  (intensity)  modulate  the 
oscilloscope.  Thus  a  pip  is  produced  at  the  instant  the  output  signal 
(Y-axis)  passes  through  frequency  •  where  f  is  the  frequency 

of  the  calibrated  signal  generator. 

With  f^  >  f ^  by  1  Mc/s,  f  can  be  set  at  1  Mc/s.  This 
would  position  the  pip  exactly  at  the  cavity  resonance  (Figure  7a). 

If  the  perturbing  bead  is  positioned  in  the  cavity  at  some  point 
z*  ,  on  the  z-axis,  the  resonant  frequency  will  be  reduced  by  a 
small  amoxmt,  df' .  This  causes  the  resonant  curve  to  shift  to 
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the  right  by  df'  (Figure  7b).  Then  by  decreasing  the  frequency  of 
the  calibrate^) signal  generator,  the  pip  can  be  adjusted  back  to  the 
zero  point  of  the  wave  form  (Figure  7c).  The  amount  the  frequency  . 
changed  is  recorded  as  df'  at  the  point  z*  . 

Next,  move  the  bead  to  the  next  point,  z"  ,  and  find  the  fre¬ 
quency  change  df"  .  This  procedure  is  continued  until  point  by  point, 
data  has  been  obtained  for  frequency  deviation  df  versus  bead 
position  z  .  This  data  is  used  to  obtain  the  normalized  plots  from 
Equation  (8)  and  Figure  4. 

D.  EXPERIMENTAL  RESULTS  AND  CONCLUSIONS 

Figure  8  shows  the  experimental  setup  as  it  was  used  to  obtain 
data  of  df  versus  z  for  several  beads.  Beads  of  various  geometries, 
sizes,  and  materials  were  used  to  perturb  the  cavity;  they  were  sus¬ 
pended  in  the  cavity  by  nylon  thread.  The  cavity,  middle  cavity  of 
SAL- 36  klystron,  and  the  bead  support  arrangement  can  be  seen  at  the 
right  side  of  the  figure. 

Some  beads  were  used  to  verify  the  theoretical  values  for  cor¬ 
rection  constants  (K^^  ,  data  for  one,  a 

spherical  brass  bead,  was  plotted  in  Figure  9:  this  can  be  compared 
to  Figure  4.  By  any  of  the  standard  methods  of  approximating  a 
curve  by  a  polynomial  (e.g. ,  the  least-square  method).  Figure  9  can 

be  represented  by  a  polynomial  in  z  and  since  the  bead  was  positioned 

•2 

along  the  z-axis,  the  polynomial  will  be  A(e)  *  That  is,  pi^(r,z)|^^Q  =  A(z). 
2 

L.  A.  MacKenzie, 'Electrolytic  Tank," Research  Report  EE  377, 

Equation  5. 
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HCLATIVE  AMPLITUDE  RELATIVE  AMPLITUDE  RELATIVE  AMPLITUDE 


Figure  7.  Three  Steps  in  Obtaining  Frequency  Shift 
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E  * 

The  experimental  value  of  im -  =  0.  758  and  the  general 

^czo 

shape  of  Figure  2  compares  well  with  that  obtained  in  the  electrolytic 
tank.  ^ 

Although  data  was  obtained  precisely,  with  error  less  than 
0.  3  per  cent,  there  are  possible  sources  of  error.  The  largest 
source  is  probably  resolution.  In  order  to  obtain  perfect  resolution, 
the  bead  size  must  be  small  enough  to  insure  a  imiform  field  across 
its  cross-section  and  must  distort  the  field  only  at  the  "point"  of 
measurement.  This,  however,  decreases  the  sensitivity  so  only 
strong  fields  can  be  measured.  Consequently,  a  compromise 
between  sensitivity  and  resolution  must  be  made. 


3 

Op.  cit.  ,  Figure  14. 
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